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Abstract 

This report deals with the quantum field theory of particle oscillations in vacuum. We first 
review the various controversies regarding quantum-mechanical derivations of the oscillation 
formula, as well as the different field-theoretical approaches proposed to settle them. We then 
clear up the contradictions between the existing field-theoretical treatments by a thorough 
study of the external wave packet model. In particular, we show that the latter includes 
stationary models as a subcase. In addition, we compute explicitly decoherence terms, which 
destroy interferences, in order to prove that the coherence length can be increased without 
bound by more accurate energy measurements. We show that decoherence originates not only 
in the width and in the separation of wave packets, but also in their spreading through space- 
time. In this review, we neither assume the relativistic limit nor the stability of oscillating 
particles, so that the oscillation formula derived with field-theoretical methods can be applied 
not only to neutrinos but also to neutral K and B mesons. Finally, we discuss oscillations 
of correlated particles in the same framework. 
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1 Introduction 



Nearly fifty years ago, the 6 puzzle led Gell-Mann and Pais |134] to discover that interaction (or 
flavor) eigenstates need not coincide with particles having well-defined masses and lifetimes. The 
concept of flavor-mixing was born. Pais and Piccioni |21£] went one step further and suggested 
that the propagation of a flavor eigenstate, as a coherent superposition of mass eigenstates, leads 
to the partial conversion of the initial flavor into other flavors. The flavor-mixing transition 
probability oscillates in space with a wavelength depending on the mass differences. The theory 
of flavor-mixing and particle oscillations was soon vindicated in a series of experiments on the 
kaon system | 190 , 114 , |148| , 171 1. Similar oscillations were observed in 1987 with B mesons 

By analogy, Pontecorvo proposed in 1957 the hypothesis of neutrino-antineutrino oscilla- 
tions p9|. 



^30| ]. Neutrino flavor-mixing, strictly speaking, was suggested somewhat later p03| ], 
and the two-flavor oscillation well as its application to solar neutrinos, was examined at 

Neutrino oscillations are, however, 



the end of the sixties by Pontecorvo and Gribov [ 231 , |14(: 
much more difficult to observe than meson oscillations, because of the small masses and the 
weak interactions of the neutrinos. They remained for thirty years a theoretical subject (see for 
example the reviews |66, ^), with the noteworthy exception of the solar neutrino experiment in 
the Homestake Mine |p8| , p7| , p^ , 34]. In the 90's, the persistent deficit in solar neutrinos (in com- 

observed at Homestake was confirmed 
|], GALLEX [||, |157|, GNO |l|] and 



parison with the predictions of solar models 40, 



by the experiments Kamiokande-II [161, |122|] , SAGE 



Super-Kamiokande 123 , |126[| . Moreover the Sudbury Neutrino Observatory recently announced 
the detection of a nonelectronic component in the solar neutrino flux j^, ^] , the total flux being 
compatible with the predictions of solar models. All these observations can be explained by 
neutrino mixing and oscillations p7|, Esl I 



127, 211] On another length scale, the observation 



of an angular dependence of the atmospheric muon-neutrino flux by the Super-Kamiokande Col- 



laboration |124| , |125|, |175| , |170[ gave a huge impetus to the neutrino oscillation hypothesis (note 
that a global muon-neutrino depletion had already been observed in other atmospheric neutrino 



experiments [16C, 121, 51, |l^, 15]). Preliminary results from the long-baseline experiment K2K 
]11, 164 ] confirm these results by showing a depletion in an artificially produced muon-neutrino 
beam. Finally, controversial evidence for neutrino oscillations has appeared in the laboratory 

29] but has not been confirmed (nor totally excluded) by the KARMEN 



experiment LSND [ 
experiment 27, 92]. 



Particle oscillations are an interference phenomenon which allows the measurement of ex- 
tremely small quantities, such as masses too small to be measurable by other means (mass 
differences between conjugate mesons ]114, 148, 13] or neutrinos |113| , |18C| ]), or CP violation 
(measured nearly forty years ago in neutral meson systems ]^, ^] and yet to be measured in 
the leptonic sector [103]). Thus, it is a privileged tool in the discovery of new physics. Nowa- 
days, B meson and neutrino oscillations have a huge importance in the determination of the 
unknown parameters of the Standard Model and its minimal extension, in particular the mixing 
matrix elements (including the CP violation phase) and the neutrino masses. Besides serving as 



4 



a delicate probe in particle physics, particle oscillations are a good testing ground for quantum 
mechanics. Einstein- Podolsky- Rosen correlations [ 162| , 100, IOC] have been studied theoretically 
|57, 58, 116, |22^ as well as experimentally [^] in the neutral K and B systems. Bell inequalities 



for neutral kaons are the subject of active research |22, 135, 76 1. Furthermore, fundamental 



issues of quantum mechanics are raised by oscillation experiments. The determination of the 
conditions under which the coherence between the different mass eigenstates is maintained is 



crucial for the phenomenology of neutrino oscillations |182, 204]. The question of the correlation 
between neutral mesons, in processes where B^B^ is included in the final state, is important for 



an estimate of the oscillation frequency |153]. 

Considering the practical importance of the oscillations, the theoretical framework describ- 
ing them should be rock-solid. Surprisingly, the validity of the standard treatment of oscillations, 
resorting to plane waves in a simple quantum-mechanical model, was not seriously questioned 
until neutrino oscillations finally gained credibility in the 90's. Two pioneering articles should 
be noted. First, a quantum field model of kaon oscillations was proposed by Sachs in 1963 p3§[| , 
but was soon forgotten, until the use of Sachs' reciprocal basis was recently revived for a rigorous 
quantum-mechanical treatment of meson oscillations 0]. Second, Kayser showed in 1981 
that neutrino oscillations are in contradiction with the plane wave assumption of the standard 
treatment of oscillations, and he proposed to modelize the oscillating particles with wave packets 



1 176]. His model, however, was not widely accepted and numerous articles have continued to 
appear to this day, discussing the ambiguities of the plane wave standard treatment, such as the 



proper choice of the frame ]201] or the equality of the energies or momenta of the different mass 



eigenstates ] 136 |. 

Although Kayser's wave packet model is a serious improvement on the standard plane wave 



treatment, it also suffers from serious defects: flavor states are ill-defined ]142]; the nonrelativistic 
limit is problematic if the mass eigenstates are not nearly degenerate ]142]; the question about 



the equality (or not) of the energy of the different mass eigenstates is not answered 1155, Ui 



the possibility remains of obtaining nonstandard oscillation formulas ]245]; the relation between 
the coherence length, beyond which oscillations vanish, and the conditions of production and 
detection is vague (problem of the determination of the wave packet size |182| ]); finally, it is not 
adequate for the oscillations of unstable particles. 

In order to solve these problems, a quantum field treatment of neutrino oscillations, similar 



to Sachs's model, has been proposed by Giunti, Kim, Lee and Lee [145]. The model consists in 
treating the oscillating particle as an internal line of a Feynman diagram, propagating between 
a source and a detector, which are represented by ingoing and outgoing external wave packets. 
This approach has been followed by other authors, but no agreement has yet been reached on 
the correct assumptions and calculation methods: 

1. There is disagreement on whether the interfering mass eigenstates have the same energy 
|15C] or whether their energy is spread; the existence of propagating wave packets, associ- 



ated to the mass eigenstates, is also at stake J 151] 
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2. Whereas oscillations vanish if the 3- momentum is measured too accurately, it is not clear 
whether the same thing happens when the energy is perfectly measured |183, 144, 139 1. 



3. Nonstandard formulas for neutrino oscillations have been derived in specific cases, either 
showing a strong dependence of the flavor-mixing transition probability on the neutrino 



mass ratios [242], or having a 'plane wave' behavior [163], or leading to oscillations of 
charged leptons in pion decay ]M^ 



Moreover, this formalism cannot be applied, as it is, to meson oscillations or to unstable 
neutrinos. The existence of neutrino mass makes indeed not only neutrino mixing possi- 
ble but also neutrino decay p 



75, 184, 192]. Although it cannot alone explain the data, 



neutrino decay has a long history as a possible explanation for the solar neutrino anomaly 
]3|, m, m, H, H, I, I, |0|, 1^, P, nil, H, l^]. Moreover this hypothesis was compatible 
until recently with atmospheric neutrino data | 



195 , 117 , 45, 202], although it now seems 



ruled out ]18C| ]. Neutrino decay has also been discussed with respect to supernovae [119, 193]. 

In this report, we intend to reconcile the different field-theoretical treatments in order to 
solve the above problems. Our first aim is to argue in favor of a field theory approach of 
oscillations. We want to show that it provides, on the one hand, the best available derivation of 
the oscillation formula and gives, on the other hand, a good explanation of the characteristics of 
the oscillation process (origin of the oscillations, observability conditions, boundary conditions). 
Our second aim is to strengthen the phenomenological basis of the neutrino oscillation formula 
through its unification with the mesonic oscillation formula. We thus extend the formalism to 
the unstable case, so that the same formula can be applied to K and B mesons, to neutrinos 
and to unstable neutrinos. 

In section ^ we discuss the quantum-mechanical approaches, so as to understand the prob- 
lems addressed in most of the articles and the different solutions adopted in the literature. The 
existing field-theoretical treatments are reviewed in section |3|, and their contradictions are sum- 
marized in five questions, which will serve as a guideline in the calculations of sections ^ to ^; 
these questions will be answered in section 7^ and |l^. In sections |^ and ^, we explain all the 
assumptions behind the external wave packet model. The relativistic case and the nonrelativistic 
(but nearly mass degenerate) cases are treated in a unified way, and the oscillating particle is 
allowed to be unstable. 

Oscillations of stable particles are analyzed in section |^. With the aim of reconciling the 
different points of view in the literature, we discuss different methods of computation, which 
allow to piece together a complete physical interpretation in section |^ We shall see that this 
formalism lends itself to a wave packet interpretation of the oscillating particle and can thus 
take up the physical insights of Kayser's intermediate wave packet model. The dispersion (or 
spreading in space-time of the amplitude) is taken into account, as it is important for mesons 
and might be relevant for solar and supernova neutrinos. We can then show that the existence 
of the coherence length is due not only to the separation of wave packets, but also to their 
dispersion. We compute the explicit dependence of the coherence length on the conditions of 
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production and detection, to check that it increases without bound with the accuracy of energy 
measurements. We also show in which sense the oscillation formula can be seen as a superposition 
of mass eigenstates with the same energy, in order to clear up the misunderstanding of stationary 
boundary conditions. In addition, we analyze in section |^ the nonstandard oscillation formulas 
derived within field-theoretical models and we disprove them. 

Oscillations of unstable particles are examined in section ^. Nonexponential corrections 
to mixing are computed for the first time in quantum field theory, and they are shown to be 
negligible. We clearly state the correspondence between the field-theoretical treatment and the 
Wigner-Weisskopf and reciprocal basis formalisms. Some more topics are discussed in section ^: 
the normalization of the transition probability, the size of the external wave packets and the in- 
stability of the source. Finally, we examine correlated oscillations (DA<1?NE, KEKB and BABAR 
experiments) in section including an analysis of the possible energy-momentum correlation 
at the source. In the same framework, oscillations of recoil particles are shown to be nonexistent. 

Let us mention four restrictions. First, we shall only consider free propagation in vacuum. 
Interactions with background matter can have a dramatic effect on mixing, as demonstrated 
by the Mikheyev-Smirnov-Wolfenstein (or MSW) effect |5|, |o|, lo^, |l8|]. In practice, the 
MSW oscillations are washed out by averaging over the energy spectrum [227, 204], so that the 
question of intrinsic decoherence discussed in the present report is not very important in that 
case. It is true that a field-theoretical formalism for oscillations with matter interactions should 
be developed to deal with dense astrophysical environments [^34| , |232 , 



I |106| , but this goes 

far beyond the limits of this review. Second, we shall not explain the nonperturbative quantum 
field theory of particle oscillations, as proposed by Blasone, Vitiello and co-workers |^3|, 
This approach is indeed of a completely different nature than the perturbative field-theoretical 
treatments discussed in this report. We shall however comment on their results in section IJ. 
Third, we shall not examine exotic mechanisms, such as oscillations of massless neutrinos due 
to a violation of the Equivalence Principle [p35|| , to resonant spin-flavor precession 
to flavor-changing (or non-universal) interactions [118| or to dissipation |19£, |5^ (comparative 
reviews can be found in [195, |220| , 132[). Finally, we consider neither the gravitational effects 
on oscillations [257, 194] nor the possibility of extra dimensions [109, p9[. 



2 Oscillations in quantum mechanics 

We first discuss the concept of particle mixing as it is a requirement for oscillations. We then 
present the standard derivation of the oscillation formula, stressing the problematic choice of 
basis for unstable particles, as well as the ambiguities arising in the conversion of an oscillation 
in time into an oscillation in space. This section ends with a review of the attempts to solve 
these problems with more sophisticated quantum-mechanical models. 
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2.1 Mixing in quantum mechanics 

The possibility of particle mixing arises from the fundamental principles of quantum theory. The 
measurement of an observable causes the initial state to jump into an eigenstate of the operator 
associated to this observable. Thus the act of measurement determines a basis of physical states, 
made up of all the eigenstates of the corresponding operator. The different bases associated to 
the different operators are linked by linear transformations; a state expanded in a given basis 
can be seen as a mixing of the eigenstates of that basis. While one gets used to the quantum 
superposition principle as long as wave functions correspond to rather abstract properties of a 
system, it is always a bit shocking to apply it to particles, which are associated in our mind 
with classical concepts. The two-slit experiment with particles is the best known example, but 
particle mixing, appearing when an interacting eigenstate is a superposition of different mass 
eigenstates, is hardly less unsettling. 

The textbook case is given by the neutral kaon system: the decay pattern of a into pious 
is well explained by the superposition of two mass eigenstates, a CP-even state decaying quickly 
into two pions, and a CP-odd state decaying more slowly into three pions. Thus a can be 
seen as a mixing of a Ks {K short) and a Kl (K long), having well-defined masses and widths 
(the so-called mass eigenstate basis). Whereas pionic decays do not allow to identify the final 
strangeness {K'^ or K'^), it is possible to do so in the decays of kaons into Treu, with the result 
that one observes a nonnegligible probability to detect a transition from the initial into a 
K^. This is the phenomenon referred to as particle oscillation. The flavor- mixing transition 
probability oscillates indeed with the propagation distance, unless the oscillation is averaged by 
some mechanism, in which case the transition probability simply follows the exponential decay 
law for each mass eigenstate. 

Particle oscillations can be compared to a 2-slit interference experiment, which requires 
firstly that several paths must be available to a particle, and secondly that these paths should 
not be distinguishable by the experimental conditions. In particle oscillations, the different mass 
eigenstates are the analogue of the different paths, since they are eigenstates of the space-time 
evolution operators P and H. Of course, their other quantum numbers should be equal. The 
equivalent of the second requirement is that the masses should not be distinguishable by the 
experimental apparatus. Under these conditions, the probability to observe the propagating 
particle in a given state, other than a mass eigenstate, oscillates in space (and in time) with 
a wavelength inversely proportional to the mass difference, because of the interference between 
the mass eigenstate amplitudes. Thus particle oscillations occur if the initial and final states are 
not mass eigenstates and if it is impossible to ascertain which mass eigenstate has contributed 
as an intermediate state. This last assumption implies that mass differences should be much 
smaller than the typical energies involved in the processes, with the result that the oscillation 
wavelength is macroscopically large and the particle oscillates over macroscopic distances. 

Of course, oscillations are only observable in processes where the propagation distance 
is of the order of, or much larger than, the oscillation wavelength. Otherwise it is perfectly 
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adequate to work with interaction (or flavor) eigenstates, as it has been done until now with 
neutrinos in accelerator experiments. Note that one often talks about particle oscillations even 
if the oscillation term is averaged by some mechanism, such as an integration over the energy 
spectrum, though this situation could be more appropriately called incoherent mixing. 

Let us be more precise about the basis pertaining to the oscillation phenomenon. The 
Hamiltonian should be such that it can be split into a propagation Hamiltonian Hpropag, de- 
scribing the free propagation of the particle, and an interaction Hamiltonian Hint, describing 
the interactions producing the propagating particle. These two subsets of the Hamiltonian can 
be distinguished if there exists a quantum number, called flavor, conserved by Hint, but violated 
by Hpropag- The appropriate basis for the production of the particle is made up of eigenstates 
of both Hint and the flavor operator. It is called the flavor basis. The basis relevant to the 
propagation of the particle is the one in which Hpropag is diagonal and is called the propagation 
basis, or mass basis. This basis does not coincide with the flavor basis, since the flavor operator 
does not commute with Hpropag- 

The simplest example is given by stable neutrinos, in which case the Hamiltonian Hpropag 
includes the mass matrix, whereas the Hamiltonian Hint includes the weak interactions of the 
neutrinos with the bosons Z and W. The flavor is here the leptonic number, which can be 
electronic, muonic or tauic and maybe sterile. This number is conserved by the weak interactions, 
but is violated by a non-degenerate mass matrix which mixes neutrinos of different flavors. 



As regards unstable particles, the Wigner-Weisskopf method |191| , |172| , |209|] is used to build 
a non-Hermitian effective Hamiltonian Hpropag, with the aim of describing both propagation 
and decay. Hint includes the interactions involved in the production of the oscillating particles. 
In the case of neutral K (or B) mesons, Hpropag includes, on the one hand, the degenerate 
effective mass matrix (due to Yukawa, strong and electromagnetic interactions) and, on the 
other hand, the weak interactions causing the decay and the flavor-mixing transitions ^ 
(or B^). Hint is an effective Hamiltonian describing the strong interactions of the mesons; 

it can be built with chiral perturbation theory. The flavor is here the strangeness (or bottomness) 
with values ±1, which is conserved by strong interactions but violated by weak interactions. The 
non-Hermiticity of Hpropag means that the neutral meson system does not evolve in isolation, 
because of its possible decay. 

The relationship between the propagation and flavor basis can be written as 

\Ua{0)>=Y,UajWj{0)>, (1) 
j 

where U is the matrix diagonalizing Hpropag, if this Hamiltonian has been expressed from the 
start in the flavor basis. The states |^'j(0)> are the propagation (or mass) eigenstates, at time 
t = 0, with well-defined masses and widths, belonging to the propagation basis, whereas the 
states |i'a(0)> are the flavor eigenstates. 

The properties of the matrix U depend on the Hamiltonian Hpropag- In the case of a stable 
oscillating particle, the matrix U arises from the diagonalization of the mass matrix and is 
unitary. In the case of an unstable oscillating particle, Hpropag is necessarily non-Hermitian, as 
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explained above, so that the time evolution is nonunitary, i.e. the oscillating particle is allowed 
to decay. The matrix U diagonalizing that Hamiltonian is generally not unitary, although it can 
be unitary if some symmetry makes the Hamiltonian normal]^ |^l| . 

The computation of an amplitude requires the knowledge of the scalar product values of 
the eigenstates. In the flavor basis, the following orthogonality properties hold: 

<M0)\ua{0)>=6^p, (2) 

since Hint is Hermitian and commutes with the flavor operator. In the propagation basis, 
Hpropag is Hermitian only if the oscillating particle is stable. Thus the scalar product of the 
flavor eigenstates must be defined by transforming the flavor scalar products (Q): 

7 

These scalar products become orthogonal if the matrix U is unitary, which is true if the par- 
ticle is stable or, more generally, if Hpropag is normal. In the example of neutral kaons, the 
approximation of CP symmetry makes Hpropag normal and U unitary. Since the violation of this 
symmetry is of the order of 10~^, the right-hand side of Eq. (^) is of the same order in the neutral 
kaon system. This lack of orthogonality implies that the states | i'i{0)> can be normalized in 
different ways. The normalization factors do not matter as long as they do not appear in the 
final result. A normalization problem arises however in the case of unstable particles if we try 
to compute an amplitude involving a mass eigenstate in the initial or final state. For example, 
the amplitude Ktttt \ Htotai \ Kl > depends on the normalization chosen for the state | >■ In 
that case, the normalization ambiguity is of the order of the direct CP violation parameter e', so 



that it cannot be neglected |61|. This problem does not arise in the field-theoretical approach, 
where the kaons are considered as intermediate states, which are not directly observed. 

It must be stressed that the orthogonality properties in the flavor basis are valid only if all 



mass eigenstates are kinematically allowed |63, 22C, M. For example, let us suppose that there 



are four stable neutrinos such that = for i = 1,2,3 and 7714 ^ 1 GeV. Under 1 GeV, the 
flavor scalar products become 

< MO) I l^a(O) >= 5a/3 - Ua4 U*pA ■ (4) 

This difficulty is a first indication of the problems arising in the definition of a flavor basis. Is it 
really possible to interpret a flavor eigenstate as an observable particle? In quantum mechanics, 
each stable particle can be associated with an irreducible representation of the Poincare group, 
with a given mass and a given spin. In that framework, the flavor eigenstates do not correspond 
to particles, unless they are degenerate in mass. The problem is not easier to solve in quantum 
field theory. Whereas creation and annihilation operators of mass eigenstates are well-defined, 
the transformation ([l|) does not yield creation and annihilation operators of flavor states that 



satisfy canonical commutation relations |142]. One could thus wonder whether flavor eigenstates 



is normal if \H,H'']=0. 
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really exist or not. A phenomenological argument in favor of their existence is that, for obvious 
practical reasons, no one has yet devised an experiment involving only one mass eigenstate of a 
system of oscillating particles^. On the contrary, oscillating particles are produced and detected 
in a superposition of mass eigenstates formally equivalent to a flavor eigenstate. In that sense, 
flavor eigenstates are observed, whereas mass eigenstates are not. 

The case of an unstable particle differs from the stable case in two respects. First, the decay 
widths of the propagating eigenstates can vary greatly, allowing the isolation of a nearly pure 
beam of the longest-lived eigenstate (so the phenomenologist is tempted to declare this state 
'observable'). Second, there is no state corresponding to an unstable particle in the Hilbert 
space of physical states. These particles are best described by S-matrix poles in quantum field 
theory (so the theoretician is tempted to declare these unstable states 'not observable'). 



2.2 Oscillations with plane waves 

Many articles and books give derivations of the oscillation formula with plane wave states (see 
for example Refs. \6G, |208|]). In this section and in the next, we examine this method in its fine 



details, emphasizing subtle points and doubtful assumptions. 

Following the rules of quantum mechanics, the transition probability between an initial 
state of flavor a at time t = and position x = and a final state of flavor (3 at time t and 
position X is given by 

^Q,^/3(t,x) =< 1^/3(0) \exp{-iHpropagt + iP ■ x)| Z^a(0)> , 

where Hpropag is the Hamiltonian operator and P is the generator of translations. The states 
I z^o,/3(0)> are taken to be energy- momentum eigenstates, so that a plane wave treatment will 
follow. This amplitude can be computed with the help of the change of basis (|l|): 

Aa-.f3{t, X) = < Up{0) I Uaj e-'^^-'<^ \ V, (0) > 

j 

= <z.;3(0)|^C/„,e-^<^^-^^C/rV,(0)>, (5) 



where the phase is given by 

(t>j = Ejt - Pj • x , (6) 

and the decay term is defined by 

The symbols Ej, pj, rrij and Tj stand for the energy, momentum, mass and width of the mass 

eigenstate 1 1/^(0) >. The factor pj • x is often dropped, either because the mass eigenstate is 

assumed to be in its rest frame (the time t is then the proper time), or because of an assumption of 

*This might have been possible for neutral kaons, had the CP symmetry been exact in that system, making 
the propagation basis coincide with the CP basis. Decay channels respecting CP {K 2it, Stt) would have then 
allowed a nonambiguous determination of the propagating eigenstate. 
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equal momenta for all mass eigenstates (the factor pj • x then does not appear in the probability) . 
These assumptions cannot be justified and neither will be assumed in the following derivation. 

If the particle is unstable, Hp^opag is an effective Hamiltonian, computed with the Wigner- 
Weisskopf method in the rest frame of the decaying particle. Hence, the phase depends on the 
proper time, (pj = rajT, so that a boost is necessary to obtain the expressions @ and (0) valid 
in the laboratory frame. However the concept of a rest frame has no meaning in the case of a 
superposition of different mass eigenstates, which may have different velocities. Thus the choice 
of the boost is ambiguous because it is not clear whether a common boost should be used for the 
different mass eigenstates, or whether each mass eigenstate is boosted by a different velocity. 
This difficulty is another indication of the problems arising in the correspondence between a 
flavor state and a particle. It will be seen in section that the choice of the boost has a big 
impact on the value of the oscillation length. 

With the help of the orthogonality property (^), the amplitude (^) can be rewritten as 

Ao^^^it, x) = ^ U^j e-"^^-^^ . (8) 
j 

In the case of K and B mesons, the oscillation formula (^) is sometimes written in the 
following way [|21C|]: 



Aa^pit, x) = 5^ Uo^j e-'^^-^^ < i^iiO) I Uj (0) > . (9) 

The equivalence of this expression with Eq. can be checked with the scalar products @. 

Another way to obtain the transition amplitude (|8|) for neutral mesons is to define a recipro- 
cal basis |23J , 110, 17, 77, p43 |. This method consists in defining two mass bases: the kets belong 



to the in basis and are the right-eigenvectors of the Hamiltonian, whereas the bras belong to the 
out basis (or reciprocal basis) and are the left-eigenvectors of the Hamiltonian. Their Hermitian 
conjugate vectors are not used to write amplitudes. The Hamiltonian can then be expressed in 
a diagonal form: 

Hpropag — ^ ^ | '^j(O) ^in '^j out 

<i/,(0)|, (10) 

i 

where the \j = mj — iTj/2 are the complex eigenvalues of the Hamiltonian in the rest frame of 
the oscillating particle (in the mass-degenerate limit). The new bases are related to the flavor 
basis by 

\l^a{0)> = Ugj I l^j{0)>in , 

j 

<'^a{0)\ = Y.out<'^jm\Uia ■ (11) 
j 

With these notations, a kind of orthogonality property is restored: 

out < M^) I >in= Sij . 
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The oscillation formula obtained with these new bases is identical to Eq. (^), since it is simply 
another way to decompose Hpropag- 

The physical interpretation of the reciprocal basis is not clear and its use is not really 
necessary, since the same result can be obtained without it. Actually, the reciprocal basis 
method is a technical trick which is best understood from a field-theoretical viewpoint. After 



all, the new bases were first introduced |238] in a quantum field theory approach as left- and 
right-eigenvectors of the propagator representing the oscillating particle. As emphasized by Enz 
and Lewis [|110(| , they should always appear through Eq. (|To|), underlining their unphysical and 
intermediate character. For example, the quantity | <TTir \Htotai \ Kl> |^ is not an observable, 
as mentioned in section p.lj . 

In the special case of a Hermitian Hpropag, U is unitary and 7^ = 0. The amplitude then 
reads 



Aa^^it,^) =Y.Uaj e-*^^*+*P^-^ U*p, 



^3 
3 

This is the expression commonly used for neutrinos, except that the term pj -x is often dropped 
for the reasons mentioned after Eq. (^). 

In the general case of a non-Hermitian Hamiltonian Hpropag, the transition probability is 
given by the squared modulus of the oscillation amplitude @ : 

V^^pit,^) = ^ Uo^^Ur^Ul^ (c/-i)*^e-(*-*^)-^-^^' . (12) 
hi 

The particle oscillation phenomenon is manifest in the oscillatory behavior of the interference 
terms in the transition probability (12), with a phase — (pj depending on the point (t,x) of 



detection. 

Regarding antiparticle oscillations, the CPT theorem applied to the amplitude (P) gives the 
relationship 

A^^pit, x) = Ai3-^a{t, x) , 

SO that the oscillation formula for antiparticles is obtained from Eq. (p^ by the exchange of a 



and p. Note that CP violation arises from terms in Eq. (12) breaking the a<-^/3 symmetry, since 
it appears through a difference between |^^^^(t,x)p and |^Q^^(t,x)p. 

2.3 Time to space conversion of the oscillating phase 
2.3.1 Standard oscillation phase 

There has been some controversy about the conversion of the oscillation phase (pi — (pj , appearing 
in Eq. (|l2|) , into a measurable quantity. Since the propagation time is not measured in oscillation 
experiment^, a prescription is needed to get rid of the time dependence in the phase difference. 

The numerous prescriptions proposed in the literature are somewhat confusing. They can be 
classified by expanding the phase around an average energy or momentum. Since the oscillating 



^However time measurements are important according to Okun |217| in an experiment performed at IHEP 
(Serpukhov) ||. 
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particle is on-shell, the energy Ej can be expressed in function of the momentum pj. The phase 
<j)j can then be expanded around an average momentum p, not very different from pj or pk, 
and an average mass m. Although the mass m can be very different from nij or in the 
ultra-relativistic case, the expansion is correct as long as the mass difference (5m| = — is 
small with respect to the energy. The momentum difference 5pj = Pj — p is expected to be of 
the same order than 5m^. In one spatial dimension, the expansion of the phase reads, to first 
order in drrij and 6pj, 



pj + m| tj - Pj X 

Et-px + — rft + {vt-x)5pj + E5tj, (13) 
IE 



with the average energy and velocity defined hy E = \/ p^ + m^ and v = p/ E, respectively. A 
different time tj has been allowed for each mass eigenstate, and expanded around an average 
time t, with 5tj = tj — t. Of course, the following arguments will only be correct to first order in 
5mj/2E, but it is useless to argue about further orders in a fiawed approach such as the plane 
wave treatment. 

All prescriptions leading to the standard oscillation formula set 8tj = 0, i.e. they impose 
that interference only takes place for equal propagation times (and lengths) for the different 



mass eigenstates. This equal time prescription has been explicitly stated |181|, |177|, 178, 201] 



in reaction against articles proposing different detection times [244, p45||. It has also been 



legitimated by an equal velocity prescription |l245| , |101[| , which is seen to be equivalent to the 



previous prescription with the help of the classical relation tj = x/ Vj . However the equal velocity 
condition leads to ^ = which is very unlikely for neutrinos [ |218|] . 

Note that imposing equal times tj = tj, in the laboratory frame, also means imposing equal 



proper times Tj = Tj if the classical relation Tj = yt^ — is used. The last relation implies 
that a boost of velocity Vj = x/tj is used to go from the rest frame of the mass eigenstate mj 
to the laboratory frame. Thus the question of the choice of the correct boost boils down to the 
question of the equality of propagation times in the laboratory frame. 

A second prescription is needed in order to obtain the standard oscillation formula. It could 
be called the classical propagation condition: vt — x = 0. It imposes that the term (vt — x)6pj, 
appearing in Eq. (13), is negligible in comparison with dmjt/lE. This condition can be weakened 



to |ft — a;| <C t. Since plane waves are delocalized in space-time, this condition cannot be justified 
without a more sophisticated treatment, for example with wave packets. 
These two prescriptions are sufficient to derive the following formula: 

(5m?- Ixl |x| , , 

2|p| L^- 

obtained from a three-dimensional generalization of Eq. (^). The oscillation length L"^^ is 
defined by 

47rlDl 
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where Sm'^j = mf — mj is assumed to be positive. The classical propagation condition also 
justifies the substitution t — > -Elxj/lpl in 7,^ (see Eq. (0)). 

The transition probability (^) together with the phase (14) form the standard oscillation 



formula used to fit the experimental data (see [ll54[| for an application to two-flavor neutrino 
oscillations). Of course, the probability ( |l2|) should first be averaged over the energy spectrum 



and over the region of production and detection (see section 9.2.1). If the oscillation is completely 



washed out by these averaging procedures, the transition probability ( |12| ) can be simplified by 
the substitution e~*('^'~'^j) 6ij. 

Although the equal time prescription and the classical propagation condition are sufficient 



to obtain the standard oscillation phase (14), additional prescriptions leading to the same result 
are commonly found in the literature: 

1. The equal momentum prescription is the most common: 5pi = 5pj = 0. As seen above, this 
assumption is not necessary. Moreover, it is impossible to impose experimental conditions 
such that the momentum uncertainty is zero, since the oscillations are destroyed by a 
momentum measurement more accurate than the mass difference. Thus, this prescription 
is groundless. 

2. The equal energy prescription has been recently advocated by Lipkin: 5Ei = SEj = 
|197]] . It has the advantage of avoiding the classical propagation condition since it leads to 



Spj = —5m'j/2p, so that the time dependence completely drops from the phase difference 
(pj ~ 4'k (at least if the equal time prescription is assumed), and the standard oscillation 



phase (14) is directly obtained. 

In principle, oscillation experiments are feasible with a zero energy uncertainty, since 
spatial oscillations are not expected to vanish in that case (note that a quantum field 
treatment is necessary to prove it). In practice, the energy uncertainty is far from being 
negligible and is often of the same order of magnitude as the momentum uncertainty. 
Thus the equal energy prescription is only justified if an extremely small uncertainty on 
the energy is imposed by the physical properties of the process itself. Do we have any 
theoretical reason to expect that this uncertainty is smaller than the mass difference? 
No convincing arguments supporting that assumption have been given until now. For 
example^, Lipkin computes the energy- momentum uncertainties SpK, of the kaon in 
the process vrp" -^AK^ [195]. If the proton is at rest in a lattice, its momentum uncertainty, 



due to the Debye temperature of the crystal, can be estimated at 5pp ~ 10 eV. Lipkin 
estimates 5pK to be of the same order of magnitude, whereas he neglects 5Ek ^ {Sppf' / nip. 
However 5Ek is still much larger than the mass difference Smx ~ 10~^ eV, which is the 
most important mass scale in the experiment. Moreover the pion momentum uncertainty 



^Lipkin gives at least three reasons for the equahty of energies: the strict energy conservation discussed here 



[L97, 198 1, the flavor-energy factorization [155, 198 1 and the stationarity resulting from a time average [197, 198 j 



The last two arguments are discussed in section 
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gives a first order contribution to dEx and should not be neglected. Even if one has shown 
that SE <^ 6p, there is a long way to go to show that 6E <C 5m. 

Another argument against the equality of energies is that it holds only in one particular 
frame. For example, if the energies of the different mass eigenstates are equal in the 
decay of a pion at rest (tt— the energy difference becomes approximately equal to the 
momentum difference if the pion is relativistic [|136[| . Thus the equal energy prescription 
should be shown to be true in the laboratory frame for any experimental conditions. 

However there is a much more reasonable way of looking at the equal energy prescription. 
It consists in seeing it as the result of a time average washing out the interference between 
wave packet components having different energies. There is then no need to prove that 



the energy uncertainty is zero. This argument will be considered in section 2.5 



3. Energy-momentum conservation at the production: first proposed by Winter | 255 |, this 
recurring prescription [75, 146, 201, 244, 245, 14C, 136, 252| allows to compute explicitly 
5pj. It has often been used to show that neither the momenta nor the energies of the 
different mass eigenstates are equal. In the example of the pion decay at rest (vr— >^z^), the 
energies and momenta of the muon and neutrino can be computed exactly if the energy- 
momentum of the pion is perfectly known. However this knowledge is usually not available: 
when the energy-momentum spread of the source is much larger than the mass difference 
Smfj, it is meaningless to compute the exact values of the energies and momenta to order 
Smfj. A more detailed examination of this question requires wave packets instead of plane 
waves, or even better, quantum field theory. 

2.3.2 Non-standard oscillation phase 

Controversial prescriptions leading to nonstandard oscillations formulas involve different propa- 
gation times 5ti / 5tj ^ or, equivalently, different proper times Tj ^ tj. Let us parametrize pj 

by 

pj=p+{p-l)^, (16) 

where /? is a dimensionless number of order unity. The corresponding energy and velocity can 
be written as 

5m? 

^ ^ 2E 



Ei = E + p — ^ 
^ 2E 



The momenta are equal if /j = 1 whereas the energies are equal if /) = 0. As explained above, it 
has also been proposed to determine p through energy-momentum conservation at production, 
leading to a value of p of order unity. In any case, the time difference is computed with the help 
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of the classical relation t 



J — x/ Vj and reads 
6t, = {E^ 



pm 



SrrijX 



Inserting this value in Eq. (13) and using t = x/v, we obtain 

1-p 



Et-px + [1 + p + 



6mjx 
2p 



(18) 



With the equal momentum prescription (p = 1), the corresponding oscillation length will be 
smaller by a factor 2 than the standard value given in Eq. dlD fig^. With the equal energy 
prescription (p = 0), the oscillation length will be smaller by a factor 1 + v^'^ than in Eq. (15). 
Thus the equal energy prescription may also lead to a nonstandard oscillation length, contrary to 
what was claimed in [197]. With the energy-momentum conservation prescription, the oscillation 
length will be smaller than in Eq. ( [l5| ) by a factor depending on the value of p, which depends 
on the energy [p44| . These formulas can also be applied to the case of correlated oscillations, 
such as 0(1020) -^K^K^ or T(4s) -^B^B^. In the center-of-mass frame of the resonance, the 
different time prescription leads to an oscillation length smaller by a factor 2 than the standard 
value, since the equality \pi\ = \pj\ valid in that frame leads to p = 1 [245]. 

Another disturbing consequence of the different time prescription is the oscillation of recoil 
particles, for example A in ■K~p^KK^ 1244], or the muon in vr— ]p 



247]. This is easily 



seen by applying Eq. (|T^) to the recoil particle. Although = (the recoil particle has only 
one mass eigenstate), (pi ^ because 5ti ^ 5tj. The oscillation of particles having only one 
mass eigenstate is unacceptable since it leads to non-conservation of the detection probability 
of this particle. This is not the case when there are several mass eigenstates, as the sum of 
the detection probabilities of the different mass eigenstates is always equal to 1 for a given 
propagation distance. 

The treatment of neutrino oscillations m ir ^ pv proposed by Field 112 ] also resorts to 
the different time prescription. This author claims that the different neutrino mass eigenstates 
are detected at the same space-time point but are produced at different space-time points. 
Following Field, the oscillation phase would not only receive a contribution from the neutrino 
path difference but also from the path difference of the source which decays at different times. 
Field computes the first contribution with the energy-momentum conservation prescription, 
which gives for a source at rest an oscillation phase larger than the standard result by a factor 
2 (see Eq. (18) with v = 1). The second contribution, which only appears in Field's article, is 
computed in the same way. The oscillation phase obtained by Field therefore differs from the 
standard result (|l^). Field's method also leads to the prediction of muon oscillations because 
of the use of the different time prescription^. 

Another example of the different time prescription can be found in Malyshkin and Kulsrud's 
analysis ]204] of the time variations of the solar neutrino flux, which leads to a result different 



^Field's criticism of wave packets models will be examined in section 2.5 
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from Parke's formula |224]. As we do not consider here oscillations in matter, this new oscillation 
formula will not be discussed further, although our comments on the different time prescription 
also apply to that case. 

The different time prescription has been strongly criticized by several authors within the 



plane wave framework [201, 181, 177, 178, 197, 198]. Their argumentjj, which can be stated 



as 'interference only occurs between states taken at the same space-time point', does not hold 



when examined in a wave packet or field-theoretical model. Kiers, Nussinov and Weiss [182, 182] 
have indeed shown that the coherent character of the detection process allows wave functions 
at different space-time points to interfere. In particular, a long coherent measurement in time 
may be used to revive oscillations, even after the mass eigenstate wave packets have completely 
separated spatially (see section |2.5| ) . The question of the correct time prescription is thus subtler 
than it seems at first sight. The crux of the matter is to take into account the production and 
detection processes. In this way, each increment in the phase associated to the propagating 
particle, due to a slightly different production (respectively detection) point, is cancelled by a 
decrement in the phase of the wave packet of the source (respectively detector) ]140| . For exam- 
ple, this cancellation can be implemented in the intermediate wave packet model (see section |2^ ) 
by computing the total amplitude as an overlap of the propagating wave packet with the source 



and detector wave packets [ |139| , 137, |138| ]. Note however that this mechanism is much more 
natural in a quantum field model such as the external wave packet model (see section |^ , where 
the amplitude is integrated over all possible microscopic production and detection points, with 
the result that the phase depends only on average (i.e. macroscopic) production and detection 
points. 

In conclusion, neither the equal time prescription nor the different time prescription can be 
justified in the plane wave approach, although the choice of the prescription has an important 
effect on value of the oscillation length. Moreover other prescriptions, such as the classical 
propagation condition and the equal energy prescription, cannot be understood within the plane 
waves formalism. A wave packet or quantum field treatment is thus inescapable. Let us also insist 
on the dubiousness of the arguments using energy-momentum conservation. In most cases, they 
are invalidated by the energy-momentum spectrum of the source: different energy-momentum 
components of the source can contribute to different mass eigenstates. 

2.4 Problems with the plane wave treatment 



In section 2.3, it has been shown that the plane wave treatment of particle oscillations cannot deal 
in a satisfactory way with the time dependence of the oscillating phase. Besides, this approach 
implies a perfectly well-known energy-momentum and an infinite uncertainty on the space-time 
localization of the oscillating particle. Oscillations are destroyed under these assumptions |176[] . 
On the one hand, the perfect knowledge of the energy-momentum allows to determine which 
mass eigenstate propagates. On the other hand, the spatial delocalization makes impossible the 



more convincing (though indirect) argument has been given by Giunti |137|: the use of the different time 
prescription in the double-slit experiment leads to a wrong interference pattern. 
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measurement of the oscillation length. A correct oscillation formula should include observability 
conditions in such a way that the oscillation term vanishes if either 



1. the energy-momentum uncertainty is smaller than the mass difference between the inter- 
fering mass eigenstates, or 

2. the oscillation length is of the same order, or smaller, than the uncertainty on the position 
of the source or of the detection point of the oscillating particle. 

Another kind of problem is not specific to the plane wave treatment, but affects all ap- 
proaches where the oscillating particle is considered to be directly observable. On the one hand, 
flavor eigenstates are ill-defined for stable particles, because we do not know how to define cre- 
ator and annihilation operators of fiavor states satisfying canonical (anti) commutation relations 



1 142]. On the other hand, mass eigenstates are ill-defined for unstable particles, since they are 
in general not orthogonal (except if the Hamiltonian is normal, see section |2.lD . The solution to 
these problems simply consists in considering the oscillating particle as an intermediate state. 
Actually this stand reflects well the experimental situation, where the oscillating particle is not 
directly observed. One rather detects the particles in interaction with the oscillating state, both 
at the source and at the detector. The flavor transition probability should thus be computed 



with observable particles as initial and flnal states |187| , |142| , |145| , |237|1 . 

Finally, we should mention two other problems regarding unstable particles. These difficul- 
ties arise because of the nonrelativistic Wigner-Weisskopf method used to compute the effective 
Hamiltonian. First, interference between different mass eigenstates is forbidden in nonrelativistic 



quantum mechanics (Bargmann superselection rule [48, 175, |133| ]). The argument is the follow- 
ing. The invariance of the Schrodinger equation under Galilean transformations determines the 
transformation law of a quantum state: it is multiplied by a phase factor depending on the mass 
and space-time position of the state. Thus different mass eigenstates transform differently, so 
that the relative phase in a superposition of such eigenstates is not conserved under Galilean 
transformations. Therefore a coherent superposition of different mass eigenstates is forbidden. 
Second, unstable particles cannot be consistently described in nonrelativistic quantum mechan- 
ics for the same reason at the origin of the Bargmann superselection rule: transitions between 



different mass eigenstates are forbidden |17S]. Unstable states cannot be considered as asymp- 



totic states. Thus they do not appear in the Hilbert space of physical states and must be treated 



in quantum field theory where they appear as complex poles of the full propagator [253]. 



2.5 Intermediate wave packets and other improvements 

Some of the problems of the plane wave treatment are solved by the intermediate wave packet 
model, in which a wave packet is associated with each propagating mass eigenstate. Note that 
this model is usually discussed with respect to neutrinos, i.e. in the relativistic limit. Nussinov 
was the first to put forward the existence of wave packets as the cause of a coherence length 
beyond which oscillations vanish [|216([ . Oscillations with wave packets were then studied in 
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Figure 1: Decoherence due to a wave packet width Ox larger than the oscillation length. 




Figure 2: Decoherence due to the non-overlapping of the different mass eigenstates wave packets. 



detail by Kayser |17(:]. The oscillation formula was later explicitly computed with Gaussian 
wave packets by Giunti, Kim and Lee | 141 , 139 , 137 1. 

The intermediate wave packet model shows that oscillations vanish if Ox ^ ^°if^ ^ 
the uncertainty over the position is larger than the oscillation length (Fig. 1). For minimal 
uncertainty wave packets, this condition can be rewritten as < 5mf^/2p^ i.e. oscillations 
are forbidden if the momentum spread of the wave packets is smaller than the mass differ- 
ence between the interfering eigenstates. Thus oscillations are destroyed by energy-momentum 
measurements aiming to determine which mass eigenstate propagates. 

The other main result of the model is the existence of a coherence length beyond which 
oscillations vanish [216, 17€, 141 1. Its usual explanation is that wave packets associated to 
different mass eigenstates have different group velocities. Hence, wave packets progressively 
separate, and interference disappears when they do not overlap anymore (Fig. 2). We shall see 
in section |7.1.3| that dispersion is also at the origin of the coherence length. 

In the intermediate wave packet model, the classical propagation condition, \vt — x\ <C i, 
is automatically implemented by the space-time localization of the wave packet. However the 
question of the equal or different time prescription is not solved, since the times of production 
and detection are not specified by the formalism. As regards the question of the equality of 
energies or momenta, it is in principle no longer relevant, since the energy-momentum is spread 
out. In fact, the oscillation length depends only on the zeroth order values (in a (5m?j -expansion) 
of the average energies and momenta, because of the cancellation between the time and space 
parts of the phase (see Eqs. (^) and (|T6| ) with pj equal to the average momentum). The time 
dependence of the probability is dealt with by averaging over time [|141| . It is thus not necessary 
to worry about the equality of energies or momenta. 

Nevertheless there is still much argument about the equality of the energy spectra, or the 
equality of the average energies, momenta or velocities. For example, Grossman and Lipkin 
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|155| ] imposed a same energy spectrum for all mass eigenstates, so that the probability of finding 
a particle with the wrong flavor vanishes for all times at the position of the somxe. It can be 
objected that this boundary condition is automatically satisfied without imposing an identical 
energy spectrum, since the probability to detect the oscillating particle at the source becomes 
negligible once the wave packet has left it. Equality of average velocities has been proposed 



by Takeuchi, Tazaki, Tsai and Yamazaki pSO] , 251], but is very unlikely for neutrinos, as noted 
by Okun and Tsukerman [218|, because the ratio of the average energies Ei/Ej should then be 
equal to the mass ratio mi/rrij. 

Actually there is another way of looking at this question. Starting with arbitrary wave 
packets, it is simple to check that interference occurs only between wave packet components with 
the same energy, at least if the oscillation process is strictly stationary. This line of thought has 



been advocated by Sudarsky et al. [^49(| (within a field-theoretical model) , by Lipkin 1 197 , 195 ] 
(with wave packets) and by Stodolsky [ |248| ] (within a density matrix formalism). Let us explain 
it with the intermediate wave packet model. In one dimension, the wave function corresponding 
to the mass eigenstate rrii can be written as 

^i{x,t)= [ d^^^(^)e-*^(*-*o)+*(^-^•o^ 



where ki = y E'^ — m? and {to, xq) are the coordinates of the source. The fiavor-mixing transition 
probability is a linear superposition of terms ipi{x, t)ip*j{x, t). If the time of emission is unknown, 
the probability is averaged over to, so that it becomes a superposition of terms like 

j dtoMx,t)ij*ix,t) = 27ry d£;^i(^)A*(S)e*(^'-^^)(^-^'«) . (19) 

This equation makes clear that interference occurs only between wave packet components with 
the same energy. The time-averaged transition probability in the wave packet model is equivalent 
to an incoherent sum over energy eigenstates ('incoherent' means a sum in the probability, not 
in the amplitude). Eq. (|l^) can thus be seen as a sum over plane waves of energy E and 
momentum ki, weighed by the factor Ai(E)A*{E). The equivalence between a wave packet and 
a plane wave decomposition of a stationary beam was already observed by Dicke and Wittke 
1 104] in connection with electron-interference experiments, and by Comsa in connection with 



neutron- interferometry. The same issue was recently raised by Kiers, Nussinov and Weiss | 182 ] 



and by Stodolsky | 248 |, this time in connection with neutrinos. 

Of course the time integral in Eq. (|l9|) only yields a delta function if the time interval 
is infinite, i.e. if the process is strictly stationary. In practice, there is always some available 
information on the time of emission. If the time integral is limited to the interval AT, the delta 
function becomes a narrow peak of width l/AT, which can be neglected as long as it is small 
in comparison with the mass difference Smfj/2E. This condition can be written as AT 3> T°j^'^, 
where Tff^ is the oscillation time. 

A few more results have been obtained with wave packets. Kiers, Nussinov and Weiss have 
used the equivalence between wave packets and plane waves, explained above, to show that it 
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is always possible to increase the coherence length by more accurate energy measurements at 



the detector [1821. After Giunti and Kim 1 139 1 showed how to take into account the detector's 



momentum uncertainty in the intermediate wave packet model, Giunti |137] did the same with 
the detector's energy uncertainty, so as to reproduce Kiers, Nussinov and Weiss' prediction. 
Wave packets have also been applied by Nauenberg to study correlated oscillations of neutrinos 
or mesons pl2(1 , and to the propagation of neutrinos in matter |143| , p25|| . Finally, it is worth 
mentioning that the wave packet approach is explained in great detail in Kim and Pevsner's 



book [|185[| , is discussed by Fukugita and Yanagida |12g|] and has been reviewed by Zrajfek p5Sf| . 

The intermediate wave packet model has been recently discussed by De Leo, Nishi and 
Rotelli ]102| ], who recover the standard oscillation formula under two conditions. The first con- 
dition, which they call 'minimal slippage', is equivalent to the existence of the coherence length 
discussed above, i.e. mass eigenstate wave packets should overlap at detection otherwise there is 
decoherence. The second condition, which they call 'non instantaneous creation' or 'pure flavor 
creation event-wise', means that interference occurs between wave packet components corre- 
sponding to different initial wave packet points, so that the standard oscillation phase (14) is 
obtained. This mechanism, which was already explained by Giunti and Kim [|140| (see end of sec- 
tion |23]2), is automatically included in quantum field treatments where the oscillating particle is 
coupled to a source and detector (see section ^). In the same article, De Leo, Nishi and Rotelli ar- 
gue that the oscillation formula might depend on the wave packet shape, for example if the wave 
packet is a superposition of two well-separated Gaussian peaks, each having its own phase fac- 
tor. In that case, there would be a succession of coherence/decoherence/coherence/decoherence 
regions as the mass eigenstates wave packets shift one with respect to the other. Moreover, the 
constant phase difference between the two Gaussian peaks crops up in the oscillation phase in 
the second coherence regime. While surprising phenomena indeed arise for such special wave 
packet shapes, it seems difficult to set up an experiment in which a specific wave packet shape is 
maintained for all the particles within the produced flux. In practice we should average over all 
wave packet shapes, with the result that only the typical width of the wave packet will have an 
effect on the final oscillation formula. All other effects due to the different shapes of the wave 
packets cancel each other out. Besides, the analysis of oscillations in a faraway region (such as 
the second coherence regime mentioned above) should take into account the fact that dispersion 
wipes out oscillations at large distances (see section 7.1.3 ). 

The representation of the oscillating particle by a minimal uncertainty wave packet has been 
criticized by Field [112|, who claims that that oscillating particles follow classical space-time 
trajectories and that the oscillation phase is due to the propagation time difference between 
different mass eigenstates (see section |2.3D . One of his main arguments is that the neutrino 
momentum spread in vr ^ /uz/ is of the order of decay width F^ of the muon, so that neutrino 
oscillations should be completed suppressed by decoherence in wave packet models. Moreover, 
the Heisenberg relation is violated if one admits that the neutrino spatial uncertainty is of the 
order of the cm. These arguments do not hold since the neutrino momentum spread refiects the 
momentum uncertainties of the pion source and of the emitted muon, which are much larger 
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than r^. Finally, Field criticizes models where the neutrino wave packet is coupled to a wave 
packet source, by claiming that the initial quantum state is not the same for all amplitudes, so 
that interference is forbidden. This criticism is not justified, since it makes no sense to cut the 
source wave packet into disjointed parts which are then interpreted as different initial states. 

Although they solve some of the problems of the plane wave approach, intermediate wave 
packets are not satisfying for many reasons: 

1. The question of the existence of a Fock space for the weak eigenstates remains. 

2. Oscillating particles are not, and cannot be, directly observed. It would be more convincing 
to write a transition probability between observable particles, involved in the production 
and detection of the oscillating particle. 

3. The wave packet shapes, sizes and normalizations are arbitrary. In particular, the ampli- 
tudes of production and detection are not taken into account. This is not a problem in 
the relativistic limit or for nearly degenerate masses, in which cases these amplitudes can 
be factorized from the sum over the mass eigenstates. However this factorization is not 
possible for nonrelativistic particles (unless they are nearly degenerate in mass). 

4. The coherence length depends on the difference between the group velocities of the wave 
packets, the exact values of which are unknown in the nonrelativistic case0. There is 
no reason to believe that the group velocities could be determined, to order 6mfj/2p, by 
energy-momentum conservation at the source, contrary to what is claimed in |141| . 

5. The coherence length increases with the precision of the energy-momentum measurements 
at the source and at the detector. Whereas the momentum precision is limited by the 
condition ap > 6mf-/2p, there is, a priori, no such limitation on the energy accuracy and. 



therefore, no bound on the coherence length. Although it has been claimed |13£] that such 
a bound exists, the intermediate wave packet model cannot settle that question, because 
the energy uncertainty has to be inserted by hand into the model. 

6. Experiments measure a particle flux, not a particle density. The transition probability 
must thus be converted to a flux density, involving inverse velocity factors. These enhance 



nonrelativistic mass eigenstates and strongly skew the probability 141 |. This effect is due 
to the neglect of the amplitudes of production and detection and to an analysis restricted 
to one dimension. 

7. Although it has not been controversial in the wave packet model, why should the times 
appearing in the different wave packets be identical? A justification of this equal time 
prescription is necessary to rule out the nonstandard oscillation formulas discussed in 



section 2.3 



''Eq. ( |l7| ) can be interpreted as giving the group velocity of a mass eigenstate wave packet. In the relativistic 
limit {v 1), the p-dependence drops from this equation, so that the velocity can be determined up to order 
5m1j/2p without needing to know p. 
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8. It is not possible to analyze in the wave packet model the corrections to the oscillation 
probability due to the possible instability of the propagating particle. 



9. The influence of the source instability on the observability of the oscillations cannot be 
studied in that framework. 

A quantum-mechanical model of neutrino oscillations solving the first three problems has 
been put forward by Rich |237]. This author considers the oscillating particle as an intermediate 



state and computes the probability transition between initial and final observable states with 
second-order time-dependent perturbation theory of standard quantum mechanics. Since this 
model takes into account the production and detection processes, it has the other advantage of 
not requiring the equal time prescription. The spirit of Rich's model is similar to the stationary 
boundary condition models resorting to field-theoretical methods (see section |3.1j) , but it has 
the disadvantage of being nonrelativistic. 

Ancochea, Bramon, Muhoz-Tapia and Nowakowski have tried to solve the difficulty of 
converting the probability density for neutral kaons into a flux [^] . In their nonrelativistic wave 
packet model, they construct flavor probability currents associated to a Schrodinger equation. 



A problem arising in this approach is that the flavor currents are not conserved |255]. Moreover, 
no one knows how to extend this method to the relativistic case. 

Sassaroli has proposed an hybrid model for neutrino oscillations, going half-way to quantum 



field theory, in which a coupled system of two Dirac equations is quantized [|239f| . However 
boundary conditions cannot be applied consistently, unless lepton flavor wave functions are 
considered as observable and the relativistic limit is taken. This difficulty was already noted by 



Giunti, Kim and Lee [142 



The review of intermediate wave packet models would not be complete without mentioning 
its latest and most sophisticated version, as proposed by Giunti |138|] . Instead of representing 
the oscillating particle by a superposition of arbitrary wave packets, Giunti computes the form 
of the intermediate wave packet from basic principles. More precisely, it consists in creating the 
intermediate wave packet from quantum field interactions between the external wave packets 
involved in the production process of the oscillating particle. The wave packet then evolves in 
space-time, before interacting with quantum field interactions at the detector. This interacting 
wave packet model is very close in method, spirit and results to the external wave packet model 



at the core of the present report. It will be further discussed in section 5.2. 



3 Oscillations in quantum field theory 
3.1 Review of the literature 

Few authors deny that the most rigorous treatment of oscillations is done in the quantum field 
theory framework. However, although the quantum field computations in the literature all 
reproduce in some limit the naive quantum- mechanical formula given by Eqs. (|l2|) and (|T^), 
there is not yet an agreement in which respect they differ from the naive formula. There are 
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two reasons for the lack of agreement between the existing quantum field derivations of the 
oscillation formula: first, the different authors use different physical assumptions, and second, 
they use different approximation schemes to compute the transition probability. 

The field-theoretical approach to particle oscillations is quite old. Already in 1963, Sachs 
1 235] applied S-matrix methods to neutral kaon interference^^, using a model developed earlier 
with Jacob |16(:] for unstable particles. Let us call this model, as applied to non-oscillating 
particles, the Jacob-Sachs model. Its application to systems of mixed particles will be called 
the external wave packet model. In this model, the particle to be studied is represented by its 
propagator; it propagates between a source and a detector, where wave packets representing the 
external particles are in interaction. Much later, Sudarsky, Fischbach, Talmadge, Aronson and 



Cheng 1 249 1 studied the influence of a spatially varying potential on the neutral kaon system. 
They resort to a one-dimensional model similar to Sachs' but do not specify the contour of 
integration, so that the finiteness of their final expressions is not guaranteed. Neither Sachs nor 
Sudarsky et al. studied the observability conditions of oscillations. Another simplified model was 
proposed by Beuthe, Lopez Castro and Pestieau |61], with the aim of modelizing experiments 
at CPLEAR and DA<^NE. Their model is not satisfying since they use external wave packets 
localized in time but not in space, which does not correspond to actual experiments. 

As regards neutrinos, Kobzarev, Martemyanov, Okun and Shchepkin analyzed neutrino 
oscillations with a bare-bones quantum field model [187]: the source and detector are infinitely 
heavy nuclei, so that the propagation distance is perfectly known, whereas the propagation time 
is left undetermined. No constraints on oscillations are discussed in that article. Next came an 



important article by Giunti, Kim and Lee [142], showing that it is impossible to build a Fock 
space for flavor states, because the mixing of the ladder operators for mass eigenstates does not 
yield flavor ladder operators satisfying canonical (anti) commutation relations. This observation 
strikes a blow to the quantum-mechanical wave packet approach, which should not be mistaken 
for the quantum field model with external wave packets, as it is sometimes the case in the 
literature. Nonrelativistic corrections to the neutrino propagation are explicitly computed in 



Ref. [142] in a few examples. 

In a pioneering article, Giunti, Kim, Lee and Lee [145[ studied neutrino oscillations within 

a Gaussian external wave packet model. They derive a localization condition (no oscillations 

if L°^^ < ax) and a coherence length beyond which oscillations vanish. These conditions agree 

with those obtained in the quantum-mechanical picture |216, 176, 141]. A later paper by Giunti, 

Kim and Lee ]144] contains essentially identical results as in Ref. ]14£], but with more generality. 

Cardall used the same model [^5|, paying greater attention to the normalization of the event 

rate and to the spin structure. A model close in spirit to the external wave packet model was 

proposed by Kiers and Weiss ] 183 ]. These authors couple the oscillating neutrino with localized 

^^Quoting Sachs, from ]238t : 'The question of whether it is the momenta or the energies of these particles that 
are to be taken equal in the treatment of interference phenomena has often been raised. [Our] method provides a 
clear answer to this question since it is based on an analysis of the phenomena in terms of wave packets'. Nearly 
40 years after, new articles continue to appear, wondering about the equality of energies or momenta. 
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source and detector, which are ideahzed by oscillators. They show how the coherence length 
increases with the energy precision at the detector. Note that they had already predicted this 
phenomenon using elementary quantum mechanics |182]. 

Another kind of model arises from the use of stationary boundary conditions, leading to a 
unique value for the energy of the oscillating particle. The simplest example is the Kobzarev 
et al. model mentioned above [187]. Going a bit further, Grimus and Stockinger proposed a 
model with external particles represented either by bound states or by plane waves. They obtain 
a localization condition but no coherence length. In two other papers, Grimus, Mohanty and 



Stockinger |151, 152] studied the influence of an unstable source on the oscillations, with the 
model of the previous paper modified by a Wigner-Weisskopf approximation for the unstable 



source. loannisian and Pilaftsis analyzed neutrino oscillations ]163] within a scalar version of 
the Grimus-Stockinger model. They claim to have found a novel form of neutrino oscillations at 
short distance, which they call 'plane wave oscillations'. Stationary boundary conditions were 
also used by Cardall and Chung j^] to study the MSW effect in quantum field theory. Note 
that a quantum field derivation of the MSW effect has been proposed by Mannheim |205[] , who 
represents the interaction with matter with an effective potential, whereas this potential has 
been derived with finite temperature field theory by Notzold and Raffelt ]214] and by Pal and 
Pham [ p23[ . 

Campagne ]|8^ studied a neutrino source decaying in flight, such as a relativistic pion, 
with field-theoretical methods. He sidesteps external wave packets by limiting arbitrarily the 
interactions regions with the help of step functions, but it is only a trick to replace complicated 
interactions at the source and detector by simple stationary boundary conditions. Decay in 
flight has also been studied by Dolgov ]|105(] , within a simplified external wave packet model 
modified by a Wigner-Weisskopf approximation. 

Another approach [source-propagator models) consists in using the propagator in configu- 
ration space coupled to a source but not to a detector. Srivastava, Widom and Sassaroli chose 
this method to modelize correlated oscillations of two kaons [p45|] and neutrino oscillations ]247]. 
Since external wave packets are absent, time-space conversion problems cannot be avoided and 
lead to nonstandard oscillation lengths or recoil oscillations (see section |2.3[ ). Shtanov also used 
a source-propagator model, and claims to have found a strong dependence of the oscillation 
formula on the neutrino masses if the source and detector are very well localized in space-time 



1 24^ ]. 

In a completely different line of thought, Blasone, Vitiello and other researchers have at- 
tempted to define a Fock space of weak eigenstates and to derive a nonperturbative oscillation 
formula. The main results of these studies are summarized in J 72]. Note that a previous 
formalism developed in |73, 14, 67, 158] lead to an oscillation formula not invariant under 
reparametrization ]12g|, 74, 130, 69] and was replaced by a more satisfying theory in ]71, |3], 



which was further developed in ]72, 167, |16^ , [70t] . This new theory aims at defining flavor cre- 
ation and annihilation operators, satisfying canonical (anti)commutation relations, by means of 
Bogoliubov transformations. As a result, new oscillation formulas are obtained for fermions and 
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bosons, with the oscination frequency depending surprisingly not only on the difference but also 
on the sum of the energies of the different mass eigenstates. Apart from the speculative nature 
of the enterprise, the drawbacks of the approach are the dependence on time, not on space, 
of the oscillation formula (Lorentz covariance is broken), as well as the lack of observability 
conditions. Although these studies are very interesting from a fundamental point of view, it is 
not obvious whether the new features of the Blasone-Vitiello oscillation formulas are observable 



in practice. Since these new oscillation formulas tends to the standard oscillation formula (12) 
in the relativistic limit or if the mass eigenstates are nearly degenerate, we can focus on the 
case of a nonrelativistic oscillating particle having very distinct mass eigenstates. In that case, 
p ~ 5m? /2E, so that either ap < 5m? /2E or p < Up. Under these conditions, the quantum 
theory of measurement says that interference between the different mass eigenstates vanishes. 
Once the oscillation terms have been averaged to zero, the Blasone-Vitiello formulas do not differ 
anymore from the standard oscillation formula (0). Therefore, the Blasone-Vitiello formalism 
does not seem to be relevant to the phenomenology of oscillations on macroscopic distances. 
This observation does not detract from the theoretical worth of that approach. 

All the above models, whether for mesons or for neutrinos, can be grouped in four categories: 
external wave packet models, stationary boundary conditions models, source-propagator models 
and Blasone-Vitiello models. In the following sections, the connection between external wave 
packet models and stationary boundary condition models will be studied and the loannisian- 
Pilaftsis and Shtanov models will be analyzed in detail^. 

Finally let us mention again the interacting wave packet model recently proposed by Giunti 



|138|] which was described at the end of section 2.5, since it bridges the gap between the inter- 



mediate wave packet model and quantum field treatments. Moreover its results are equivalent 
to those obtained in the external wave packet model. 

3.2 Five questions 

The contradictions between the existing quantum field derivations of the oscillation formula can 
be summarized into five questions: 

1. Whereas the external wave packet model allows to associate intermediate 'wave packets' 
to the oscillating particle, it is not possible to do so in models using stationary boundary 
conditions. Moreover, the coherence length is finite in the former case, but infinite in the 
latter. Is it possible to see the oscillation formula derived in the models using stationary 



boundary conditions |18?, 15C, 151, 152, 163 , 87], as a particular case of the oscillation 
formula derived in the models using external wave packets [145, 144, 35 1? 



2. Kiers, Nussinov and Weiss [182] have shown in a quantum-mechanical model that the value 



of the coherence length depends on the accuracy of the energy-momentum measurements 
at the detector. This effect has been confirmed by quantum field theory calculations 



Some of our results have been published in Ref. @. 
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1 144 , 185 ], but it is not clear whether oscihations survive a perfect measurement of the 



energy [144, 139 1. Does a perfect knowledge of the energy lead to an infinite coherence 
length or is there decoherence anyway? 

3. For a source strongly localized in space and time, Shtanov has derived an oscillation formula 
where each oscillating exponential exp{—i4>j) is multiplied by a prefactor depending on the 
mass rrij. Unless the masses nij are nearly degenerate, these prefactors strongly modify 
the standard oscillation formula given by Eqs. (|l^)-([T^. Is this result correct? 



4. loannisian and Pilaftsis |163| ] claim to have found a novel form of neutrino oscillations 
('plane wave oscillations'), if the spatial spread of the source and detector is of macroscopic 
size. Does such behavior exist? 



5. As in the quantum- mechanical treatment of section |2.3| , the correlation between an oscil- 
lating particle and other particles (or recoil particles) at the source has been said to modify 
the oscillation length with respect to Eq. (|T^ ) and to bring about oscillations of the recoil 



particles |245, 247]. Is this assertion true? 



In order to answer these questions, we shall use a model such that all specific models used in 
the articles cited above can be recovered in some limit. Moreover, to extend the usual treatment 
of a stable relativistic oscillating particle, our model will also allow the oscillating particle to be 
nonrelativistic and/or unstable and will take into account the dispersion. Answers to questions 1 



to 4 can be found in section 7.3 while question 5 is treated in section IC. 



4 The external wave packet model 

The numerous problems arising in the plane wave approach of oscillations (see section |2.4D and in 



the intermediate wave packet method (see section 2^) show that the oscillating particle cannot 
be treated in isolation. The oscillation process must be considered globally: the oscillating 
states become intermediate states, not directly observed, which propagate between a source 
and a detector. This idea is easily implemented in quantum field theory, where intermediate 
states are represented by internal lines of Feynman diagrams. Quantum field theory has the 
advantage of providing a relativistic treatment from the start, which is required to study the 
mixings of relativistic and nonrelativistic particles. It also allows to describe unstable particles 
in a consistent way. The oscillating particle is represented by a relativistic propagator, which 
determines the space-time evolution and the possible decay of the particle. Boosts and the 
consideration of specific rest frames become pointless. Since interactions are included in the 
amplitude, equal or different time prescriptions are not needed anymore. Particles interacting 
with the oscillating particles at the source and at the detector are described by wave packets. 



This external wave packet model was first proposed by Sachs [238] for kaons and by Giunti, Kim, 
Lee and Lee [11451] for neutrinos. 



28 



As appears clearly from the review of section |3.1| , there is an ongoing controversy in the 
literature about whether the boundary conditions (i.e. the source and the detector) should 
be taken as stationary or not. In other words, can the time independence of most oscillation 
experiments be translated into the assumption that the energy of the process is perfectly defined? 
Although it is true that most sources are stationary from a macroscopic point of view (for 
example, the flux of solar neutrinos is steady), there is no reason to think that it should be 
the case from a microscopic point of view. Whereas it can be reasonable to make a stationary 
approximation in a quantum-mechanical model where we do not have any information on the 
microscopic processes (for example a density matrix model as in p48| ] ) , the same approximation 
is very dubious in quantum field theory, which describes the interactions of individual particles. 
After all, perturbative quantum field theory applies to one-particle propagation processes, in 
which a stationary source and detector are the exception, rather than the rule. As emphasized 
by Cardall |85], the Sun is certainly not stationary at the atomic scale, and neither is a detector 
composed of bound state particles. For example, a water Cerenkov detector sees charged lepton 
wave packets with finite energy and time spread. The finite character of the spread is partly due 
to the limited coherence time of the bound state particle that has interacted with the incoming 
neutrino. Boundary conditions can be considered as stationary at the microscopic level when the 
energy uncertainty at the source or detector is smaller than the inverse propagation distance, i.e. 
(^ep,D ^^/L. This constraint is extremely stringent. In the example of atmospheric neutrinos, 
the process can be considered as stationary if aep^o ^ 10~^^MeV, which is not satisfied in 
current oscillation experiments. 

Although it seems difficult to argue that the energy uncertainty is smaller than the mass 
difference, it is possible to take up another stand regarding stationary boundary conditions. In 
the same way as in section 2^ , it consists in arguing that interference occurs between wave packet 
components with the same energy, because of the time average on the transition probability. As 
it is unrealistic to consider in isolation an interference for a given energy, we should take care to 
integrate the probability (computed with stationary boundary conditions) over the wave packet 
energy width. At that point a question arises: is this energy width determined by the source 
or by the detector? There is thus information to be gained by working with non-stationary 
boundary conditions. Moreover, we shall show that stationary boundary conditions can be 
imposed in the external wave packet model by assigning zero velocities to some states, whereas 
other states are represented by plane waves. 



In section 4.1, we describe the Jacob-sachs model which is the prototype for the external 



wave packet model. We then examine the diagonalization of mixed propagators in section 4.2. 
Next, we compute the transition amplitude of the process associated to particle oscillations in 
section 4.3. In section 4.4, we discuss the simplest model in which oscillations can be consistently 



described. Finally, we present in section 4.5 the external wave packet model in its Gaussian 
version. 
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Figure 3: Propagation of a particle between a source and a detector, centered in Xp and x 



4.1 The Jacob-Sachs model 



The first version of the external wave packet model was developed by Jacob and Sachs [166] to 
describe the propagation of a non-oscillating unstable particle. 

4.1.1 The process 

The propagating process of a particle between a source and a detector (extended in space-time 
and indicated by dotted circles) is symbolized by Fig. 3. 

The arrows indicate the momentum flow. Pj represents the set of incoming particles, of 
total momentum q, arriving in the production region (or source), which is centered around the 
point (tp,Xp). Pp represents the set of outgoing particles, of total momentum k, coming from 
the production region, with the exception of the intermediate particle whose propagation is 
studied. Dj, Dp and (tojXo) are defined similarly, but apply to the detection process. If the 
intermediate particle decays, Dj(q') can be seen as an outgoing state. The interaction points at 
production and detection are noted x and x' , respectively. 

The external particles are assumed to be stable, whereas the intermediate particle is allowed 
to be unstable. Although it is possible to describe external unstable particles by considering 
a more global process where all unstable particles are considered as intermediate states, the 
technical difficulties involved in computing the corresponding amplitudes are enormous. All 
computations will be carried out for a scalar propagating particle since it can be argued (see 



section 4.3) that the spin has no impact on the characteristics of oscillations. 

With the aim of localizing the production region around {tp,Xp), the in- and outgoing 
particles at point x are represented by wave packets which overlap (in configuration space) 
around (tp,Xp). As the energy-momenta are often measured, the wave packets are supposed to 
be well-localized in momentum space, around their average momentum. The detection region is 
localized in the same way. Note that there is no difficulty in replacing certain wave packets by 
bound states or by plane waves, if needed. For example, it is practical to use plane waves if an 



emitted particle is not observed [150]. 



If the dotted circles are removed, the above picture can be interpreted as a Feynman di- 
agram. The internal line represents either a particle or an antiparticle. The experimental 
conditions are chosen such that a quasi-real particle propagates on a macroscopic distance. 
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transferring positive energy from x to x'. The energy- momenta of the initial and final states are 
such that the S-matrix element, corresponding to the above process, is evaluated at the singu- 
larity of the propagator of the intermediate particle. If x'^ — x*^ is a macroscopic time, it will be 
seen that the intermediate particle (positive energy component) contributes to the propagation, 
but not the intermediate antiparticle (negative energy component). 



4.1.2 Wave packets 

A wave packet representing a given state | ^> of mass m can be expressed in momentum space 
|22|] as 

I V>= j [(ik]V'(k)|k>, 

where V'(k) is the wave function in momentum space taken at time t = 0, and | k> belongs to 
the basis of one-particle states of momentum k, also taken at time t = 0. In the interaction 
picture, the time-dependence of the basis and of the wave function cancel each other so that the 
state I "0 > is independent of time. That is why we may take tp(k) and | k> at time t = 0. The 
following notation is used 



(27r)3 ' 

where -E'(k) = \/k^ -|- m^. The normalization of free states is chosen as follows [226]: 

<k|p>= 2E{k) {27rf 6^^\p-k) , 

so that 

<V|^>=1 if I ^|^(k)|2 = l. (20) 
In configuration space, the wave function reads 

It is here necessary to use the time-dependent wave function, i.e. tp{k) exp{—iE{k)t), since the 
wave function is considered alone. 

If the wave packet represents a particle with an approximately known momentum K, 
the wave function V'(k) is sharply peaked at k = K. The configuration space wave function 
'i/'(x, t = 0) has a maximum at the point x = if tpCK+k') = tpCK — k'). On this assumption, the 
wave function will be noted ilj{k, K). 

Wave packets centered in xg at time to built with the help of the space-time translation 
operator exp(zP-xo), where xo = {tQ,xo). If a wave packet in momentum space is given by 

^^(k,K,xo,to) = V'(k,K)e^^«*o-*^-''o, (22) 

the corresponding wave packet in configuration space will be peaked at the point xq at time to^ 

$(x,t,K,xo,to)= |^^(k,K)e-^('^)(*-*»)+^'^-(™). 
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Without loss of generality, let us choose to work with only one particle in Pi{q), in Ppik), 
in Di(q') and in Dpik'). The extension to a larger number of initial and final particles is 
straightforward and would only complicate the notation. The wave packets are built such that 
those involved in the production are centered at Xp at time tp, whereas those involved in the 
detection are centered at Xp, at time tp,. They are noted | 145 , 144 1 

\Pl> = yM^-p, (q,Q,Xp,tp)|P/(q)> 

\Pf> = j[dk]^p^{k,K,Xp,tp)\PF{k)> 

\Di> = j[dq']^fDj{q',Q',^n,tn)\Dj{q')> 

\Df> = y[(ik']^'z5p(k',K',Xp,tp)|Z)p(k')> . 



Gaussian wave packets are discussed in section 4.5.1. 



4.1.3 The propagation amplitude 

The general formula of the connected amplitude corresponding to Fig. 3 is given by 
A=<PF,DF\f (exp(-i I d^xHi]] -l\Pi,Di>, 



where TCi is the interaction Lagrangian for the intermediate particle and T is the time ordering 
operator. Let g be the coupling constant in Tii. Expanding the amplitude to order g'^, and 
inserting the wave packets expressions, we obtain 

A = J [dq] ^p, J [dk] ^*p^ J [dq'] ^D, I [dk'] Api anewave 

{q,k,q',k') (23) 



with 



Apianen^aveiq, K q\ k') = J d^x Mp{q, k) e~'^^~^)- J d'x' M^(g', k') e-'^"' G{x' - x) (24) 

where Mp{q, k) and MF){q', k') are the interaction amplitudes at production and detection. The 
propagator G{x' — x) has been obtained by field contraction and reads 

" = / (0^ e-^"-^'"'-^ G{p') . (25) 

where G{p'^) = i{p'^ — Mq +ie)~^ is the free propagator in momentum space and Mq is the bare 
mass of the propagating particle. It is assumed that renormalization has been carried out so 
that Mq can be replaced by the physical mass m. The external particles are on shell: 



qO = EpM) = Jq'+nil. 



"Pi 

and so on. 
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We have supposed that the particle (p" > 0) propagates from x to x' and the antiparticle 
{p^ <0) from x' to X. If interactions at the source and at the detector are such that the particle 
propagates from x' to x and the antiparticle from x to x', the contraction of the fields would 
have yielded the propagator 

^^"^'"^^^ / (^^"''■^'''"''^ 

In that case, the different phase sign would select the antiparticle pole in the subsequent contour 
integral. 

If the propagating particle is unstable, the complete propagator should replace the free 
propagator in the amplitude. For this purpose, the amplitude is expanded to all orders and one 
sums over all diagrams corresponding to one-particle- irreducible insertions (IPI self-energy) in 
the propagator | }78| , 226]. This sum is a geometrical series in the IPI self-energy — in(p^) and 
yields the same formula as above, except that G(p^) represents now the complete propagator: 



It is convenient to proceed to a change of variables in Eq. p^): 

X ^ x + Xp and 
where Xp = {tp,'x.p) and x^ = {tp,,:x.p,). The amplitude (|23| ) becomes 

-4 = / P) G(/) e-P-(-«--^) , (27) 

where the overlap function ip{p^,p) represents the overlap of the incoming and outgoing wave 
packets, both at the source and at the detector. It is defined by 

i^ip', P) = I e^^- I d^x' e-'P-' I [dq] ^^p^^ (q, Q) e'^'^^ j [dk] ^p^^^ (k, K) e^^^ 

X j [dc^]i,,Jc^Mle-^'''-' l[dk']rn^Jk',K')e''^'^' Mp{q,k)Mn{q',k'). (28) 

All external particles are on-shell. Note that the overlap function is independent of Xp and Xq. 
The integrals over x and x' in Eq. (^) yield delta functions, which impose energy-momentum 
conservation at the source and the detector. Remark also that the overlap function depends on 
the directions of incoming and outgoing momenta. 

Most computations in field theory apply to microscopic processes, where in- and outgoing 
states can be approximated by plane waves. In that case, the wave functions V'P/, V'Pf' V'-D/ and 
ipDp become delta functions, so that the overlap function ip{p^,p) is proportional to 

V'(/, p) ~ 5^*^ (K + K'-Q- Q') {p-Q + K) Mp{Q, K) Md{Q', K') , 



where Q = ^yQ^ -|-mp^,Qj and so on. The momentum integral in Eq. ( p7|) is then easy to 
perform and yields 



A ~ {K + K' -Q-Q')G {{Q - Kf) e"' 



{Q-K)-{xo-xp) 
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This expression cannot lead to oscillations in a superposition of several amplitudes, since the 
phase of the exponential does not depend on the mass of the intermediate state. Thus, oscilla- 
tions vanish in the plane wave limit. 

In the case of an amplitude corresponding to a macroscopic process, such as the propagation 
on a macroscopic distance, the overlap function is not a delta function so that the momentum 



integral in Eq. (27) is not trivial (see section 4.4 for a simple exact model). The main contribution 
comes from the pole associated to the propagating particle, so that the phase of the exponential 
will depend on the mass. 

4.2 Mixed propagators 



In section 2T, mixing was defined in quantum mechanics as the result of the noncoincidence of 
the flavor basis and the mass basis. The definition of mixing is similar in field theory, except 
that it applies to fields, not to physical states. This difference allows to bypass the problems 
arising in the definition of flavor and mass bases. 

The complete Lagrangian is split in a propagation Lag rangian Lpropag and an interaction 
Lagrangian Lint- These two subsets of the Lagrangian can be distinguished if there is a flavor 
transformation which is a symmetry of Lint but not of Lpropag- Particle mixing occurs if the 
propagator built from Lpropag, and representing the creation of a particle of flavor a at point 
X and the annihilation of a particle of flavor /? at point x' , is not diagonal, i.e. not zero for 
a ^ (3. The Lagrangian Lpropag always includes the kinetic and the mass terms. If the particle 
is unstable, Lpropag also contains the interaction causing the decay. 

In the simplest case, illustrated by the light neutrinos Ve, Vfj, and Ur, the propagating particle 
is stable and the form of its fundamental Lagrangian is known (although the mass matrix might 
have its origin in an effective Lagrangian). Lpropag contains the mass matrix, generated by 
Yukawa interactions, whereas Lint includes the weak interactions. In the flavor basis (called the 
gauge eigenstates basis for neutrinos), the mixed propagator is nondiagonal. It is well-known 



1 205, 64 1 that Lpropag can be diagonalized by a unitary transformation on the fields^: 

^" = E^i^i' (29) 

3 

where Greek and Latin indices refer to the flavor and to the mass basis, respectively. It could 
however turn out that only a subset of V is used to describe the mixing in an oscillation 
experiment. One reason could be that some of the mixed states are too heavy to be produced 
(see Eq. (^)). In this case, the mixing matrix for the light states can be nonunitary |63, 22C, p^ . 



Another reason could be that some of the mixed states are sterile, i.e. they only interact through 
mass mixing. In that case, the mixing matrix for the non-sterile states is rectangular and does 



not satisfy unitarity [24C] 



The propagator is defined by the time-ordered two-point function: 

Gpo.{x' -x)=<Q\T{vp{x')vl{x))\Q> . 



^The matrix V corresponds to the matrix 17* of section 
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Since Wick's theorem applies to fields in the mass basis, we must substitute Eq. (^) in this 
equation. In the scalar case, field contraction yields 

Gpa{x - a^) = X] GD,jj{x - x) Vja , (30) 

i 

where G^jjix' — x) is the free propagator (with D for 'diagonal') of a scalar particle of mass 



GD,jj{x' - x) = I -^-^ e 



-ip{x'—x) 



(27r)^ — rrij + ie 

The derivation of the mixed propagator for a spin one-half particle is similar. 

The treatment of unstable particles is more involved. The complete Lagrangian is again 
split in Lint and Lp^opag, but the interactions causing the decay should be included in Lpropag- 
The free propagator is then replaced by the complete propagator, obtained by an infinite sum 
on the self-energy — in(p^) p^ ]: 

iG~\p^)=p^l-M^-Uip^), (31) 

in an obvious matrix notation. The propagation eigenstates are characterized by the complex 
poles of the complete propagator Zj = mj — irujTj or, equivalently, by the zeros of the inverse 
propagator. It is always possible to diagonalize the propagator if its eigenvalues are distinct: 

G(/) = y-i(/)Gz5(p2)F(p2), (32) 
where G^ip^) is the diagonal matrix the elements of which are given by 

— nij + iniji j 

The matrix V is not necessarily unitary, since the self-energy matrix for unstable particles is 
usually not normal (see section Moreover, it depends on the energy. If the mass eigenstates 
are nearly degenerate, the self-energy can be approximated by its value at = and the 
matrix V becomes constant. In section |8.1| , we show that the energy dependence of the self- 
energy generates corrections to the amplitude in inverse powers of the propagation distance L. 
These corrections are due to production thresholds of multi-particles states [241, |166|, 78 1. 



4.3 The flavor-mixing amplitude 

In order to derive the flavor-mixing amplitude, the Jacob-Sachs model of section |4.1| has to be 
modified to take into account the different flavors. The process corresponding to the production 
at the source of a particle of flavor a, and the detection of a particle of flavor /? at a detector 
placed at a macroscopic distance from the source, can still be symbolized by the Fig. 3 of section 
[4.1.1| , except that the intermediate particle should be considered as a superposition of different 
mass eigenstates. The initial flavor a is tagged, for example, by the outgoing state Ppik), 
whereas the final fiavor /? can be tagged by the outgoing state Dpik'). If it is impossible to 
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identify the flavor at the detector (ex: K^,K^ tt+tt""), one should sum over the different 
flavors. 

If the energy dependence of the matrix V diagonalizing the propagator can be neglected 
(see Eq. (^)), the amplitude corresponding to the global process can be expressed as a linear 
combination of amplitudes Aj corresponding to the propagation of different mass eigenstates: 

A{a^(3,T,L) = Y, V^' Aj Vj^ , (34) 

j 

where the average propagation time T is defined hy T = — and the average propagation 
distance by L = x^, — Xp. The partial amplitude Aj has the same form as the propagation 
amplitude of an isolated particle, given in Eq. (27): 



^^ = jW)'' ^""'^^"'^ e-^^^+^P-^ . (35) 

The overlap function '(/'(p'^, p) is defined by Eq. ( p8| ) and the propagator for the jth mass eigen- 



state by Eq. (33) 



As noted in section [4.1.1| , it will be enough for our purpose to work with a scalar oscillating 
particle. When the spin is taken into account, three cases must be distinguished: 

1 . If the particle is relativistic and the interactions at the source and detector are chiral, only 
one helicity eigenstate contributes to the propagation. The spin structure can then be 
factorized and the computation can proceed with a scalar propagator. The two mentioned 
requirements are satisfied for light neutrino oscillations, since the energy threshold of the 
detectors is much higher than the light neutrino mass scale, and because of the chirality 
of the Standard Model flavor-changing interactions. 

2. If the particle is nonrelativistic, the spin structure is approximately equal for the different 
helicities provided that the mass eigenstates are nearly degenerate: rrii = rrij = m. The 
helicity factors can then be factorized from the sum over the mass eigenstates: 



^r",(T,L) ~ \y\H{s,m)\^ \A%^iW{TM 



where H{s, m) includes all helicity dependent factors (with the index s referring to the spins 
of the external particles). However the source and detector contributions to \H{s^ m)p 
cannot be disentangled if the particle is not relativistic, so that the whole process cannot 
be factorized into a product of source/propagation/detection probabilities (this point was 
emphasized Refs. |145,|87]). As the mass eigenstates are nearly degenerate, it seems natural 



to define the oscillation probability as the modification introduced by the mass difference: 



'^0,^/3(^5 '^)\mi=mj 

so that the calculation of the oscillation can be done as if the oscillating particle were 
scalar. The discussion is very similar for a relativistic particle with nonchiral interactions. 
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3. If the particle is nonrelativistic and the mass eigenstates have very different masses, each 
specific process has to be computed separately, with the influences of the amplitudes of 
production and detection carefully taken into account for each mass eigenstate. This situa- 
tion is not really relevant to oscillations, but rather to incoherent mixing, since oscillations 
are averaged to zero in such experimental conditions. In other words, the contributions of 
the Feynman diagrams corresponding to the different mass eigenstates are summed in the 
probability, not in the amplitude. 

In this report, we always assume the condition (5m?j <C -E^, under which oscillations are observ- 
able, so that the third case is not considered. Some explicit examples of the third case have 
been given by Giunti, Kim and Lee |142| ], and by Kiers and Weiss |183]. The conditions of 



factorization of the transition probability have been discussed by Cardall and Chung [^]. 
4.4 The simplest consistent model 

A simple model, in which the flavor-changing amplitude (Eqs. (^)-(|35|)) can be computed 
exactly, has been proposed by Kobzarev, Martemyanov, Okun and Shchepkin |187]. In this 



model, a charged lepton, represented by a plane wave, collides with an infinitely heavy nucleus, 
situated in Xp. The neutrino produced at that point has a definite energy equal to the energy 
of the incident lepton. At point x^,, the neutrino collides with a second infinitely heavy nucleus 
and, as a result of this collision, a charged lepton is emitted with an energy equal to the neutrino 
energy. 



Let us see what these assumptions mean in the wave packet notation of section 4.1.2 . There 
are two ingoing states at the source, the lepton and the nucleus. The lepton is represented by a 
plane wave, so that 

The uncertainty on the momentum of the nucleus is infinite, so that the ingoing and outgoing 
momentum wave functions are constant: 

'^Pi,nucii^n, Qn) ~ const and VPF,„„d(kn, K„) ~ const . 

The situation is similar at the detector, except that there are one ingoing and two outgoing 
states. 

The overlap function ( p8| ) is easily computed: 

^Pip^, p) ~ - E,n) <5(pO - Eout) , (36) 

where Ein is the energy of the incoming lepton at the source and Eout is the energy of the 
outgoing lepton at the detector. The partial amplitude (|35|) becomes 



Aj ~ S{E,n - Eout) J d^P GD,n{E,n,p) e''' ^ 
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where L = x^, — Xp. The Fourier transform of the propagator with respect to its momentum 
can be computed by contour integration and yields 



A- ~ X 5{E,r, - Eout) e'P^'^ , (37) 



where pj = \J Efn ~ ""^i ^'^d L = |L|. 

The transition probabihty between an initial state of flavor a and a final state of flavor (3 
is given by the squared modulus of the amplitude ( |34|) with Aj given by Eq. (^7|). Oscillations 
between mass eigenstates rrii and rrij arise from interference terms AiA*^ ~ 6*^^*"^^-*^. Thus, the 
zj-interference term oscillates with a frequency equal to 

{p,-pj)L^-^L, (38) 



where Smfj = mf — m| and pm = \J ~ with m the mass in the degenerate limit. The 
phase ( |3^ ) is equal to the standard oscillation phase ( [T^ ) calculated in the quantum-mechanical 
plane wave treatment. 

This very simple model has the advantage of being consistent, unlike the plane wave ap- 
proach in quantum mechanics. It could be recommended as a pedagogical tool. It is however not 
sufficient for a thorough study of the oscillation phenomenon. First, it cannot describe unstable 
particles, since they cannot decay into an infinitely heavy state. Moreover, the approximation of 
plane waves and stationary states is too strong to allow the study of the observability conditions 
of oscillations. Another drawback of this model is that the amplitude is independent of the 
direction of L. Last, a spatial localization of the source or the detector more precise than the 
Compton wave length of the lightest external particle is unphysical |165| ] . 

4.5 The overlap function for Gaussian wave packets 

In most cases, it is not possible to compute exactly the flavor-changing amplitude (p^). The 
technical difficulty lies in the energy- momentum integration in Eq. (|35|). The two exceptions are, 
on the one hand, the Kobzarev et al. model of section |4.4| and, and the other hand, the limit in 
which all the external wave packets are plane waves (in which case oscillations vanish as shown 
in section 4.1.3| ). To proceed further, we have to be more specific about the shape of the overlap 



function (^8|). In order to answer the questions on the coherence length, it is useful to work 
with an overlap function which depends explicitly on the energy and 3-momentum uncertainties 
at the source and at the detector. This can be done by approximating the in- and outgoing 



particles with Gaussian wave packets, as first proposed by Giunti, Kim, Lee and Lee |145, 144]. 
The general case of arbitrary wave packets can then be analyzed as a superposition of Gaussian 
wave packets. Note that the Gaussian approximation is used by most authors since it allows 
analytical integrations. 
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4.5.1 Gaussian wave packets 



General wave packets were defined in section 4.1.2. The wave function of Gaussian wave packet 
can be written at time t = in momentum space as 

We chose to set the initial phase to zero. This wave packet is thus centered in x = at time 
t = (see Eq. (|2^)). Recall that the normalization is given by Eq. (20). The width ap^ is the 
momentum uncertainty of the wave packet, as can be checked by computing the variance of the 
operator p"". 

<{6p^f>^=<x\{p^ -Plf\x>=^l^, 
with the average momentum given by 

Px =<x\p''\x>=Px- 

It is also useful to define a width in configuration space by Cpx^^^x = 1/2- 

The wave function is given in configuration space by Eq. (21). This integral is Gaussian if 
the energy is expanded to second order around the average momentum p^-: 

p' + rnl ^ + vx • (p - Px) + ^ ((p - Px)' - (v^ • (p - Px))') , (40) 



where Ey, = y P^ + aiid = p^/Ey.. The wave packet is thus given in configuration space 
by 

(27ra2^)-3/4 / (x-Vxt)S-i(x-Vxt)~ 



^x(x, t) = -^=1^ exp ^-E,t + .px • X - ^ ^ J , (41) 

where the matrix notation is implicit and with the matrix S defined by 

2ita?, 

S*:' = 5'y + (5*^-44)^. (42) 

The spatial uncertainty of the wave packet can be computed in configuration space with Eq. (|4l|): 

<i5x'')^>^=<x\ix^ - xlf\x>= [{Re^-'rT ' (43) 
where the average position is given by 

x« =<x\x''\x>=vp. 
If the z axis is chosen along Vx, we obtain 

[{ReE-YT = [{ReJ^-')-T = ^ + ^-^^ (44) 
[{Re^-'rT = 1 + ^^^^^- (45) 
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These equations show that the wave packet begins to spread in directions transverse to at t = 
E-^/2ap^, whereas the spreading along begins at the later time t = E'^/2m'^ap^. In this report, 
we refer to the time-dependence of Eq. ( |4^ ) as the dispersion of the wave packet. Although 
dispersion will be neglected in the external wave packets, we shall see that the propagation 
amplitude associated to the oscillating particle spreads in the same way. 

In order to study the stationary limit 0, it is of great interest to compute the energy 

uncertainty of the wave packet (as proposed in Ref. [|138|| ). The average energy E-^^ and the 
average squared energy should both be computed to order ap. In momentum space, the average 
energy is given by 

2 

E^ =< x\E\x >= E^ + i3- v2)^ . (46) 

The second term comes from the second order term in the expansion of the energy operator 
around the average momentum. Remarkably the average energy E^^ is different from E^ = 



p2 + . The average squared energy can be computed exactly: 

<x\E^\x>=El + 3al, (47) 
so that the squared energy uncertainty is given to order ap by 

< i5Ef >^=<x\{E - E^f\x>= 4 . (48) 
In configuration space, the same results are obtained only if the dispersion (i.e. time spread) 

.2 

p 



of the wave packet is taken into account. We compute to order (additional terms appear at 



higher order) the average energy and squared energy: 

i„ fdT. 



< x\E\x > = E^--Tr[ — ] , (49) 



4 \ dt 

iEy _ /5S\ 



<X\E'\X> = ^x + ^^^l-^Tr^^y (50) 

With the help of Eq. (||), we can check that Eqs. (||) and are equal to Eqs. (|6|) and (^), 
respectively (but only if dispersion is not neglected). Thus the squared energy uncertainty is the 
same in both representations as expected (similar computations for the momentum uncertainty 
show that its value is the same in momentum and configuration space even if dispersion is 
neglected) . 

Note that the energy uncertainty coincides at first order with the naive expectation: 



5E^ ^(Px + (5p)2+m2--^p2 +m2^vx-(5p. (51) 

It could be objected that the squared energy uncertainty does not vanish to order ap/E^ in 
the limit 0. However external wave packet models are always built with non-spreading 

wave packets, i.e. one sets S*-' = 6^^ . If this approximation is taken seriously, we can forget 
about the original expression in momentum space. In that case, the average energy and squared 
energy computed in configuration space are shifted to E^ and £"2 -|- 'v'^CTp respectively with 
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no correction of higher order in 0"^ (see Eqs.(|49|)-(|50|)), so that the energy uncertainty is stih 
given by Eq. (^). Thus the energy uncertainty of non-spreading wave packets vanishes in the 
stationary hmit — > to all orders in cip. 

This result is of great interest to study the stationary limit of the external wave packet 



model. For example, the bound state wave function used by Grimus and Stockinger |15C] 
coincides with the wave function of a non-spreading wave packet at rest, i.e. Eq. ( pl| ) with 
v-^ = and S = 1. Besides, plane waves can be obtained by taking ^ oo in Eq. (|4l|). 

4.5.2 Gaussian overlap function 

We now compute the overlap function (Est) with non-spreading external wave packets as done in 



Ref. |144]. The factors Mp{q, k) and M^iq' , k') multiplying the exponential vary slowly over the 
width of the wave packet and can be approximated by their value at the average momentum. 
They can be factorized outside the sum over the mass eigenstates since it is assumed that the 
neutrinos are either relativistic or nearly degenerate in mass. The wave packets ipp. and 



tpDi„ have the form given in Eq. (39). Expanding the energy to first order around the average 
momenta, the momentum integrations in Eq. ( p8|) give results similar to Eq. (|4^) but with E = 1 
and additional factors \/2E coming from the integration measure. For example we obtain 



/ 



[dq] Vp„ (q, Q) e-'^^ = iVp,„ exp ( -lEp^^ (Q)t + iQ . x - ^ J^'"^^' ) , (52) 



where Np.^ = (Ina'^p, ) ^/^(2£^p.^ (Q)) is a normalization constant. 

Doing the same for the other wave packets, we can write the overlap function as 

^(/,p) =iVVp(/,p)V';(p°,p), (53) 

with 

WpKV ,P) = J d xexp\^i[p - Ep)t-i{p-pp) ■:>i ^ j , 

where Ep = Ep.^ — Ep^^^, pp = Q — K. The function ipoip^^p) is defined in the same way, 
with the index P replaced by D, except for the energy- momentum which is defined so as to be 
positive: E^ = Ep,^^^ — Ep,.^ and po = K' — Q'. The constant N includes the normalization 
constants as well as the factors Mp^u evaluated at the maxima of the wave packets. 
A new width cjpp has been defined by CTppaxp = 1/2, with 

1 1 1 



+ 



0-2 0-2 0-2 

XP XPj^^ ^Pout 



Gpp can be interpreted as the momentum uncertainty at the source. The width (Jxp is mainly 
determined by the external particle with the smallest space width. This is expected since the 
production region depends on the overlap in space-time of the external wave packets. 
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The symbol Vp is defined by 



r, / Vp, Vp , 

I in I -* out 



cr^-p — ;5 h 



It can be interpreted as the velocity of the production region, approximately equal to the velocity 
of the particle with the smallest spatial spread (unless the velocities of the different in- and 
outgoing particles are very different). 

The symbol Sp, satisfying < Sp < 1, is defined by 



/ v2 v2 



0-2 0-2 



Recall in all above definitions that there might be more than one in- and outgoing state in the 
production and detection process. 

Integrating over x and t, we obtain 

^p{jp\ p) = ttV^p^ a-pi exp (-/p(/, p)) , (54) 

with 

;,(,».P)=(P^+(|!!Z^1Z(PZP£L^, (55) 

where 

(tIp = (Sp - v^) < . (56) 

The quantity Uep can be interpreted as the energy uncertainty at the source, or also as the inverse 
of the time of overlap of wave packets during the production process. Indeed, we can show that 



where the sum is over all wave packets involved in the production process. This sum is dominated 
by the term including the two smallest wave packets in configuration space (unless their velocities 
are nearly equal). If Uxi is the smallest width and ax2 the second smallest, we obtain 

cr ^ ~'^2| ^ 1 .g^x 
0-^2 rpoverlap ' ^ 

where Tp^^^''"'^ is defined as the duration of the production process. Thus, crgp can be interpreted 
as the energy uncertainty at the source, since it is proportional to the inverse of the time of 
overlap of the external wave packets at the source. The quantities axo, CFpo, Vd, ^d, ctcd, 
rpoveriap ^^^^^ similar definitions and properties. 

Note that stationary boundary conditions are recovered by setting Vp ^ = and aep,D = 0, 
with app,D different from zero. Besides, we shall need to know the stationary limit of \vp^\/aep,D- 
Since Vp is the velocity of the production region and a^p is the duration of the production 
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process, the ratio |vp|/o"ep is bounded by the size Sp of the macroscopic region of production. 
The argument is similar for the ratio |v£,|/(TeD- Thus we assume the following constraints: 



(58) 



eP.D 



These bounds are very conservative, since we shall see that stationary models such as those 
found in Refs. |15C, 163] are recovered by setting |vpo|/o"ep.D = 0. In the example of the 



Grimus-Stockinger model [150|, an initial stationary neutron (vp,^_„ = 0) decays into a stationary 
proton i^p^^^,pr = 0), a 'plane-wave' electron {axp^^^,ci = oo) and the intermediate antineutrino. 
At detection, the antineutrino collides with a stationary electron (v^,.^^ = 0) and the outgoing 
antineutrino and electron are represented as plane waves [(TxDout.v = (^xoout.ei — argued 
at the end of section 4.5.1 , all these stationary states are limiting cases of wave packets. 

As the propagation distance is macroscopic, only processes satisfying global conservation of 
energy-momentum have a nonnegligible probability of occurring. Since our aim is not to prove 
this well-known fact, we impose that 



Pp = Pa = po 



and 



Er 



En 



(59) 



but we still have Vp 7^ Vj,. These approximations allow to do expansions around po and Eq. An 
associated velocity can be defined by vq = Po/^o • 



5 Propagation amplitude and dispersion 
5.1 Introduction 

In section ^, we described the external wave packet model and we applied to flavor-mixing 
transitions. All the information relevant to the source and detector was included in what we 
called the overlap function (Eq. (^)). We computed this overlap function, on the one hand, in 
the pedagogical case of infinitely heavy source and detector (Eq. (^)) and, on the other hand, 
in the more general case of Gaussian external wave packets (Eqs. (|53[)-(p5|)). It is now possible 
to evaluate the transition amplitude. 

Although the flavor-mixing amplitude, as defined in Eqs. ( |3^ ) -(|35|), is the common starting 
point of most field-theoretical models found in the literature, the obtained oscillation formulas 
do not always agree. In section ^, several conflicting results were mentioned, namely stationary 
states versus non-stationary states, the possible existence of plane wave oscillations and of mass 
prefactors. Actually, it will be seen in sections |6| and |^ that these apparent contradictions vanish 
if we proceed to a careful evaluation of the integrals present in Eq. (^5|). In the present section, 
the possible methods of evaluating the amplitude (|3^) are compared. We shall see that the choice 
of the integration method depends on the propagation distance, so that the distance range can 
be divided in three regimes. We then evaluate the amplitude in the three scenarios. 
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5.2 Pole integration and wave packet correspondence 

Since the experimental conditions are such that the propagating particle is on-shell, the main 
contribution to the transition amplitude (|35| ) comes from the pole of the propagator. Although 
integrals on and p can, in principle, be done in any order, it is practical to integrate first on the 
energy, i.e. to do the pole integration on the energy, in order to make possible an interpretation 



of the amplitude in terms of propagating wave packets. In section 6.4, we shall discuss another 
method of integration (for stable oscillating particles) that has the advantages of being shorter 
and of clearing up the meaning of the stationary limit. However, this last method neither lends 
itself to a wave packet picture nor can be easily extended to unstable oscillating particles. 

The crucial step in the pole integration is the choice of an appropriate contour in the 
complex plane. The contour should be carefully chosen as the analytic continuation of most 
overlap functions diverges at infinity in the complex plane. The integration on the energy can 



be performed with the help of the Jacob-Sachs theorem \lQt\. This theorem is based on the 
assumption that the energy spectrum of all incident particles is limited to a finite range. Thus 
the overlap function tp{E, p) is distinct from zero only for = E'^ — p^ within certain bounds 
(with E > 0). On this interval, 'ilj{E,p) is taken to be infinitely differentiable. In that case, 
the Jacob-Sachs theorem says that the asymptotic value (T ^ oo) of the energy integral in 
Eq. (35) is given by its residue at the pole below the real axis. Thus, the evaluation of the 
partial amplitude Aj with the Jacob-Sachs theorem yields 

where 



Ej{p) = ^m2 + p2 (61) 
<Ai(p) = Ej{p)T-p-L (62) 

For Gaussian external wave packets, the overlap function is given by Eqs. (|5|)-(||). In 
principle, this function should be cut off outside the energy range determined by experimental 
conditions so as to satisfy the conditions of the Jacob-Sachs theorem. However these corrections 
are very small and will be neglected in the computations. The reader is referred to the Appendix 
for more details on the validity conditions of this theorem. 

It is interesting to note that the amplitude ( |60[ ) is mathematically equivalent to the ampli- 
tude obtained in the intermediate wave packet model |176(| , in which the mass eigenstates are 
directly represented by wave packets. The overlap function ip{Ej{p), p) corresponds to the wave 
function of the jth mass eigenstate. Thus it makes sense, in an external wave packet model, to 
talk about mass eigenstate 'wave packets' associated with the propagating particle. 

This wave packet picture brings up a problem, as the overlap function takes into account 
not only the properties of the source, but also of the detector. This is unusual for a wave 
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packet interpretation and not satisfying for causality. This point was recently clarified in an 



interesting article by Giunti |138| which was already mentioned in section 2.5. In this paper, 
Giunti proposes a sophisticated version of the intermediate wave packet model, in which the 
wave packet of the oscillating particle is explicitly computed with field-theoretical methods in 
terms of external wave packets. The intermediate wave packet depends only on the function 
ijjp{p^,p) (see Eq. (|5^)). The translation of this wave packet to the detection point yields 



a factor exp[—i(pj{p)]. At the detector, its interaction with external wave packets gives rise 



to the function ipi:,(p^,p). Remarkably, the total amplitude is exactly the same as Eq. (6C). 
However the propagating wave packet now depends only on the properties of the source through 
the function ipp^p^jp). The only drawbacks of this model are that it needs the equal time 
prescription and that it is not applicable to unstable oscillating particles. 

5.3 Three propagation regimes 

5.3.1 Laplace's and stationary phase methods 

Unfortunately, the integration over the 3-momentum in Eq. ( |60| ) cannot be done analytically. 
Resorting to the explicit form (|53|)-(|55[) of the overlap function valid for Gaussian external 
wave packets, we see that the integral (|60|) can be approximated by means of an asymptotic 



expansion for which two kinds of large parameters can be used. On the one hand, CTppj^ and 



c7ep,o are large parameters appearing in the overlap function (p^)-(|55|). They suggest a second 
order expansion of the integrand around the maximum pj of the overlap function, followed by 



a Gaussian integration: this is called Laplace's method [111, 53 1. On the other hand, T and 
L are large parameters appearing in the phase. They suggest a second order expansion of the 
integrand around the stationary point Pdj of the phase, followed by a Gaussian integration: 



this is called the method of stationary phase [111, p3[]. The competition between these two 



asymptotic behaviors implies a detailed study of the oscillation of the phase around the average 
momentum pj (Fig. 4). The expansion of the phase in Eq. ( |60| ) should be compared with the 
expansion of the overlap function. Although both methods are expected to lead roughly to the 
same answer in the case of the propagation of a single particle, it should be checked whether 
the delicate compensation mechanism resulting in the oscillation phase is independent of the 
method chosen. 

In section |5.3.2| , it will be shown that Laplace's method is preferable if the dispersion (i.e. 
spreading in time) of the amplitude is negligible. This is certainly the case sufficiently near the 



source but is not true anymore beyond some threshold. From the discussion of section [4.5.1 
it is expected that the spreading of the amplitude first begins in directions transversal to the 
propagation axis, before becoming significant along this axis. It is thus appropriate to divide 
the propagation range into three regimes: the no-dispersion regime, the transversal-dispersion 
regime and the longitudinal- dispersion regime. Laplace's method will be used in the first regime 
in all directions and in the second regime along the propagation axis, whereas the method of 
stationary phase will be used in the second regime along transversal directions and in the third 
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overlap function: i|) (Ej(p),p) 




Figure 4: Laplace's method is preferable if the phase varies slowly (upper drawing), whereas 
the method of stationary phase is preferable if the phase varies rapidly (lower drawing). The 
overlap function (resp. the cosine of the phase) is indicated on both drawings by a continuous 
(resp. dotted) line. 

regime in all directions. 

As the phase of an amplitude corresponding to a propagating particle oscillates in pL, the 
interference between two amplitudes corresponding to two mass eigenstates provides a phase of 
the form epL, where e = 5mfj/2p'^, with 5mfj = m? — m^. When the conditions of production 
and detection of the oscillating particle are studied with the help of Gaussian wave packets, a 
new parameter is introduced, namely cTp. Therefore, the parameter e can also appear with other 
dimensionless coefficients like {p/cTp)^ and (apL)'^. These new coefficients are smaller than the 
original coefficient pL if the distance L is macroscopic. The new contributions will be neglected 
in the phase, since they are much smaller than the oscillating phase epL. On the other hand, 
they should not be neglected in the modulus of the amplitude, where they are the leading terms 
in the e-expansion and yield observability conditions of oscillations. 

We shall thus calculate the transition probability to 0{e'^) in the real part of the argument 
of the exponential and to 0{e) in the phase. Gaussian integrations will be consistent with the e- 
expansion if the extremum is computed to 0{e), the real part of the argument of the exponential 
to O(e^), the phase to 0{e), the first derivatives to 0(e) and the second derivatives to 0{e^). 

5.3.2 No-dispersion regime 

As discussed above, the integration over the 3-momentum in Eq. (^) is done by Laplace's 
method if the phase varies slowly over the width of the overlap function, i.e. if the particle 
is sufficiently near the source. The aim of this section is to determine the extension of this 
no-dispersion regime and to compute the amplitude Aj in this range. 
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The first step consists in analyzing the form of the overlap function (|53|)-(|55D. The influence 
of the decay term 7j(p) (Eq. (63)) will be neglected in this section since we are near the source. 
The study of the overlap function amounts to study the argument of the exponential, i.e. of the 
function fj{p) defined by 

/,-(p) = ME,ip),p) + UEj{p),p) , (64) 

where fp^o are defined by Eq. (|55|) (with the constraint (^)) and Ej{p) by Eq. (|6l|). It is 
convenient to expand all results in small mass differences around (£"0, po, rno, vq), where 
ttiq = Eq — Pq and vq = Pq/Eq. The expansion parameter is noted e and refers collectively to 
all (5m^/2£'g, where drrij = mj — tuq. The value pj minimizing fj{p) is given to 0{e) by 

Pj = Po + (aup + P^d) tttt- , (65) 



where 



2En 

Vq - Vp.D 

20"r p n 



The associated energy Ej = \/Pj + "i^ can be expanded to 0{e) as 



Sm'i 

Ej = Eo + p—^, (66) 



where p = 1 + avQ ■ Up + /?vo • u^,. The values of the dimensionless coefficients a and /3 can be 
computed but their explicit expressions will not be needed. It is sufficient to know that p — > 
in the stationary limit. A velocity Vj = Pj/Ej is also defined for future use. 

We are now going to approximate the overlap function as a Gaussian and compute its three 
characteristic widths. At the extremum pj, the Hessian matrix of fj{p) reads to ©(e*^) 

yab ^ i_^!/z_|'n-^ 

§ab ^ ^ 

— 2 ~l~ UpUp + '^D^D 1 

where the momentum width dp is defined by 

1 _ 1 1 

9 9 ^ 9 ' 

(7 (7 (1 

u U p '^pjj 



(67) 



and is approximately equal to the smallest width among the production and detection momen- 
tum widths. The associated width ax in configuration space is defined by apax = 1/2. The 
matrix S"^ determines the range of p values for which the overlap function il){Ej{p),p) is not 
negligible. As S"'' is symmetric, it can be diagonalized by an orthogonal coordinate transforma- 
tion. The eigenvalues of S"^* are 

4^2 ' 



al^ = + J(4 + u2)±lj(u2 +u2)2_4(upXU^)2 



4^2 ' 2 ' 2 
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The eigenvector associated with o"^ is in the direction of Up x u^,, whereas the eigenvectors 
associated with a'^^. belong to the plane defined by Up and u^,. In the limit |up| » \uo\ (resp. 
I Up I <C |u£,|), the eigenvalues cr^ and ct^_ become degenerate and the eigenvector associated 
with (T^^ becomes aligned with Up (resp. u^,). This is also the case in the limit of parallel 
Up and Up,. These limits are relevant to the case of stationary boundary conditions which are 
examined below. 

Let us choose coordinate axes {ex,ey,ez) coinciding with the normalized eigenvectors asso- 
ciated with ((T^, (T^_, o"^_,_) respectively. The quantities {(7p,ap_,ap_^) (with ap±ax± = 1/2) can 
be interpreted as the momentum widths of the overlap function, since they give constraints on 
the range of p values for which the overlap function is non-negligible: 







< 


Up 




-Pj\ 


< 


Up- 






< 





The case of the stationary limit is of special interest since this assumption is used by several 



authors |187, |150| , 151 , p,63| , 87]. Recall that stationary boundary conditions are obtained in the 
external wave packet model by taking Vpp, — > and aep,D with |vp_p,| < (Tep^oSp^D (see 
Eq. (|58[)). In this limit, the velocity vq becomes aligned with the axis e^. 



,l^y ^ ^-.^ ^ , (70) 

and two eigenvalues become degenerate whereas the third diverges: 

2 2 

2 -'^ I 2 I 2 /ryi \ 

CT^+ — ^ + Up + U^ ^ DO , (71) 

In other words, the transversal widths (i.e. in directions orthogonal to po) are given by Up in the 
stationary limit, whereas the longitudinal width (i.e. in the direction of po) is given by 0"^+ 0. 

The second step consists in comparing the expansion of the overlap function with the 
expansion of the phase </'j(p) around pj, which reads 

,^,(p) - 4>Avj) + (v,r - L)(p - p,) + Pj)R''\p' - p]) , (72) 

where R^^ = 5°"^ — VqV^. The second derivatives have been evaluated to 0{e^). 

Laplace's method will be preferred to the method of stationary phase if the phase <pjip) 
varies slowly over the 'bump' of the overlap function. As Up > ap- > the variation of the 
phase will become important in the first place along the direction p^, then in the direction and 
finally in the direction p^. The criterion for the use of Laplace's method in all three directions 
pX,y,z T^iii thus be determined by considering the largest momentum width ap. The insertion of 
the constraints (^) into the phase ( [72|) yields first order conditions for a slowly varying phase, 

\v]T-L-\ap < 1, 
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b|r-L^|ap_ < 1, 

\v'^T-L^\ap+ < 1, (73) 



as well as a second order condition, 



T 



<^l<^, (74) 







where we have used the property q°'R°'^q^ < and the fact that <Tp is the largest width. For a 
given T satisfying Eq. ([74|), it is always possible to find a range of L values so that conditions 
(|7^ ) are satisfied. For other L values, the amplitude is negligible as will be checked on the result 
(see Eq. (^)). Thus the criterion allowing to choose the integration method is determined by 
Eq. ([7^): the integration on p will be done by Laplace's method ii T < Eo/ap or, equivalently 
(with L = |L| and po = |po|), if 

L<^, (75) 

since conditions (|73|) impose the relation L = vqT as long as (Tx>+ ^ L. 

For ax+ ^ L (stationary limit), we now show that condition (^) is directly obtained without 
going through condition (IT^). The overlap function imposes in that limit that |p| = -^JEq — m^, 
so that we are left with an angular integration with the angular part of the integrand given by 

Condition (|75|) shows that the angular variation of the phase in ( [76|) is slow with respect to the 
angular variation of the overlap function, in which case Laplace's method will give good results. 



Therefore condition (75) constitutes a good criterion for the use of Laplace's method whether 
the stationary limit is taken or not. 

Condition ( |75D is usually not verified in oscillation experiments, because L/ax ^ Po/cTp in 



most cases. This condition is the same than the one under which loannisian and Pilaftsis |16S] 
obtain 'plane wave' oscillations. 

If condition ( |75|) is satisfied, the evaluation of the integral ( |60|) can be done by Laplace's 
method and yields 



Aj = Napap_ap+ exp - iEjT + ipj ■ L - fjipj) - Fj{T) - 7i(Pi)j , (77) 
where A'" absorbs numerical constants. The function Fj(T) is defined by 

F,{T) = i(v,r - L)* + i^i?) (v,r - L) , (78) 

where S"* and R'^'' are considered as matrices. In the framework of the wave packet interpretation 
developed after Eq. (|6C|), the function exp(— Fj(T)) plays the part of the space-time envelope of 
the wave packet associated with the jth mass eigenstate. The elements of the matrix Re{T, + 
Z2^-R)~^ constrain the extent of the wave packet envelope in space-time. As T increases, the 
wave packet spreads because of the i-^R term. Thus the dispersion of the wave packet is due 
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to the second order term in the expansion of the phase (j)j{pj). Therefore, condition (|74| ) or, 
equivalently, condition (75) means that dispersion has not yet begun in any direction, transversal 
or longitudinal. For that reason, the range of L values defined hj L < po/a^ is called the no- 
dispersion regime. Note that this threshold was already met in the analysis of Gaussian wave 
packets in section 4.5.1 (see Eq. (^)). Of course, the above interpretation is not valid for 



ax+ ^ L, in which case the propagation time T becomes indeterminate and dispersion loses its 
meaning. 

Now that the origin of dispersion has been clarified, the term in Eq. (|7^ ) including R can 
be neglected with respect to S. Moreover, we choose to approximate Vj by vq in Fj{T). In 
comparison with the oscillation formulas that will be derived in the transversal- and longitudinal- 
dispersion regimes, this approximation will lead to the absence of the coherence-length term, 
since this term exclusively arises, when the dispersion is neglected, from the velocity difference 
Vj — Vj. Dropping the index j, the wave packet envelope in Eq. (^) can then be written in the 
coordinate system diagonalizing S as 

F(T) = ^ ° „ ' + ^^-^ ^ + ^^-^ , 79 

which shows that the conditions (^) assumed for Laplace's method are required to obtain a 
non-negligible amplitude Aj. In other words, po and L should be nearly parallel in order to 
have an amplitude significantly different from zero. We shall see that this constraint explicitly 
appears in the final oscillation formula (Eq. ( |111| )). 

5.3.3 Transversal-dispersion regime 



Let us now assume that condition (75) is not satisfied, i.e. L > po/cTp. In that case, Laplace's 
method cannot be used to integrate on all three components p^''^'^ in the amplitude ( |60| ) since 
dispersion becomes significant. However, the spreading of the amplitude is not identical in all 
directions. More specifically, the onset of dispersion in the direction po can be delayed by two 
factors. First, the matrix element R°'^ present in Eq. ( |72[ ) leads to a relativistic contraction 
(of 1 — Vq) in the direction po of the dispersion of the amplitude (see Eq. (|78|)). Second, the 
momentum width along po is given for aep,D ^ o'x (i.e. in the stationary limit (|7T|)) by a 
vanishing (Tp+. Thus Laplace's method is valid for a longer time T in the direction po than in 
directions transverse to this vector. For this reason, the choice of the integration method in the 
direction po will be postponed for a short while, whereas the method of stationary phase will 
be preferred for momentum integrations in directions transverse to po. 



Let the z axis be along L, i.e. L = L e^. As in section 5.3.2, the examination of the amplitude 
shows that the quick variation of the phase averages the amplitude to zero unless po and 
L are nearly parallel. The method of stationary phase can thus be applied in directions p^'f , 
the stationary points of which are given by p^ = pV = 0. The result of the method of stationary 
phase for the transverse momenta in (pOf) can be written as follows: 



A' = / dp exp {-i(l)j{p) - fj{p) - lj{p)) , (80) 
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where p = and 



<Pj{p) = ^p^ + m^T-pL, (81) 
fjip) = f,p{p) + fjoip) , (82) 



with 



2 



fAp) = '^'j^ + ^ (83) 

^"pp ep 

where Vp = v% and Po = Pq- The average energy Eq has been redefined so as to absorb a factor 
PqVp +PqV%. The definition of fjoip) is similar. Finally, 7j(p) expresses the possible decay of 
the oscillating particle: 

7.(P)= 7^'^ ■ (84) 



The function ^((l) comes from the transversal part of the overlap function. It expresses the 
geometrical constraint between the direction of observation 1 = L/L and the momentum po: 

«m--(-n^)- <-> 

where ap is defined by Eq. (|67|). The function g{l) restricts the propagation to a cone of axis po 
and angle arcsin{ap/po) (see Fig. 6 in section |7.lD . 

Numerical constants have been included in A^, as well as 4-volume factors appa^p and 
(TpD<^eD} s-iid an energy factor Eq. The constant ^ = Eo/la^ comes from the overlap function and 
acts as a cut-off for small T. As T is macroscopic, the prefactor 1/(T — i/i) can be approximated 
by 1/T so as to give a prefactor 1/L^ in the transition probability. This expected geometrical 
decrease is seen to originate in the transverse dispersion of the wave packet corresponding to 



the oscillating particle. As in section p.3.2| , it is convenient to define a reference mass mo and a 
velocity vq by 

ml = EI-pI, (86) 
vo = po/Eo. (87) 



As before, the choice of the method to perform the longitudinal momentum integration 
is done by comparing the expansions of the phase and of the overlap function around the value 
Pj for which fj{p) + 7j(p) is extremal (the decay term is not neglected in this section). The 
explicit value of pj will be computed in section ^is] for the stable case and in section for the 
unstable case. 

The expansions of fj (p) and {p) are given by 

m = /.(p.) + ^S^' 

m?T 

Hp) = Mpj) + i^jT -L){p- Pj) + -^{p - p,f , 
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where Ej = \JVj and Vj = Pj/Ej. 

The effective width (Jp^ff and (Txeff are defined by 



1 

(^peffCTxeff = ^ ■ (89) 

Note that o"pe// has the same form as the asymptotic value of cjp+ (see Eq. ([Til)). 

With the help of the wave packet correspondence discussed at the end of section 5^, the 
effective width can be interpreted as the energy-momentum width of the oscillation process, 
since it is the width of the overlap function. It is dominated by the smallest among the energy 
uncertainties (recall that agp.o ^ <^pp.d)- The effective width <Jxeff is then approximately equal 
either to the production or to the detection time uncertainty, depending on which one is the 
largest. 

Laplace's method is preferable if the phase (|88| ) varies slowly over the width o"pe//, i.e. if 
the two following conditions are satisfied 

\vjT -L\2apeff < 1, (90) 

^^^leff < 1- (91) 



As in section ^.3.2 , the first order constraint ( pO[ ) will be included in the result of Laplace's 



method. Thus the criterion allowing to choose between Laplace and stationary phase methods 
is given by Eq. (pT]). In other words, it is better to use Laplace's method if T is smaller than a 
dispersion time Tj^**^ defined by 

Tf'P = fo . (92) 
^"^ypeff 

The term 'dispersion time' is justified by the fact that it is the time at which the longitudinal 
dispersion of the amplitude becomes important, more precisely twice the initial size. A dispersion 
length L^*''^ can be defined by L^**^ = vqTj^^^. The distance range Po/o'p ^ L < -^qj i^q 

called the tranversal-dispersion regime. For L > j^^^^p ^ stationary phase method is more 
accurate: this distance range will be called the longitudinal- dispersion regime. Note that this 
threshold was already met in the analysis of Gaussian wave packets in section 4.5. 1| (see Eq. (p5[)). 



Various estimates of the dispersion length are discussed in section 5.4.2 , showing that the concept 
of dispersion length is relevant to nonrelativistic particles such as K and B mesons, as well as 
to supernova neutrinos, and possibly to solar neutrinos. 

In the transversal-dispersion regime, the evaluation of the amplitude (80) as a Gaussian 
integral around pj gives 



TJl + iT/Tj^ V 1 + iT/Tj ^ 4cj^.e// 



(93) 
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where absorbs numerical constants and pj is the value for which fj{p) + Jjip) is extremal. 
The amplitude (93) behaves as a wave packet of group velocity Vj and space-time extent 
axeff- If the longitudinal dispersion is neglected (T^^-'^f = 00), the amplitude 
is similar to Eq. (18) of Ref. |l44|. 



5.3.4 Longitudinal-dispersion regime 



At sufficiently large distance, dispersion becomes significant and all neutrinos propagating freely 



enter into the longitudinal-dispersion regime. In this regime (L > L'^ 
that the integral ([ 
point of the phase 



disp\ 



it has been argued 



) should be evaluated with the method of stationary phase. The stationary 
'j{p) is given by 

Vcl 



Pd,j 



(94) 



1 



where v^i = L/T. It can be interpreted as the classical momentum of a particle of mass rrij, 
travelling at the classical velocity Vd- Of course, a stationary point exists only for T > L. 
Otherwise Laplace's method must be used, but the amplitude is nearly zero in that case anyway. 
The evaluation of the amplitude (^0|) as a Gaussian integral around pdj gives 

/ 

exp -iruj y^T^ - - fj {p^j ) - 7j {Pd,j ) + 



^5'(l)f^pe// 



TJl + iT/T 



(lisp 



2 

Peff 



V 



1 + iT/Tl 



disp 



(95) 

where fj{p) is defined by Eq. (^), f'jip) refers to its derivative and N absorbs numerical con- 



stants. The wave packet interpretation of the amplitude (95) is not obvious but the shape of the 
associated wave packet can be studied by an expansion around the maximum of the amplitude. 



5.4 Estimates of the dispersion thresholds 

In section |5.3| , it was shown that the propagation range could be divided into three regimes 
separated by two thresholds (Fig. 5). The first threshold is determined by the no-dispersion 
condition (|75|) and separates the no-dispersion regime from the transversal-dispersion regime 
at L"'°'^^'^P = po/ap. The second threshold (see Eq. (p^)) is given by the dispersion length 
jl^disp _ pQj7;2^2mj<Tpgjj and separates the transversal-dispersion regime from the longitudinal- 
dispersion regime. It is not straightforward to estimate the value of these thresholds, mainly 
because of our ignorance of the external wave packet sizes. 



5.4.1 No-dispersion condition 

The first threshold is rarely of the order of a macroscopic distance. A far-fetched example where 
the propagation distance could be smaller than the no-dispersion threshold L"'"'^^^^ = po/<Jp is 
given by loannisian and Pilaftsis [163|: atmospheric neutrinos are detected in the no-dispersion 
regime if ax ^ 10^^ cm (with L~ 1000km and p^l GeV). Although such a macroscopic size of 
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Figure 5: The propagation range can be divided into three regimes according to the dispersion. 



the neutrino 'wave packet' is not totally excluded (see section 5.4.2| ), it is unlikely to be verified. 



However this case will be examined in detail in section 6.1 both for reason of completeness 



and because it was argued in Ref. [163] that the corresponding oscillation formula would differ 



significantly from the standard result. 
5.4.2 Dispersion length 

This threshold is a priori more interesting than the first, as it is of the order of a macroscopic 
distance. As said above, the main difficulty lies in the lack of knowledge of the external wave 
packet sizes. In the case of neutrinos, another uncertainty arises because of the dependence of the 
dispersion length on the absolute mass scale, instead of a dependence on the mass differences. 
Some examples are given below, partly based on the wave packet size discussion of section 
|9.2.2| . We should always remember that o"pe// is determined by the smallest width among the 
production and detection energy or momentum widths. 

First, the neutrinos in the experiment LSND ||28| , [29|] have a momentum around 30MeV. 
The neutrino mass can be taken to be 1 eV, and the production and detection widths should 
be around 0.01 MeV [152|. The corresponding dispersion length, around 10'' m, is much longer 
than the 30 m length scale of the experiment. Second, atmospheric neutrinos have a momentum 
around 1 GeV and let us take a mass of 1 eV as a bound. An estimate of the source width 
can be obtained from the inverse pion and muon lifetimes and is given by app ~ 10~^^ MeV 



'pp 

i-io 



in the case of a pion source. The detector width can be estimated at axo ~ 10 m [237] 



i.e. CTpo ~ 10 MeV, and is dominated by the source width. The corresponding dispersion 



length, larger than 10'^^ km, is completely irrelevant to the experiment. It is however doubtful 
that the source width is truly estimated by the inverse source lifetime (see the discussion of 



section 3.2.2). Third, solar neutrinos have a momentum around IMeV and let us take a mass 
of 1 eV as a bound. Various estimates of the source width have been given, between 10~^ MeV 
and 10^3 MeV |T|, |l8|, |200|, [l8|, |8|, |0|. The detector width can be estimated as above at 



Cpu 10^'^ MeV. The dispersion length, varying between 10^ km and lO^'' km for a 1 eV neutrino 
mass (depending on the chosen source width), could be relevant to the experiment, but is much 
larger for smaller neutrino masses. 

As regards mesons, kaons in the CPLEAR experiment have a momentum around 550 MeV 
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|120|]. If the width is guessed to be 0.01 MeV, the dispersion length, around 10~ m, is much 
smaher than the length scale of the experiment. More generally, the ratio between the dispersion 
length and the decay length (defined by L'^'^^'^v = p/mT) is given by 

j^disp Y 



Idecay (j^ „ 2m ' 
peff 

SO that dispersion occurs for all quasi-stable particles (F <^ m), and in particular for B mesons. 

In conclusion, if the estimates above are truthful, dispersion cannot be ignored for K and 
B mesons, and might be important for solar neutrinos (and certainly for supernova neutrinos). 
It is then necessary to develop a formalism taking into account the dispersion, in which case the 
use of the stationary phase method cannot be avoided. 

6 Oscillations of stable particles 

In section ^, the choice of the method of evaluation of the amplitude (35) was seen to depend 
on the distance range within which the oscillating particle detected. Three distance ranges were 
distinguished, depending on whether the amplitude had begun to spread in the transverse and 
longitudinal directions (with respect to the direction of propagation). The amplitude was then 
evaluated in the three possible cases, yielding formulas (|77|), (|93|) and (|95|) corresponding to the 
no-dispersion, transversal- and longitudinal-dispersion regimes respectively. 

In this sectionPI, we compute the transition probability for stable oscillating particles. The 
transition probability is then obtained by squaring the amplitude and averaging over the macro- 
scopic propagation time T, which is not observed in experiments: 

Va-.(3{^-Yl^^aVp-'V*,V^-/* J dT Ai A* . (96) 



We shall come back to the proportionality factor in sections ^ and 9T 



The oscillation formula is derived in the three different propagation regimes in sections 3.1 



6.2 and |6.3|. A shorter method of integration, valid for stable oscillating particles in the second 



and third regimes, is then discussed in section 6.4 



6.1 Oscillations without dispersion 
6.1.1 Amplitude 

For L ^ Po/Cp, the momentum integration in Eq. (|60|) is done by Laplace's method, yielding for 
stable oscillating particles (see Eq. ( [77[ ) with 7j(pj) = 0) the following amplitude: 




^''The reader not interested in all computational details can skip this section. Its results are summarized in 
section ItI. 
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where pj is given by Eq. (|65|), and 



{v?iT-L^f {vlT-LVf {v'T-L'f , , 



with o"a;, (Tx± defined by Eqs. (65). The parameter am, which has the dimension of a width, is 
related to fj{Pj) (see Eq. (|6^) by 

fjiPj) 



AamEo 



6.1.2 Probability 

The last step towards the oscillation formula consists in computing the time average of the 
transition probability, which is a Gaussian integral on T: 

where Tq is the solution of F'(Tq) = 0, while the phase (/)jj(To) is given by 

(^ij (To) = (Fi - Ej)To - (pi - Pj) • L . (100) 
Non-exponential factors are included in Ng. The function F(T) can be rewritten as 



(voT - L 



i2 



F(T) = ^ - ^ , (101) 



where 



Vo = cr^VS Wo, 
L = o-^\/s^L, 

with E^-*^ = diag{a^'^, cr^^, cr'?). It is then easy to compute 



To = (102) 



(vo X L) 



T"(To) = ^. (104) 

Actually, the time interval of integration AT is finite. The approximation AT ^ oo is good if 
AT is larger than the width of the overlap function given by Eq. ( |104| ): 

|vo|AT>a,e//, (105) 

where 

O'xeff = , , 

I Vol 



^xeff = CTx ■ (106) 



56 



The insertion of Eq. ( |103| ) in Eq. ( |99D shows that F{Tq) plays the role of a directional 
constraint since / dTAiA* is non-negligible if L is nearly parallel to vq- In order to express this 
constraint in terms of non-tilded quantities, two cases should be distinguished: ax+ ~ cr^ and 
ax+ > CTx- 

For ax+ ^ ax, J dTAiA* is non- negligible if L is nearly parallel to vq: 

L = ^L + 0{ax+). (107) 



't^ijiTo) = ^{L + 0{ax+)) , (108) 



With the substitution (107), the value of Tq given by Eq. (|102| ) becomes Tq = L/|vo| + 0{a, 
so that the phase ( 100 ) becomes 

which is equal to the standard oscillation phase (|l4|) if UxJr ^ L. 

If ax+ ^ (Tx, Vq ^-nd tend to zero so that F{To) should directly be studied as a function 
of Vq and L. General conclusions for arbitrary ax± can be drawn from the study of the quadratic 
form in (L^ , L"^ , L^) associated with F{Tq) = 1. This analysis shows that there is an eigenvalue 
S3 = corresponding to an eigenvector along vq. The two other eigenvalues si^2 are positive 
(with S2 < si) , so that the surface F{Tq) = 1 in (L^, L^, L^)-space is a cylinder of elliptical 
section with an axis along vg. This geometrical picture can be interpreted as imposing that the 
components of L orthogonal to vq should be smaller than y^l/s2, whereas there is no constraint 
at all on the component of L along vq. In other words, the probability is maximal within a 
cylinder of axis vq. 

In the stationary limit {crx+ oo), the lengthy expressions of the non-zero eigenvalues 



become 

si 



1 

4^ 



The properties of the stationary limit, namely Eqs. (|58|), ( [70| ) and (|7l|), lead to the bound 
■\f\fs2 ^ Sp^u- In the stationary limit, the components L^'^ (which are then orthogonal to vq) 
should thus be smaller than Sp^^. The time-averaged probability is therefore always negligible 
in directions other than along the average momentum po- Thus the evaluation of the phase 
(/)ij(ro) yields 

^^A'r^) = -^{F^o{Sp.^)), (110) 

which is equal to the standard oscillation phase (|l4|) if Sp p, <^ L. 

With the results (H), (|lO|), (|o|) and (|T08|) (or (|l|)), the flavor-mixing transition prob- 
ability ( p6| ) can be written for a propagation distance satisfying Sp^p, <^ L < po/ap as 



P„_^(L) ~ exp ( -i^g^ ] Vtp V;^ V, 
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X exp —2711 — tt^ 2tt — , (HI) 

y I j^osc Wa^E^ I rose I I ' V ; 



where p and &xeff are defined by Eqs. (36) and ( |106| ), respectively. Ng includes all pref actors 
independent of L as well as normalization factors. The oscillation length L°J^ is defined by 



As discussed above, the probability is maximal within a cylinder of axis vq and radius smaller 
than Sp^o. There is no decrease in since the spatial spreading of the 'wave packet' has not 
yet begun. 

6.2 Oscillations with transversal dispersion 
6.2.1 Amplitude 

For Po/(^p ^ L < L^*''^, the momentum integration in Eq. (|60| ) is done by the method of station- 
ary phase in the transverse directions and by Laplace's method in the longitudinal direction, 
yielding in the stable case (see Eq. ( |93[ ) with 'Jjipj) = 0) 



TJl + iT Tf'P V l + iT/T. ^ 4cj^e// / 



with fj{p) and o"pe// defined by Eqs. (|8^) and (p9[), respectively. We had not yet determined 
the value pj for which fj{p) is extremal. The computation of pj to 0{e) yields 

Pj=Po + {p-l)^, (114) 
where the dimensionless number p is defined by 

P = ''peff [ Z2 -2 -2 ) ' (115) 



with (Tp defined by Eq. (pTf ). The associated energy Ej = ypj + mj and velocity vj = Pj/Ej 
are given to 0(e) by 

Ej = Eo + p—^, (116) 
V, = -o + (.(l-.g)-l)^. (117) 



The parameter p has been defined so as to be in correspondence with the notation of Ref. |144]. 
The symbol oj appearing in that article is related to our notation by a; = ap/a^^jj:. Note that 



the authors of Ref. [144| do not compute p explicitly and also take the relativistic limit vq = 1. 
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The exphcit value of p is very interesting to know, since /? = in the case of stationary boundary 
conditions, in which case all mass eigenstates have the same energy Eq. 
The value of fj{p) at order is given by 

2 



1 / StdP' 



where am is defined by 



\ = <// {-A^ + ^]+ ^^^^1=^1 • (119) 
(Tm will be called the mass width, as it imposes a constraint between the masses mjj and the 



ingoing energy-momentum (see section 7.1.5). 



6.2.2 Probability 

As in section |6.1.2| , the last step towards the oscillation formula consists in computing the 
time average of the transition probability. Actually, the computation can be simplified by 



considering separately two cases. The examination of Eq. (118) shows indeed that the amplitude 
Aj (Eq. ( |113[) ) is negligible unless \6m?-/2EQ\ < a^- Together with the same condition on Ai, it 
puts a constraint on the mass difference: 

/ Smf- V 2 2 2 2 

I I ^ 2cT^ where Sm^j =171^- ruj . 

There are two possible cases: 

1. the masses are nearly degenerate; 

2. if not, the masses must be very small in comparison with the energy Eq, i.e. the oscillating 
particles are relativistic. 

We first consider the case of nearly degenerate masses. If |mj — mj\ <C mi,mj, it is possible 
to make the approximation Tf^'^^ = Tj^'^^ and to work with only one dispersion time defined by 
rpdisp _ £'3^2mgO"pgjj, with mo the mass in the degenerate limit. The integrand in ( |9^ ) can be 
written as 

with the phase given by 
^,,(r, L) = (E, - E,)T - ip, - p,)L - ^ ^ ^ (^^'^,,,,^2 ^'^^"^r,2"j'^'^"^^' ' (120) 



and the function fij{T,L) given by 



(v^T - Lf + (vjT - Lf 



fij{T,L)-^^ 4^2 • (121) 
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The prefactor A^^ absorbs numerical constants. 

The time integral in Eq. ( |96|) can be evaluated with Laplace's method. We want to compute 
the transition probability to order in the real part of argument of the exponential, and to 
order e in the phase. This again implies computing the minimum Tij of fij(T,L) to order e, 
fij{Tij,L) to order e^, its second derivative to order e*^, and the first derivative of the phase to 
order e. The minimum of fij{T, L) is reached to order e*^ for T = L/vq, and to order e for 

Tij = -(l- ^ ^ ^ + 0{e^) , (122) 



vq V 2^0 



If the expansion (117) is used, the minimum of fij{T,L) reads 

f- (T - L) - 



2 



- ^Tl^^(.a-*-if(|fj (1.3, 

where £ = L / (vqT'^^^'P) and 1 — = m^/EQ. The velocity vq has been replaced by vq in Vij 
because \p{vq — Vq)\ <^ 1. The value of the second derivative of fij{T,L) reads 

Y^(T,.L) = '-^ + 0(.). (124) 

The value of the phase at T = Tij is 



5m: 



2 

b,^{Tij,L) = —^L + Oie^). (125) 



The derivative of the phase is given to order e by 

2£ 



— {Tij,L) = E,- Ej - ^-p^ ap^ff (vi - vj) L 

+ 0(e2), (126) 



1 + ^2 2m^ 



''0 

where the expansions (|116| )- (|117D have been used. As above, the velocity vo has been replaced 
by Vq in Vij because \p{vq — Vq)\ <^ 1. The second derivative of the phase is of order e and thus 
does not contribute to the transition probability to order e^, at least in the argument of the 
exponential. 

The approximation of the time average integration by Laplace's method can now be done 
and yields 

^7^(1) 



J dTAiA* = voN'^apeff 



L2 
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where A^^ absorbs numerical constants. The oscillation length Lff'^ for the masses rrii and m,- is 
given by 



Without loss of generality, Smfj is taken to be positive. The coherence length L'^°^ is defined by 



L-^ = . (129) 

V2tT (Tpeff 



It is assumed that the time interval AT used to average is larger than the width of the overlap 
function: 

voAT > a^.ff (£ « 1) . (130) 

Moreover, the time separation T'^^'p between the wave packet peaks is supposed to be smaller 
than the time interval AT. This assumption is true as long as the distance is not hugely larger 
than the oscillation length, since Eq. (|122| ) shows that T*'^^ is of the order of 

Tsep ^^r^u _j,n^^ (131) 
L°j vopo 

where I/vqPq is of the order of a microscopic time. The violation of the condition AT > j^^^p 
gives rise to the interesting possibility of detecting separate pulses, corresponding to the different 
mass eigenstates |176, |236| ]. For example, a time-dependence of the neutrino burst from the 
supernova SN1987A has been searched for (see |]3^ , |234| ] for reviews), but only upper mass limits 
have been derived. 

Let us now consider relativistic particles with very different masses. Without loss of gen- 
erality, we suppose that nii ^> nij. As Tf^''^ <C Tj^"^^, the approximation of taking only one 
dispersion time for the two mass eigenstates is not valid anymore. Therefore, we shall suppose 
that L < j^dtsp show that the decoherence sets in before the dispersion length ^^'■^p ig 
reached. The calculation of the time average proceeds as in the nearly mass degenerate case, 
with the approximations £f 1 (with ii = L/hf''^^) and p ^ 7q (the particle is relativistic), 
so that the same result is obtained (see Eq. ( |127| )). Now, the coherence length Lf-^ is shorter 
than the dispersion length x"'*''^ if 5mf-/mf > apejf/EQ. This condition is always true for very 
different masses, since in that case 6mfj/mf = 1. Therefore, the interference becomes negligible 
before the dispersion of the heaviest mass eigenstate begins, so that the result (127) is also valid 
for very different masses. 

In conclusion, the flavor-mixing transition probability ( p^ ) for a propagation distance L = 
LI satisfying Po/cTp < L < min{Lf^^,L'^^^^) has the same form whatever the mass values and 
reads 
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6.3 Oscillations with longitudinal dispersion 
6.3.1 Amplitude 

For L > jJ^^'^P ^ the momentum integration in Eq. ( |60|) is done by the method of stationary phase 
in all directions, yielding in the stable case (see Eq. (^) with Jj{pd,j) = 0): 

/ 



Ng{\)apeff 
tJi + iT/Tf'^ 



exp 



Vt^ - - fj{Pci,j) + ct: 



2 

peff 



f'jiPdj) 



-idisp 



(133) 



v 



where fj{p) and CTpeff are defined by Eqs. ( |82D and 
given by Eq. (||): 

rrijVci _ 



1 + iT/Tf 
respectively. The stationary point is 



Pclj 



rrijL 



1 



The corresponding energy is given by 



VT2 - L2 



(134) 



VT2 - L2 

The value of fj{pd,j) is given by fj{pd,j) = fjp{Pd,j) + fjoiPdj) with 



(135) 



fjp.oiP, 



cl,j) 



{Pd,j - Pof , {Ed,j - Eq - {pclj - Po)Vp,dY 



40-2 



+ 



4a2 



(136) 



while the value of fj{pd,j) is given by fj{pd,j) = f'jpiPdj) + fjoiPdj), with 



6.3.2 Probability 



Pd,j - Po 
2rj2 

^'^pP,D 



Ed,j - Eq - {pclj - Po)Vp,D 



2^x2 



As in sections |6.1.2 and 5.2.2 , the last step towards the oscillation formula consists in computing 
the time average of the transition probability. Examination of the term fj{pd,j) in Eq. ( |133| ) 
shows that the amplitude Aj is nearly zero unless \pd,j — Pol ^ fpp.o (see Eq. ( |136| )). This 



condition means that the interference term AiA* will be negligible unless \pd,i — Pd,j\ ^ (^pp.o-, 
that is \5mij\/m < (Tpp^a/p (with m referring to rrii or rrij). In other words, the interference 
term is negligible if the masses are not nearly degenerate. This result is in agreement with the 



conclusion of section 3.2.2: if the masses rrii and rrij are very different, decoherence sets in before 
the dispersion length is reached. 



It is thus possible to work with only one dispersion time T'^'^'^^ = T- 



idisp 



Tf'P. It is 



understood that for noninterference terms, i.e. f dTAjA*, the dispersion time T'^^'^^ will be 



taken to be T'^'^^^ . The integral to be computed is / dTAiA*, with 



AiAj 



exp -i(t)ij{T, L) - fij{T, L] 



(137) 
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with the phase given by 



T 



4>i, (T, L) = 6m,j Vt^-L^ + al,jf ^ ^'^ (T/T^i^P)^ ' ' ^^^^^ 

and the function fij (T, L) defined by 



fijiT,L) = fi{pcl,i) + fjiPdj) - CTpeff 



We would hke to evaluate J dTAiA* by Laplace's method, in the same way as in section 
S.2.2 . This involves expanding the argument of the exponential (137) in powers of small mass 
differences. However, we should take care not to expand nrij around mo, as there is no guarantee, 
in the relativistic case, that 5m j = nij — uiq is much smaller than mo- This is linked to the 
extreme sensitivity on T of the factor \/T'^ — in the relativistic case. Thus the masses rrii 
and rrij should be expanded around rfiQ = (?n.j + mj)/2. A new expansion parameter is defined 
by Sfhj = rrij — rfiQ. The parameters 5mj are said to be of order e. 

Unfortunately, the value T = Tq minimizing fij{T,L) to order e*^ cannot be computed 
exactly. An approximate solution would be any T satisfying approximately pdj — Po and the 
exact solution can be computed by perturbation around it. If T = EqL/pq is chosen as the 
approximate solution, with Eq = \/p'q + friQ (this solution satisfies pdj = Po to order e'^), the 
value of To, to order 6mQ = — rriQ, reads 

Po ^PqEo 

where p is defined by Eq. (115). The expansions in Sfh^ around mo and in 5m j around mo will 
be performed to the same order. It can be checked that the value of T minimizing fij(T, L) to 
order e is still given by T = To (the reason is that mo is the average of mi and mj). 
The value of fij(T, L) at its minimum reads 

^'^^^'^ + TTT- - ''^^ - [-1^) ^^^'^^ 



(139) 



where £ = L / [vqT'^'^'^p) as before, is defined by Eq. (|119| ) and vq = pq/Eq. 
The value of the second derivative of fij{T,L) with respect to T reads 



i^(%i) = ^^+0(^). (140) 

This equation shows that the spatial width (linked to the time width by a factor vq) of the wave 

packet associated to the oscillating particle increases linearly with L for 1 (i.e. at a distance 
much larger than the dispersion length): 

m"^ 

spatial width = axeff V I + ~ i ^ _ ^i^i^ 

Po Eq 
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This expression agrees with the quantum- mechanical result |176| , |185| , derived by observing that 
the relation 6t = ^^L6p, obtained from t = ^ (with L kept fixed), leads to 5L = ^^L. 
The value of the phase at T = Tq is 



kin,L) = ^^L = plL + Oie^), (142) 



where as before dm^j = mf — rnj. The derivative of the phase is given by 



The second derivative of the phase is of order {drhj)"^ and can thus be neglected. Eqs. ( |139| ), 
(pep , (p^ ) and ( p3|) can be compared with Eqs. Q, ([l24|) , O and ( |l2l ). The relation 



6mij = 6mfj/2rho, valid for nearly degenerate masses, may be used. From that comparison, it 
is clear that the final result will be the same as Eq. ( |132| ). 

In conclusion, the flavor-mixing transition probability ( p^ ) for a propagation distance L = 
L 1 satisfying L > L'^'^^p reads 

where the relation 2m^{5mij)'^ + 2{5rn^)'^ = [5mf)'^ + (Sm'j)'^, valid for nearly degenerate masses, 
has been used, and with L°j'^, Ll"^ defined respectively by Eqs. (|128| ) and ( |129| ). Recall that 
we can set Jj^^^^ = = j^disp ^j^jg Yegime^ since decoherence occurs in the transversal- 

dispersion regime if the mass eigenstates are not nearly degenerate. 




As in section 3.2. 2, the time interval AT used to average is assumed to be large enough. At 



first sight, we should have foAT > ^CTxeff in the limit i ^ 1 but this condition is always violated 
at some large distance. Actually, it is sufficient to suppose that AT > (7 = ^ is the usual 



Lorentz factor) so that the phase can freely oscillate around its average value (see Eq. ( |143| ) in 
the limit i ^ 1). This last condition is easily satisfied. Putting together this condition with the 
conditions on AT valid for the no-dispersion and transversal-dispersion regime (Eqs. ( |105| ) and 
(|l3C|)), we obtain (with a^^eff ^ CTxeff)- 

('J^OSC \ 
^,^1. (145) 



The relative strength of the constraints present in Eq. (145) is linked to the relative values of 
the dispersion and coherence lengths: 



rpOSC 

'^xeff ij ^ ^ jdisp ^ jcoh 



Vq t 

It is striking that two different methods of approximation, a priori valid in different regimes, 
give the same oscillation formula, Eqs. (Il32| ) and (144). The dispersion length L'^'^^^ does not 
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play any special role in the final result. Each method is thus accurate enough to be extended to 
the whole range of distances. However, it will be seen in section |^ that the physical interpretation 



disp 



depends on the relative values of L and L- ■ 



6.4 A shortcut 



We now explain another method to derive the oscillation formula (in the transversal- and longitu- 
dinal dispersion regime) which does not require the notion of a dispersion length. Although this 
method is shorter, it is not obvious how to apply it to unstable oscillating particles. Moreover 
it does not make the wave packet picture clearly apparent. 

Returning to Eq. (^), we can choose to integrate on the 3-momentum before doing the 
energy integral. In particular, the integration on the 3-momentum can be done with the help of 
the Grimus-Stockinger theorem JlSO] ]. Let V'(p) be a 3 times continuously differentiable function 
on such that itself and all its first and second derivatives decrease at least like 1/p^ for 
IpI — > oo. Then, for a real number A > 0, 

j'^'^T^^e -¥^(^^l)e^^VO(L-3/2), 

where L = |L| and 1 = L/L. For j4 < 0, the integral decreases like L~^. 

The remaining energy integral in the amplitude (|35|) can be done by a saddle-point approx- 



imation |144|. However, it is quicker to perform first the time average in the probability 



which yields a delta function, and makes one of the energy integrations trivial: 

j dTAiA] = ^ I dEiPiE, qi 1) r{E, qj 1) e^^"'"'?^)^ , (146) 



where ip{E, p) is the overlap function defined by Eqs. (55)-([55|) and qj = J E"^ — m?,. N"^ absorbs 



numerical constants. Actually the time interval of integration AT is finite, so that the the delta 
function is only an approximation certainly valid|^ for AT > T°-^'^. 

Eq. ( |146| ) shows that the transition probability can be interpreted as an incoherent sum 
(i.e. occurring in the probability) over energy eigenstates: interference occurs only between the 
components of -0(^1?) having the same energy [ p49|| . In this way, the correspondence between 
models with and without stationary boundary conditions is obvious: the time-integrated non- 
stationary probability is equivalent to the energy-integrated stationary probability. For example, 
the oscillation formula obtained by Grimus and Stockinger with stationary boundary conditions 



| 15C ] has the form of the integrand in the right-hand side of Eq. (|146| ). This equivalence con- 
firms that the stationary case can be obtained from the more general nonstationary case in the 
limit of a vanishing energy width. Note however that the stationary limit cannot be realized in 
experiments and that the oscillation formula is always averaged over the energy spectrum. 



In section 5.2, we have shown how a wave packet can be associated with the amplitude 



Ai, so that oscillations can be seen, like in the quantum-mechanical treatment, as the result 



As seen from Eq. (|l45|), this condition can be weakened to AT > max{vQ ^a^ejf,^ '^Tl 



^3 
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of an interference between propagating wave packets. Thus the equivalence (|146D shows that 
this physical picture still holds in the case of stationary boundary conditions, provided that 
an incoherent sum over the energy is performed, contrary to what was claimed in Ref. |151]. 
Therefore a wave packet picture can always be associated with the oscillation formula (once the 
incoherent energy average has been done), though this physical picture is well hidden in the 
formalism using stationary boundary conditions. 

This equivalence between the time-integrated nonstationary probability and the energy- 
integrated stationary probability is similar to the equivalence we have met, in quantum- mechanical 
models, between intermediate wave packet models and stationary approaches (see Eq. (|l9|)). 
However, the question of the equality of propagation times does not crop up in the quantum 
field theory formalism. Recall that integrations over microscopic space-time variables are in- 
cluded in the overlap function, with the result that the phase depends only on the average time 
T and L. 

Note that the time average on the probability in sections |6.2| and |6.3| could also have been 
done before the longitudinal-momentum integration in the amplitude, apparently yielding a delta 
function which makes one of the momentum integrals trivial. However this method is spoilt in 3 
dimensions by the prefactor \T — ifi\~'^- Instead of a delta function, one obtains a delta function 
look-alike of width ~ ap/Eo, introducing an additional momentum uncertainty which is 
larger than the mass difference 6m? /E since L°^'^ ^ Po/o'p- For this reason, it was preferable to 
avoid this shortcut in sections |6.2| and (though it yields the same final answer as given by 
the following method). Moreover, it was interesting for the physical interpretation to postpone 
the time average, so as to obtain the explicit dependence of the amplitude on time and distance 
as shown in Eqs. ( |113|) and (|133| ). 

If the coordinate system is chosen so that L is oriented along a coordinate axis, it is easy 
to rewrite the integral ( |146| ) as 

j dTAiA* = j dE e'{n)L-ME)-f,{E) ^ ^^^^^ 

with the definitions fj{E) = fjp{E) + fjjy{E) and 

f'-'-'^^ = 14, + 3< • 

where Vp and po are the components of Vp and po along L, while Eq has been redefined so as 
to absorb the transversal part of po • Vp. The definition of fjo{E) is similar. The geometrical 
constraint 5 (I) is defined by Eq. (|85|). 

Since the phase in Eq. (|147] ) has no stationary point, there is no problem in using Laplace's 
method to integrate over the energy E. The integrand is maximal for 

dm? + 5m?- ^ „^ 

Eij=EQ + p '^^^ ' + 0{e^), 

where the dimensionless number p is defined by Eq. (Ill5| ), and e ~ 5m] ~ 5m'j. 
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The value of fi{E) + fj{E) at the extremum reads 



2 



where axeff,o'm and L°j'^ are defined by Eqs. (^), (|11S| ) and ( |128D , respectively. The value of 
the second derivative of fi{E) + fj{E) at the extremum reads 

The expansion of the phase around the extremum reads 

Stti^ ■ Stti^ ■ 

ME) = --^L + L{E- Ei,) . (150) 

The second derivative of the phase is of order e and can be neglected with respect to the second 



derivative of + fj. The approximation of the integral (147) by Laplace's method is now 
straightforward and yields 



j dTAiA* = voN^apeff^ 



The linear superposition of the different partial transition probabilities J dT AiA'^ gives the 
same result as those obtained in the transverse- and longitudinal-dispersion regimes (Eqs. ( |132| ) 
and ( |144D ). Though the computation is shorter, the origin of the decoherence and localization 
terms is not as clear (see section ^), because 'wave packets' in configuration space do not appear 
explicitly at any stage of the calculation. Moreover the oscillation formula ( p. 11] ) valid in the 
no-dispersion regime cannot be exactly reproduced. 



7 Analysis of the probability in the stable case 

In section ^, the oscillation probability was computed in the three distance ranges studied in 
section |5|, yielding the formulas ( 111 ), (|132| ) and (144). In this section, we analyze the different 



terms appearing in these formulas and we explain their origin. This analysis will allow us to 
answer most of the questions raised in section ^. 

Before analyzing the oscillation formulas, let us sum up the assumptions used in its deriva- 
tion. The oscillating particle is stable, and propagates in vacuum over a macroscopic distance L, 
i.e. L ^ 1/pO) where pq is the average momentum of the particle. The velocity of the oscillating 
particle is arbitrary, but the condition Smjj <^ Eq is assumed to be satisfied. This means that 
nonrelativistic particles are supposed to have nearly degenerate masses (if it were not the case, 
oscillations would vanish anyway). The oscillation formula has been derived for a scalar particle. 
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but this assumption is not very restrictive, since the spin structure factorizes from the sum on 
the mass eigenstates as long as 5m'^j <C Eq. 

In order to obtain a time-independent formula, the transition probability has been averaged 
over a time interval AT satisfying AT > max{yQ^axeffil~'^T°?^). We have also supposed that 



the time separation T**^p between the wave packets, given by Eq. (131), is smaller than the time 
interval AT. The dispersion has been taken into account. 

The oscillation formula valid at large distance (i.e. in the second and third regimes) will be 
studied first since it is the generic experimental case. 

7.1 Oscillations at large distance 

7.1.1 Oscillation formula at large distance 

At large distance (L > po/Cp)) the flavor- mixing transition probability for a stable particle of 
arbitrary velocity propagating over a distance L = LI, with dispersion taken into account, is 
given in a very good approximation by Eqs. (|132| ), ( |144| ) or ( [LSlI) : 



' iff lf>alF?„ \ / I I r A 




where ^(l) is the geometrical factor defined by Eq. (p5[): 

(PO X 1)2 



g{\) = exp 



4cj2 



with Op defined by Eq. (|67|). Recall that the function g{\) restricts the propagation to a cone of 
axis po and angle arcsin{(Tp/pQ) (see Fig. 6). A'^ is a normalization constant determined by the 
conservation of the probability, in the case of a stable oscillating particle: 



/3 



From Eq. ( |152D , we can see that J^fs'Pa^pi^) independent of L, since V^^ = V"^ for stable 
particles (unless the indices a, (5 are restricted to active flavors in the presence of sterile flavors). 
Note that the unitarity relation is only verified to order e because of the approximations made 
in the computations. Thus the normalization constant Ng is fixed by 

Ng j dng^{l) = 1. (153) 

The question of the normalization of the probability will be addressed in more detail in section 



9.1. The oscillation and coherence lengths are defined by Eqs. ( |128| ) and (129), respectively: 



47r^ and T?°'^ = ^ ^^L°f . (154) 
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Figure 6: The probability of detection is maximal within a cone at large distance and within a 
cylinder at short distance. 



The effective width axeff is defined by <ypeff<^xeff = 1/2, with Cpe// given by Eq. (|89|): 



The dimensionless parameter p is defined by Eq. ([1151): 



■z VpivQ-Vp) VDivo-VD) \ , . 

P = ""peff [^ + —2 -2 -2 ) • (156) 



pP pD ^ fi=> ^ f 



The mass width am is defined by Eq. (|119D: 



T - ""leff (-2 ( -2" + ) + ^ ^2 ) • (157) 



As will be made clear in the following sections, the oscillation formula ( 152 ) reduces to the 



standard formulas (|l^) and ( |l4|) (with the additional property of 1/L^ geometrical decrease) if 
the observability conditions |po x 1| < Up, L <^ L1°^ and L°^^ » 5p.£, are satisfied. 

7.1.2 Oscillation phase 



Each term of the sum appearing in Eq. (152) oscillates with a phase given by 2ttL/ L°^'^. This 
phase is independent of the details of production and detection up to order e (except the de- 
pendence on the momentum po, of course). Its origin can be better understood by going back 
to the expression of the probability before the time average is done. 



1. In the transversal-dispersion regime, the phase is given by Eq. ( 120 ). In the limit T <^ 



T'^'^'P, the third term of (Q can be neglected and the phase reads 

0,j(r, L) ^ [Ei - Ej)T - {p, - pj)L . 
Using the expansions ( |114| ) and ( |116[ ), we can write 



5m'i 

EjT - pjL ^ EoT -p^L+—^{L + p{v^T - L)) . (158) 

2po 
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The exphcit expression for p is not useful, as we shah see below. The oscillating phase 
(t)ij (T, L) can now be written 

<P,,iT, L) - 27r— ^ + 2vrp "^^—^ . (159) 



Since \voT — L\ < dxeff (see Eq. ( [113D ), the second term of the right-hand side of Eq. (|159|) 
is negligible if 

\P\ ^ « 1 ■ (160) 

Such a constraint is indeed present in the third exponential term of the probability ( |152| ) , 
so that either the second term of the right-hand side of Eq. ( |159| ) does not contribute to 
the phase, or the corresponding interference term in the probability is negligible. 

Since the leading term of the oscillating phase does not depend on p, the phase is in- 
dependent to order e of the exact values of Eij and pij. In other words the phase is 
independent of the conditions of production and detection. The two fundamental reasons 
for this independence are clearly seen in Eq. (158): 



(a) the particle is on-shell: Ej = y + ; 

(b) the particle is well localized in space-time: {vqT — L\ < o'xeff- 



2. In the longitudinal-dispersion regime, the phase is given by Eq. (138). In the limit T ^ 
rpdisp^ the second term of the right-hand side of Eq. ( |138| ) can be neglected. The phase 
now reads 



4>^J{T,L)^6miJ^/W^ . 

Using the expansions ( |142| ) and ( |143D around the average propagation time Tq, the phase 
can be written in the limit ^ ^ 1 as 

0,,(r,L) - 2vr-4, + (T - To) . (161) 

L^f mo 

The second term of this equation is small if the spatial spread of the wave packet is 
smaller than L°^'^ /^"^ , where 7 is the usual Lorentz factor. Since the possible time range is 
constrained by the width of the overlap function, i.e. |T — To| < l(yxeff/vo (see Eq. ( p.40D ), 
the second term of the right-hand side of Eq. ( [L6ID is negligible if 

i.e. if L ^ Lf-^. Such a constraint is indeed present in the last exponential term of the 
probability (152). If it is satisfied, the oscillating phase is equal to 27rL/L°|'^, as usual. If 
not, the ij interference term vanishes and decoherence occurs. 
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7.1.3 Coherence length 

The last exponential term in the probability ( |152| ) shows that the ij-interference term vanishes 
at a distance larger than the coherence length L'^'-^. This decoherence, predicted by Nussinov 



|21£], has two possible origins, since both the overlap function and the derivative of the phase 



contribute to exp{—{L/Lf-^)'^). The physical explanation depends on whether the coherence 
length is larger than the dispersion length or notf^. 

1. If the coherence length is smaller than the dispersion length (here L'^'^^^ = min{L'f^'^^ , L^*^^)), 
decoherence takes place in the transversal-dispersion regime. In that case, the decoherence 
term comes mainly from the function fij(Tij,L), i.e. from the overlap function. This is 
clear from Eq. ( |123| ), with the dispersion neglected for simplicity: £ <C 1. This decoher- 
ence arises from the progressive separation of the wave packets, due to the different group 
velocities Vi and vj: if T is large enough, both terms \viT — L\ and \vjT — L\ cannot remain 
small. 

2. If the coherence length is larger than the dispersion length (here L'^***' = ^ L^*'^^), 
decoherence takes place in the longitudinal-dispersion regime. In that case, the decoherence 
term does not come anymore from the separation of the wave packets since the packets 
spread out beyond the dispersion length as quickly as they separate (see Eq. (141)). Still, 



decoherence arises from the variation of the phase over the width of the overlap function 



(see Eq. (161)). The interference term is averaged to zero by the time integral when the 
spatial spread of the 'wave packet' becomes larger than L'^^'^/j'^. A similar mechanism has 



been observed in connection with neutron interferometry [186|. 



All in all, the coherence length arises not only from the separation of wave packets, as it is 
usually explained in the literature: in the case of nearly degenerate masses, it can also originate 
in a too large dispersion of the wave packet in comparison with the oscillation length (Fig. 7). 

Remark that the phenomena of separation of wave packets and dispersion contribute to the 
coherence length in both regimes: if approximations such as £ <^ 1 (or ^ ^ 1) are not made, 
terms proportional to come from fij(Tij,L) (or fij{TQ,L)) and from the squared derivative 
of the phase, divided by the second derivative of fij (or of fij). Furthermore, the transition 
between the two regimes is not visible in the oscillation formula, because the two decoherence 
mechanisms smoothly blend. 

Besides the above explanations in configuration space, the coherence length can be also 
understood in energy-momentum space, where it is due to the variation of the phase over the 
'wave packet' width cxpe// (see Eq. ( |150| )). The explanation in momentum space is thus simpler, 
but the two mechanisms at the origin of the coherence length in configuration space are more 
intuitive. 



^^As noted in section ^.2.2| , the condition L""'' > L"^'^*" can be written as ^ < which is equivalent to 

say that the masses are nearly degenerate. 
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III 




Figure 7: In the transversal-dispersion regime (II) , the coherence length is due to the non-overlap 
of wave packets, whereas it is due in the longitudinal-dispersion regime (III) to their dispersion 
(the coherence length was not computed in the no-dispersion regime). 



As dispersion is usually neglected in the literature, it is interesting to derive the oscillation 
formula with this approximation (rj^*^^' = cxd). Such a computation in the transversal-dispersion 



regime yields the same result as when dispersion is taken into account (Eq. (152)), except that 
the following substitution has to be made: 



L 



coh 



T coh 



(FALSE) . 



(163) 



The incorrect multiplying factor has its origin in Eq. (117): in the limit Tj^^^ = oo, the coherence 
length term arises only from the difference between the group velocities vi and Vj. However the 
^) — 1| tends to 1 in the relativistic limit, so that the substitution (|63|) becomes 



factor \p(l — vl 

trivial in that limit. This observation explains why our result (|152|) coincides with Eq. (26) of 



Ref. [144 1, as the authors of this article, while neglecting dispersion, consider only relativistic 
neutrinos. Note however that even relativistic neutrinos spread at large distances so that a 



calculation neglecting dispersion such as in Ref. [144| is only valid for L < L'^*'*^. 



It is interesting to observe that the coherence length increases when a long coherent measure- 
ment in time is performed at the detector, even if the oscillating 'wave packets' have separated 
spatially |182| , |183| ]. In that case, the energy uncertainty at detection goes to zero, aeo — > 0, so 
that the effective width also goes to zero, (Tpejf 0, and the coherence length becomes infinite, 
L^f ^ oo. 

After all these theoretical considerations, it must be said that the decoherence at the level of 
the wave packet is irrelevant in most experiments, since it is usually dominated by decoherence 
effects originating in the energy spread of the beam. Different situations are discussed in section 
9.2.2| . Finally, let us remark that the exponential decrease in L^ of this decoherence term is 
model-dependent. It results here from the Gaussian approximation. However, the definition of 
the coherence length is model-independent, apart from a multiplying constant. 
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7.1.4 Localization 

A third kind of term appearing in the flavor-mixing transition probability (152) are localization 



terms, that is, observability constraints imposing that the oscillation length should be larger 
than the space-time uncertainty: L"?'^ > (T^eff- This condition can be rewritten as Smfj/po < 
Cpeff, stating that oscillations vanish if the energy-momentum measurements allow to distinguish 
between the different mass eigenstates. 

Two localization terms appear in the probability ( |152| ), while one more is implicitly assumed 
when applying the Jacob-Sachs theorem: 

1. The term containing am can be rewritten 

malEl 32alEl 32alEi 



The first term of the right-hand side of Eq. ( |164D can be written as a localization term: 

L!r>^., (165) 
as o"m ~ vqc^p whether the stationary limit is taken or not. The second term of the right- 



hand side of Eq. ( |164| ) is not a localization term and will be discussed in section 7.1.5 . 

2. The term containing p is also a localization term and imposes that 

L'ir^\p\^xeff- (166) 

It is not obvious whether this constraint is stronger than the previous one, i.e. whether 
it is possible that \p\ CTxeff S> o'x- This situation might arise from an energy uncertainty 



much smaller than the momentum uncertainty. In that case, the definition (155) of the 
effective width shows that c^e// ^ o'x- Say, for example, that the energy uncertainty at 
the detection goes to zero, i.e. aeo —>■ 0. With the help of Eqs. (|155|) , (|156|) and (|58|), we 
see that 

\p\ Cfxeff ""-^^ — < Sd , (167) 

O'eD 

where So is the size of the macroscopic detection region. Thus the localization term does 
not give a stronger constraint than L°j^ > Sd- This constraint is always satisfied, as it 
is equivalent to the constraint obtained by averaging the transition probability over the 



production region (see section 9.2.1 ). Therefore, the coherence length can be increased 



without bound by more accurate energy measurements, contrary to what was claimed in 



Refs. [144, 139|. Note that this is not true if the accuracy of the 3-momentum measure- 
ments is increased, as the localization makes the oscillations vanish when the corresponding 
spatial uncertainty becomes larger than the oscillation length. Of course the opposite con- 
clusions would be reached if experiments measured time, not distance. 

3. The contour integral used in the Jacob-Sachs theorem yields a third localization constraint. 
The poles corresponding to the mass eigenstates mi and rrij cannot be both included in 
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no overlap no overlap 



Figure 8: In the no-dispersion and transversal-dispersion regimes (I and II), the localization 
condition is due to the initial spread of the wave packet, whereas it is due in the longitudinal- 
dispersion regime (III) to their non-overlap. 



Table 1: Origins of the coherence length and localization terms. 





Transversal-dispersion regime 


Longitudinal-dispersion regime 


Coherence length 
Localization 


decreasing overlap of wave packets 
initial spread of each wave packet 


increasing dispersion of each wave packet 
constant overlap of wave packets 



the same contour integration unless \5mfj\/pQ < ap. This third constraint is not stronger 
than the previous ones and can be ignored. 

As was the case for the coherence length, there are two origins for the localization term 
containing p (Fig. 8). 

1. In the limit £ <C 1 (transversal-dispersion regime), this term comes from the variation of 



the phase over the width of the wave packet (see Eq. ( 159 )). If the oscillation length is 
smaller than the width of the wave packet, the interference term is averaged to zero by 
the time integration. 

2. In the limit £ ^ 1 (longitudinal-dispersion regime), this term arises from the overlap 
function (see Eq. ( |139| )). Thus it comes from the separation of the wave packets, which 
remains constant in the longitudinal-dispersion regime. 

Recall that the coherence length was shown to arise, in the transversal-dispersion regime, 
from the wave packet separation and, in the longitudinal-dispersion regime, from the variation 
of the phase. Therefore, the coherence-length and localization conditions in the transversal- 
dispersion regime transform respectively into localization and coherence-length conditions in the 
longitudinal- dispersion regime. The situation is summarized in Table 1. 
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7.1.5 Energy-momentum conservation 

Tlie last term of the oscillation probability ( |152| ) to be discussed is the second term in the 
decomposition of Eq. ( [1641) . This term could be a matter of concern since it does not vanish 
in the limit m-j = rrij. Note that this kind of term is not specific to our computation. For 
example, it would appear in the oscillation formula (26) of Ref. |144U if the terms Sa{Ea) present 
in Eq. (22) of that article had been expanded beyond zeroth order in m'^/E'^. 
The second term in Eq. ( 1641 ) imposes that 



|m? + mj — 2mQ\ ^ 



(168) 



77T III' 

As mo is related to the average in- and outgoing momentum by mg = Eq — Pq, condition 
( |16^ ) means that the mass eigenstates have to be on-shell with respect to {Eq,pq) within the 
uncertainty am- For example, this constraint is impossible to satisfy in the mixing of relativistic 
and nonrelativistic neutrinos if the available energy-momentum is such that only the lightest 
neutrino can be produced. However it has no effect on the oscillations in the two cases considered 
in the present report, namely relativistic particles or nearly mass degenerate particles. Condition 
(|16^ ) should simply be considered as expressing the conservation of energy-momentum. Such 
kinematical constraints are usually not included in the oscillation formula, though they rightly 
belong to it. A complete computation of the transition probability should not only include this 
energy-momentum constraint, but also similar terms arising from the interaction amplitudes 
Mp^D, from the prefactor resulting from the Gaussian integration and from the spin structure 
of the propagator. Actually the masses of the oscillating particle should be expected to appear 
not only through mass differences but also through their absolute values. 

7.2 Oscillations at short distance 

At short distance {L ^ Po/Cp), the oscillation formula is given by Eq. (Ill): 



V^^piL) = Ng exp I E 



X exp I -2..^ l^' ^'''^ 



where a^eff and are defined by Eq. ( |106|) and Eq. ( 112 ), respectively. Ng includes all 
prefactors independent of L as well as normalization factors. 



The normalization of the probability can be done in the same way as in Eq. ( 153 ): 



As in the case of oscillations at large distance, it can be argued that the two last exponential 
terms in Eq. (|169|) yield an observability condition not stronger than L°j'^ > Sp^d (see section 
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7.1.4 ) . It could also be interesting to understand the origin of the oscillation phase by going back 
to the amplitude, before doing the time average. This discussion can be found in Ref. ( ||60|] ). 
Actually, the oscillation formula ( |169| ) is very similar to the formula valid at large dis- 



tance (Eq. (152)). There are three main differences. First there is no coherence-length term in 
Eq. ( |169D , which is due to the neglect of the terms beyond 0{€^) in the evaluation of F{Tq) (see 
Eq. ([79|)). Second, the geometrical decrease in is lacking in Eq. ( |16S| ), which is explained 
by the fact that the dispersion of the oscillating particle 'wave packet' is not yet significant. 



Finally, the directional constraint present in Eq. (1169) confines the propagation to a cylinder of 



radius r (with ctx r < Sp^o), whereas the directional constraint present in Eq. (152) restricts 
the oscillating particle propagation to a cone of angle arcsin{ap/po) (Fig. 6 in section [7.1.1| ). 
This different behavior is also a result of the absence of dispersion for L <pQ/a'i. 



In conclusion, the oscillation formula ( 169 ) reduces to the standard formulas (IT^ ) and ([T^ ) 



if the observability conditions |vo x L| < Ivolux and L°^^ » Sp^o are satisfied. 
7.3 Answers at last 

In section ^, our review of field-theoretical treatments brought to the fore several questions: 

• Is there a limit in which the oscillation formula, derived with stationary boundary condi- 
tions, can be obtained from the oscillation formula derived with external wave packets? 

• Does an infinitely precise measurement of the energy lead to an infinite coherence length, 
or is there decoherence anyway? This question is inseparable from the first one. 

• Does a strong localization of the source lead to a dependence of oscillation formula on 
mass prefactors? 

• Do 'plane waves oscillations' exist? 

The external wave packet model used in this report is of sufficient generality to allows us to 
clarify and answer these questions. 



7.3.1 Existence of a stationary limit 

Two important results concerning stationary boundary conditions have been demonstrated. 
First, there is no contradiction between models with stationary boundary conditions and those 
with external wave packets. The former type of model can be obtained from the latter in a 
smooth limit. Let us state again the gist of the problem. Stationary boundary conditions, 
given by Vp^p, = and aep,D = 0, lead to an infinite effective width (Jxeff and thus to an 
infinite coherence length L'l'-^ (see Eqs. ( |154| )- (|155D ). Hence, the localization term depending on 
paxeff/L"j^ seems to diverge in Eq. ( |152| ), with the result that the ij- interference term seems 
to vanish in the stationary limit (Jxeff — >oo. If it were true, it would be impossible to increase 
without limit the coherence length by performing long coherent measurements. This would be in 
contradiction with stationary boundary condition models, which have a zero energy uncertainty 
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but an infinite coherence length. However the product paxeff remains finite, as shown by 



Eq. (167), so that the interference term does not vanish. For example, the Grimus-Stockinger 



model |15C] is obtained in the limit p = 0, CTxeff — > oo with pcTxeff = 0. The latter condition 
means that this model can be recovered from the external wave packet model if the stationary 
limit has the property Vp^^/aepj^ = (see discussion after Eq. ([s^)). The only localization 



condition comes from the momentum spreads of the source and detector (see Eq. (165)). Note 
also that the amplitude Aj (see Eq. (|60[)) takes in the stationary limit the form of a plane wave, 
with a well-defined energy and momentum (see Eqs. (|ll4 )-( |TT6 ) with p 0). 



The second important result is the equivalence ( |146| ) between the time-integrated nonsta- 
tionary probability and the energy-integrated stationary probability. Since oscillations can be 
described in the former model by interfering 'wave packets', this interpretation is also valid in the 



latter, contrary to what was claimed in |151]. Thus the incoherent superposition of stationary 
amplitudes, which have a plane wave form, reproduces the characteristics of wave packet propa- 
gation. Of course, this is not true if the probability is not integrated over the energy spectrum, 
but we have already argued that a source and a detector with zero energy uncertainties do not 
constitute realistic boundary conditions. 

7.3.2 Unbounded coherence length 



As shown by Kiers, Nussinov and Weiss 182|, a more precise measurement of the energy at 



detection increases the coherence length. This feature is reproduced within the external wave 
packet model: L^°^ ~ '^xeff ~ l/ceo (see Eqs. ( |154| )-( p!5^ )). In the limit aeo —>■ 0, the external 
wave packet model becomes equivalent to a stationary model (even if the source is nonsta- 
tionary), as the detector picks up wave packet's components having exactly the same energy. 
Contrary to what was asserted in Refs. |144, |139|] , oscillations do not vanish in that limit since 



we argued above that paxeff remains finite. 



7.3.3 No mass prefactors at large distance 



Contrary to what was claimed by Shtanov [ p42[ | , mass prefactors never appear in front of the 
oscillating exponentials present in the transition probability (152). Let us first explain what 



these 'mass prefactors' mean. Shtanov derives oscillation formulas for neutrinos within a source- 
propagator model in configuration space. Shtanov computes the convolution of the neutrino 
propagator with a source, keeping only the phase in the oscillating exponentials but neglecting 
the contributions depending on the width ax of the source. With this last approximation, the 
space-time localization (i.e. the 'wave packet envelope') of the amplitude is lost. 

Let us consider his model will be considered in a scalar version for an easier comparison 
with our results. The propagator for a scalar particle propagating from x to x' is asymptotically 
{mj\/ {x' — xY ^1) given in configuration space by 

A, -im,^{x'-x)^ /-j^y^N 
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where {x' — x)^ is the Lorentz interval. The convolution of Eq. (|l7lD with a monochromatic 
source yields without problem the standard oscillation formula (^). However, convolution with 
a strongly localized source [(Tx,t ^ ^/E) leads to an amplitude that keeps its dependence on the 
mass prefactor ^yrrij. In that case, the transition probability is not equivalent to the oscillation 
formula ([l^), unless the masses niij are nearly degenerate. 

We now proceed to show that mass prefactors only appear at an intermediate stage of the 
calculation, but not in the final result. Note first that the amplitude ( |171 ) is in correspondence 



with our amplitude (^), computed with the stationary phase method, since the prefactor in 
Eq. ( |133| ) can be rewritten, for T » j^j^^^P^ g^g 



cr. 



peff 



where the definition ( |9^ ) of the dispersion time has been used. This prefactor coincides with the 
prefactor in Eq. (|17lD . However, the subsequent time average of the transition probability com- 
pletely cancels this dependence on the mass, yielding Eq. (|144|) . This can be seen by expanding 
the argument of the exponential in Eq. (133) around the average propagation time. The width 



with respect to T of the amplitude Aj is found to be equal to "^"^^^^ (where 7c; is the Lorentz 
factor associated with the velocity Vd = L/T), thereby providing a mj-dependent factor that 
cancels the rrij prefactor in the Gaussian integration. 

Shtanov does not obtain such a result, since he does not perform any time average on the 
probability. Besides, mass prefactors would remain even if a time average were performed on 
the probability: no 'wave packet envelope' appears indeed in Shtanov's amplitude, so that the 
cancellation mechanism explained in the previous paragraph is not at work. For that reason, 
Shtanov wrongly concludes that the mass prefactor remains if the source is strongly localized. 
Another non-standard oscillation formula, derived by Shtanov for mass eigenstates satisfying 
7712 (x' — x)2 <C 1 <C va\-sj (x' — x)^, reduces to the standard formula when decoherence mech- 
anisms are taken into account. In conclusion, Shtanov's computations in configuration space do 
not lead to new oscillation formulas. 



7.3.4 No plane waves 



loannisian and Pilaftsis |16S] derive a formula for neutrino oscillations which exhibits a plane 
wave behavior if the condition L ^ Po/(^p is satisfied (with ax <C L"?'^, as usual). Their term 
'plane wave oscillations' means that the oscillation amplitude Aj does not decrease as the inverse 
of the distance, that its phase depends on the direction L as pj -L, and that there is no geometrical 
constraint linking the average 3-momentum and the direction of observation. If it were true, 
the oscillation length measured by a specific detector would depend on the direction of the total 
momentum of the initial particles. The plane wave condition is satisfied, for example, by taking 
ax ~ 1cm, L ~ 1000 km and p ~ 1 GeV. These conditions might be verified for atmospheric 
neutrinos. 
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In section 5.1.2 , we have proved that the detection probabihty of an oscillating particle, 
derived under the condition L ^ po/Cp, is negligible in directions other than the average mo- 
mentum whether the stationary limit is taken or not. Therefore the oscillation phase has the 
standard form given in Eq. (|T^ ) and no 'plane wave oscillations' can be observed, contrary to 
the loannisian and Pilaftsis' claim. Also, the absence of the 1/T (or 1/L) factor in Eq. (|77|), as 
noted by the same authors, is easily understood by noting that the absence of dispersion entails 
that the amplitude does not decrease with the distance (the probability is maximal within a 
cylinder as pictured on Fig. 6 of section 7.1.1 ). 

It is interesting to understand why a directional constraint is missing for L < Po/cTp in 
loannisian and Pilaftsis' result, as this fact explains their 'plane wave oscillations' prediction. At 
the end of their computations, loannisian and Pilaftsis obtain an amplitude Aj whose dominant 
term depends on exp{ipj\Ij'\), where |L'| = VT/^ is the complex 'norm' of a complex vector 
L' = L — 2icr^po (their results are translated into our notation through the correspondence 
Qj — > pj, k — > Po, 5P — > 4(7^, L ^ L', / — > L). For L <^ Po/cTp, the quantity |L'| can be expanded 
and the argument of the exponential reads to second order 

(P0-L)2- 



2 I 2 Po ■ L 

Po 



Aalpo 



pI 



The two last terms were neglected in Ref. [163|, though L ^ ax- They lead to the following 
directional constraint: 

(voxL)2- 



exp 



(172) 



where the factor 4 becomes a factor 2 when the amplitude is squared. Therefore loannisian ad 
Pilaftsis' result includes a directional constraint forbidding plane wave oscillations. 

Note that the stationary limit assumption |vp_£,|/(7ep ^ = leads to si = S2 (see Eq. (|109D), 



so that the constraint ( |172| ) becomes equal to the one present in Eq. (111). Thus the condition 
Wp,D\/o'ep,D = seems generic for models with stationary boundary conditions, since it was also 
applied in the case of the Grimus-Stockinger model (see discussion after Eq. ( ]58| ) and at the end 
of section 4.5.2| ) . 

In conclusion, the analysis of the transition probability (111) derived in the external wave 
packet model under condition (|75|) has shown that 'plane wave oscillations' do not exist. This 



negative result was confirmed by a reexamination of the formula derived in Ref. |163|. 



7.4 Conclusion 

Through the discussions of this section, the external wave packet model has shown its power 
and generality. This model can either reproduce in some limit the field-theoretical oscillation 
formulas found in the literature, or allows to understand why they are wrong. Thus the three 
categories of models described in section |3?l| (external wave packet models, stationary boundary 
condition models, source-propagator models) really make one. Furthermore, the external wave 
packet model has the advantage of associating a clear physical picture to the different stages 
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of the calculations. For example, the role played by the dispersion in the observability condi- 



tions (coherence length and localization) could be cleared up in section 7.1, while the threshold 
L = po/ap identified by loannisian and Pilaftsis could be interpreted as indicating the onset of 
dispersion. 



8 Oscillations of unstable particles 

First, we study nonexponential corrections specific to mixed unstable particles, and then examine 
in which respect the oscillation formula obtained in the stable case is modified by the decay term. 
Next, we establish the relation between our formalism and the Wigner-Weisskopf-Lee-Oehme- 
Yang model. Finally, we apply the formulas derived in this section to the B^B^ system. 



8.1 Nonexponential corrections due to mixing 

If the oscillating particle is unstable, the flavor-mixing amplitude (|3) is a superposition of 



mass-eigenstate amplitudes Aj given by Eq. (35). However these formulas were derived on 



the assumption that the diagonalization matrices V appearing in Eq. ( p^ are constant. This 
hypothesis is only true in the limit of a negligible renormalized self-energy. This approximation 
cannot be made for unstable particles, since their decay arises from the imaginary part of 
this self-energy. Actually, nonexponential corrections are always important at large time for 
unstable particles (though until now not observable) , so it is worth checking whether the energy- 
momentum dependence of the diagonalization matrix has a measurable effect on the transition 
probability. 

Nonexponential corrections to the propagation of an unmixed particle of well-defined mass 
are studied in the Appendix, in the context of the derivation of the Jacob-Sachs theorem. 
These corrections have two origins. First, they are due to the bounded character of the energy- 
momentum spectrum of the particle. Second, they are due to multi-particle thresholds, if these 
are included in the energy-momentum spectrum. 

The first type of correction is easily calculated for mixed propagators, since it has noth- 
ing to do with the diagonalization matrices. The corrections can be evaluated separately for 
each partial amplitude Ai, and are negligible (see Appendix). In order to compute the second 
type of corrections, the energy-dependent diagonalization matrices are kept inside the energy- 
momentum integral (^). Let us consider the case of p^p^ mixing in the limit of no CP violation 
(P is a i^T or a -B meson). Starting from Eq. (^), the inverse propagator for the neutral meson 
system can be parametrized by |23^ ] 

where m is the renormalized mass in the degenerate limit and —ifapip^) are the renormalized 
complex self-energies. The diagonal elements are equal because of CPT symmetry. In the limit 

CP 

of CP symmetry, the off-diagonal can be taken to be equal: /oo(p^) = fmiP^)- 



80 



This propagator can be diagonalized as follows: 



-U^2, ( iP^-ml-hip^))~^ 



Thresholds of multi-particle states appear in the renormalized self-energies fj{p'^) and in the 



diagonalization matrix V{p 

Let us define the function If^a{T) by 



.2\ 



Ipa{T) = JdE^P{E,p) Gpa{p^)e- 



iET 



where G^aip^) is the nondiagonal propagator in flavor space and ip{E, p) is the overlap function 
defined by Eq. (p^). The Fourier transform of this function with respect to the 3- momentum 
yields the amplitude A{a /3, T, L). Since we are interested by nonexponential corrections 
to mixing, let us study the effect of a threshold = 6^ on the matrix element Iqq(T). The 
corresponding propagator matrix element is 

'-^OOKP ) — ~2 2 f I 2\ ^ 2 2 f ( 2\ ' 

p'^ — m| — Ji[p ) p — m2 — J2\P ) 

The function Ij3a{T) is computed with the method explained in the Appendix. The integration 
on the energy is done by a contour integration including the two poles zi^2 appearing in the 
right-hand side of Eq. (Il74| ). The function Ip^iT) is equal to the sum of the poles' contribution 
J, the contributions Ji,2 of the energy-momentum thresholds of the overlap function p), 
as well as the contribution of the multi-particle threshold: 

im{T) = J + Ji + J2 + Jb- 

J is given by the residues of the two poles Zj = — irrijVj: 

Since CP violation is neglected, the propagation eigenstates are the CP eigenstates i^i,2 ~ 
lb , so that the matrix V is given at the pole by 

y(.o . Kfe) - ( ; ) ■ 

We shall take into account, in the evaluation of J5, that these expressions are not valid far from 
the poles. An estimate of J is given by 



m 



where m is the mass in the degenerate limit, Am = mi — m2, AF = Fi — F2 and we have 
approximated \/ m? + p^ = m. 
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The contributions Ji and J2, due to the bounded spectrum of the overlap function, can 



be computed separately on each term of the right-hand side of Eq. (174), exactly as in the 
one-particle case. The result is likewise negligible. 

Let us examine in detail the contribution of Eq. ( |216D of the Appendix shows that 
depends on the difference between the complete propagator G^aiz) and its analytic continuation 
Gpa,ii{z): 

J^^ = _i^-^V^'T r duji;{z{u;),p) {Go-o,jj{z)-Goo{z)) e"-^ . (177) 
Jo 

The difference Gn{z) — G{z) can be computed with the parametrization (173). To first order 
in /a/3) we have 

Goojiiz) - Gooiz) = -i {foo,iiiz) - /oo(^)) Go^iz) Gqq jj{z) . (178) 

For large T, the dominant contribution to Jj, comes from the co values near zero, that is z = 6^, 
because of the decreasing exponential. It is thus enough to know how the integrand tends to 
zero near the threshold z = b'^. Just under the real axis {z = x — ie), the spectral representation 
of the self-energy given by Eq. ( |213| ) entails that 

foo,iiiz) - fooiz) = fooix + ie) - /oo(x - ie) = 2iXm/oo(x + ie) . 

The value at the pole of the function fQo{x+ie) can be related to experimentally known quantities 
with the help of Eqs. (|l7|) and ( |T75| ): 

Im fooim'^ + ie) ^ . 

As the main decay channels of and are two-pion decays, the nondiagonal self-energy 
2rTi/oo(a; -|- ie) has a two-particle threshold behavior: 



m AT I X — b"^ 
Imfoo{x + te) ^ — 



Eq. ( |178| ) can then be rewritten as 



Goo, 1 1 (z) -Gooiz) ^ -mAT ^^-^^^—J^GQo{z) Gqqij{z). 

Substituting y = ojT in Eq. ( |177| ) and expanding the integrand in 1/T with the help of the z{uj) 
parametrization given by Eq. (|217|), we obtain 



1 2 o ■ y 

z = r ^" 

It follows that 



V^^+p^ + o 



[m? — 62)5/2 



where we have approximated h"^ + = m? and used the fact that the self-energy is negligible 
at the threshold. Thus Jf, is of the order of 

J,-{QT)-"'^-^l^{b\v). (179) 
Q m 



It is now possible to compare the estimates of J and Jh, given by Eqs. ( 176 ) and (|179| ), respec- 
tively. 
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At small times: Ji, J if 

I Ami T 



Q 



AmT 
sin 



It is known that |Ar| < |Am| (from the experimental data in the case of the K mesons jl56| ] 
and from theoretical predictions in the case of B mesons |115|). Thus, J;, ~ J if T ~ Q~^. 



Since Q is of the order of 0.2 — 1 GeV, for the K and B mesons, nonexponential corrections are 
dominant at small times for T < 10"^*^ s. This time range is not observable. 



At large times: J ~ if 



I ATI I Arl v'^/'^ 



where Tl = min(Ti,T2), that is if 



rLT-3 1n(rLT)~5 1n( ^ ) -21n' 



For kaons, Q~220MeV so that T^T ~ 190. The amplitude is negligible at such large times so 
that nonexponential corrections are not observable. For B mesons, (5~lGeV and AF/F^ can 
be estimated at 4 x 10"^ for the Sj] and at 10~^ for the 5° so that F^T ~ 168 in the 

former case and T^T ^ 162 in the latter. Thus, nonexponential corrections are not observable 
in either case. 

The estimate ( |179| ) of nonexponential corrections to the propagation of mixed neutral 
mesons is in agreement with the theoretical formulas obtained by Chiu and Sudarshan (Eq. (3.30) 
of ||8^), and by Wang and Sanda (Eq. (59) of [ |254| ]). The authors of these two papers com- 
pute nonexponential corrections in a quantum-mechanical framework and extend the Wigner- 



Weisskopf-Lee-Oehme-Yang formalism |191, 172, p09| (see also [|215|] ) 



In conclusion, nonexponential corrections can be neglected in the propagation of quasi- 
stable mixed particles. Therefore, the flavor-mixing amplitude ( p4[ ) with constant matrices V 
can be used not only in the stable case, but also for an unstable oscillating particle. 

8.2 Oscillation formula for unstable particles 

The approximation of the amplitude Aj and of the time average with Laplace's or the stationary 
phase methods can be transposed for an unstable oscillating particle. The decay term 7j(p) 
introduces a new small parameter, namely rrijTj/p'^. This parameter is of order e = 5mf-/2p'^, 
or smaller, since the mean decay length jj^^'^°-y = p/mjFj should be of the order, or larger, 
than the oscillation length L°j'^ = A7rp/6mfj, or else the oscillations would not be observable. 
In parallel with the discussion of section p. 3.1 , the influence of this new parameter on the 



transition probability shows up in the argument of the exponential, through combinations with 
large dimensionless coefficients. 
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Let us inspect the possible combinations of e with the available large parameters, where e 
refers either to {pL"^'^)~^ or to (pL'^^'^"^)~^. First, if the conditions of production and detec- 
tion are not considered, only one such coefficient is available, i.e. pL. This combination gives 
the oscillation term exp{—2TiiL / L"'^^) and the usual exponential decrease of the probability in 
exp{—L/L'^'^'^°'y). Next, the analysis of the propagation with external wave packets adds a new 
parameter, i.e. cTpeff- The combinations (ep/upeff)'^ and (eapeffL)^ give the following contribu- 
tions to the argument of the exponential (with minus signs omitted): 

1. Decay-independent terms, like the decoherence term, (j^)"^, and the localization term, 
C^i^oil )^, which were already obtained in the stable case. 

2. Mixing-independent terms, like and {j0i^)'^- The first is completely neg- 
ligible in comparison with the exponential decrease at order e. The second imposes that 
j^decay ^ a^eff, which is truc for quasi-stable particles. 

3. Mixing-decay cross terms, like ^o^ldLay and ^oscLdecay ■ 

With the help of the condition 0{L'^'""^y) > C'(L°^'=), it can be seen that the new contributions, 
with respect to the stable case, do not impose new constraints for the observability of the 
oscillations. They are thus neglected in the transition probability. 

We have yet to check whether the oscillation phase is modified or not by the widths Tj. 
Let us consider first the transversal-dispersion regime (the no-dispersion regime can be treated 
similarly). The position pj of the maximum is shifted by the decay term ^j{p). The new 
maximum is the solution of 

/,'(p)+7;(p) = o, 



where fj{p) and 7j(p) are defined by Eqs. ( |82D and (l84|) , respectively. It can be computed as 



before to first order in e (e referring to Smj/2EQ or to rrijTj /2Eq) and is equal to 

Pj=pf+6pf\ (180) 



where pj is the solution for Tj = 0, given by Eq. ( |114D , whereas 6pj is the contribution from 
the decay term, given by 

2 

ov^ ^ Po ■ (181) 

The decay time is defined by Tj'^'^"'^ = Eq/ rrijTj. Although the expansion parameter mjVj /2Eq 

(r) 

is multiplied in the exponential by the large parameter T, the shift 6pj is small and the 



expansion ( 180 ) is valid as long as the propagation time T is not much larger than the decay 
time . xhe corresponding energy Ej = \Jp] + j and velocity Vj = Pj/Ej are given by 

Ej = Ef+5Ef\ (182) 
Vj = vf + 5vf\ (183) 
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where E^^^ and v^^^ are the solutions for Tj = 0, given by Eqs. (|116D and ( 117 ), respectively, 



whereas the contributions from the decay term are equal to 6E, 



(r) 



Apart from the modified values of pj, Ej and Vj, the longitudinal momentum integration 
with Laplace's method yields the same results as in section |6.2.1| . The phase of the interference 
term AiA* is thus given, before the time average, by Eq. ( |12C| ), but with the modified values 



|182D and ( |183| ). It can be split in a T— independent part 4>f^{T,L), and a F— dependent 



part (/>j^P (r, L) . The latter is equal to 



4P(r,L) = (^;or-L) {6p\ 



(r) 



5p 



(^)^ 



1+^2 



(184) 



where i = T/T'^'^'P and the definition (H) of T^^^p has been used. Using the condition |foT—L| < 



(Txeff (see Eq. ( |113| )) and the values of pfj given by Eqs. ( |180D and ( p.8lD , we have 



< VQ 



En 



T 



t: 



decay 



ndecay 



which is much smaller than 1 in the observable domain 0{T) < Qi^^decayy Therefore, the 
oscillation phase is not modified by the width: 



Let us now turn our attention to the longitudinal-dispersion regime. Of course, the station- 
ary point Pel J is not shifted by the decay term ^j{p), but the phase receives a contribution from 
the derivative of the decay term. More precisely, the phase (|138|) becomes 



where £ = T/T'^^^p as above. The derivative of the decay term is equal to 

Po L 



(185) 



2E 



Y decay 
d,j J-j 



The terms depending on Tij are proportional to fypeffl'j- They will be small in the observable 
domain 0{T) < ©(T'^^™?^), so that they can be neglected in the phase ( |185[ ). 

Thus the oscillation phase is not modified by the decay term at any distance. The result 
( |152| ), obtained for a stable oscillating particle, is only modified by the usual exponential decrease 
exp{-L/L'fj''°'^), with 

J- decay _ '^PO n 8«^ 

In conclusion, the flavor-mixing transition probability for a long-lived unstable particle (i.e. 
satisfying (Jxeff ^ L'^j'^"'^), having a relativistic or a nonrelativistic velocity (but always with 
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Smfj <^ Eq), and with dispersion taken into account, is given for L > po/ap by 

The definitions of the characteristic lengths and parameters appearing in this formula are gath- 
ered below Eq. ( |152| ). For L < po/^'pj the oscillation formula for stable oscillating particles is 
similarly modified by the insertion of a decay term into Eq. 

8.3 Wigner-Weisskopf effective mass matrix 

Oscillating neutral mesons are usually described in quantum mechanics with the Wigner-Weiss- 
kopf method as extended by Lee, Oehme and Yang [191, |172| , 209 1. In this framework, the 



evolution of a two-meson state satisfies a Schrodinger-like equation with a complex (and non- 
Hermitian) effective mass matrix which can be split into a Hermitian mass matrix M and a 
Hermitian decay matrix T, H = M — iV. It would be interesting to establish a correspondence 
between this method and the field-theoretical formula ( |187| ). 

First, let us simplify the oscillation formula ( 1871) . Note that the coherence length can 



be neglected for particles such as the neutral K or B mesons, since the oscillation length is 
comparable with the decay length: 

r coh T osc 



Y decay ^ t t T decay 

Next, the quasi-degeneracy of the mass eigenstates makes it impossible to set up experimental 
conditions which destroy oscillations. Thus the localization terms (containing the parameters 



Gm and p) drop from Eq. ( 187 ). These approximations lead to the following simplified formula: 



P_ML) = ^^5:y.„Vv;*,y^>*expf--^ . (188) 



The comparison of Eq. (|188| ) with the Wigner-Weisskopf-Lee-Oehme-Yang theory is more easily 
done at the level of the amplitude. Whereas the true amplitude depends on both T and L, the 
probability (|188|) can, equivalently, be generated from the following effective amplitude 



A.jf{a - /3, L) = VMII ^ y^.i exp [-^ {rn, - ^) V,^ , (189) 

where the masses are assumed to be nearly degenerate (rrij = rrij = mo). If M is a diagonal 
matrix the diagonal terms of which are given by rrij — the effective amplitude can be 
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written as 

Aeff{a ^ P,Jj) = — {V 'expl-tM 

= r -'^'^i flavor 



Y) 



13a 



13a 



where Mfiavor = V^'^MV. The mass matrix Mjiavor corresponds to the effective Hamiltonian 
in the Lee- Oehme- Yang theory. We recognize the factor rriQL/pQ as being the classical proper 
time of propagation, common to the two mass eigenstates. However, let us again emphasize that 



this observation does not justify the equal time prescription, which was explained in section 2.2. 
It must not be forgotten that T and L are the average propagation time and length, and that 
interference takes place in a time and space range defined by the width of the effective 'wave 
packet'. 

The effective oscillation amplitude confirms the soundness of the reciprocal basis treatment 



in quantum mechanics explained in section 2^. Recalling that the diagonalization matrix V is 
related to the matrix U diagonalizing the states through V = U^, the Hamiltonian (|To| ) can be 
written as 

Hpropag = \vp> e"*^^* Yjoc <Va\, (190) 

i 

where \j = mj—iTj/2. If the time t is taken to be rrioL/po, the matrix element < \Hpropag\ i^a > 
is in correspondence with the amplitude Aeffia^ (3, L) given by Eq. ( p. 891 ). The field-theoretical 
method will thus give the same results as the reciprocal basis formalism in the case of nearly 
degenerate unstable states. Both formalisms are easier to use that the cumbersome nonorthog- 
onal basis @, since any transition amplitude can be simply computed as a matrix product. 
Since the reciprocal basis formalism has already been applied to the study of CP violation (see 



p38| , |110| , 17, |7^, |243[ |), we shall only give one example in section 3.4 of the application of formula 
(|18S| ) to the computation of a CP asymmetry. 

Let us end this section by an remark on the observability of intermediate unstable particles. 
It is tempting to associate a well-defined mass eigenstate with the exponential behavior in e"*^^* 



in Eqs. (189) or (|190D, but normalized mass eigenstates cannot be defined if the mass matrix 



Mfiavor is not normal (see section 2.1). The reason is that a final state can be produced by all 
mass channels. In practice, it is sometimes possible to maximize the probability of a channel 
(since the lifetimes Tj can be very different). In that case, one talks about decays into a mass 
eigenstate represented by its reciprocal basis vector out<t^j{0) \ [^ ]. 

8.4 Oscillations in the B^B^ system 

B'^B^ oscillations were first measured in 1987 |l^, but the discovery of CP violation in this 
system is very recent jH, We compute here as an example the following CP asymmetry 



A fpo : f r. ^iB^AL)^f)-T{B^,{L)^f) 
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where the index / refers to a CP eigenstate. The propagation eigenstates have nearly equal decay 
widths |115], and must be distinguished by their different masses. Let us define Am^ = mH—rriL, 



where the indices refer to Bl, for Light, and Bh-, for Heavy. 

In the limit of a constant self-energy matrix f^p, the diagonalizing matrix of the mixed 



propagator (|173|) can be parametrized by 

V = ^{'' M and y-i = ^ ( ^ M , (192) 





where o"^ = foo/ foo- Note that the phase of this parameter depends on the phase convention 
chosen for flavor states. The physics should be invariant under 



|5°)^e^^|B°) and \B°) ^ e'^ \B^) . (193) 

Under the ket rephasing ( |193| ), a transforms as a — > e*'^'^~^V. This parameter is related to the 
notations of through a = —qb/pb- The deviation of from 1 parametrizes the amount of 
CP violation in mixing, called 'indirect CP violation'. The constant v in Eq. (|192| ) is an arbitrary 
normalization, underlining the fact that the meson cannot be considered as an asymptotic state. 
Its decay amplitude should always be included in the full amplitude. 

Theoretical calculations (|77] and references therein) show that |(t| — 1 ~ 0(10"^), which is 
much smaller than the present experimental uncertainty, so that only a small CP violation occurs 
in the mixing. On the other hand, a much larger CP violation is expected in the decays ('direct 
CP violation'), since three quark generations are involved in processes such as B^ J/^Ks- 
Thus, CP violation in mixing is usually neglected for B mesons (i.e. \(t\ = 1), which is the 
opposite of what occurs with neutral kaons. 

The computation of the asymmetry ( 191 ) requires the knowledge of the amplitude Tf{L), 



corresponding to the propagation of an initial B^. Starting from the effective oscillation ampli- 
tude (|189|) , we can write Tf{L) as 

TfiL) ~ {MiB^^f) M{B^^f)) V^^ exp|-iM^| V 

where the effective mass matrix M is equal to diag{mL — iT^/l, rriB — iTn/^). It will be useful 
to define the parameter ///by 




Under the rephasing (193) of the kets, fij ^ e '^^ fif. If \fXf \ ^ 1, 'direct CP violation' occurs 



in decay amplitudes to the state /. Whereas the quantities a and fif are not invariant under the 



rephasing (193), the product afij is invariant. This quantity is related to the notations of [[77|] 
through a^if = —Xf. If the final state is a CP eigenstate and the transition dominated by a single 
CKM amplitude (for example B^^ J/iJjKs), crfij can be expressed in terms of the elements of 
the quark mixing matrix, so that its measurement is very important for the determination the 
CKM matrix elements. 
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The amplitude Tf{L) can be written as 
Tf{L) ^—^ — \ {l + aiif)e «) ^ +(i_cr^j)e to " j ^ (195) 

where pq is the modulus of the total momentum of the final pious. 

We also need to compute the amplitude Tf{L), corresponding to the propagation of an 
initial 5°. The same method yields 

Tf{L) — [{l+afif)e to ^ - (^i - (jfij) e to ^po " j . (196) 



With the insertion of Eqs. ( 195[ )-( p^6[) and the approximations \a\ = 1 and Th = ^l, the 
asymmetry ( |191[ ) becomes 

moL 
Po 



where direct CP violation, defined by 



has been separated from CP violation coming from the interference between the mixing and the 
decay: 

.inter/, „o 2/rnA/ _ 2Im (afif) 

A^p ^B,^f)-^^j^- i^i^^p • 

Three different kinds of CP violation have been met in this section: indirect, direct and inter- 
ference CP violations, signaled by \a\ / 1, \iif\ / 1 and Im{a^f) 7^ 0, respectively (see [^] 



p. 78). If the final state is J/ip Ks, a theoretical study (see for example [|115(| ) shows that the 
CP violation parameter is in a very good approximation equal to Xj/^k^ = e~^*^, where f3 is 
one of the angles of the unitary triangle. The CP violation in the neutral kaon system that is 
included in the final state can be neglected. The direct and interference asymmetries become 

^direct ^ j/^ ^ jsinterf ^^0 ^ j/^ ^ _ ^fj . 



The quantities sin 2/3 and |A| were recently measured by the BABAR H, ||] and Belle g, |] 
collaborations: 

s\n2(3BABAR = 0.741 ± 0.067 {stat) ± 0.033 (sysi) , 
sin2/3Be«e = 0.719 ± 0.074 {stat) ± 0.035 (syst) 
IAbabakI = 0.948 ± 0.051 (siat) ± 0.017 (syst) . 

The value of the interference asymmetry is an indication of CP violation in the B system, 
whereas the value of |A| is consistent with no direct CP violation. 
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9 Selected topics 



In section |9.1| , we examine the question of the normahzation of the transition probabiHty. We 



then discuss in section 9.2 the relation between coherent and incoherent effects, as well as the 
different estimates of the wave packet size. Finally, the case of an unstable source is reviewed 
in section Wl 



9.1 Normalization of the transition probability 

Until now, we have not explained why the time-averaged squared modulus of the amplitude 
(see Eq. (|9^)) can be interpreted as the flavor- mixing transition probability. Actually, the 



interpretation of the expression (152) as a transition probability was a bit of a guess. As a matter 
of fact, the steps between the amplitude A{a — > /3) = V^j^AjVja and a normalized event rate, 
containing an expression which can be interpreted as a flavor-mixing transition probability, are 
not straightforward. In the case of the large-distance oscillation formula, the constant prefactors 
voN'^apeff appearing in Eq. ( |132| ) (or in Eqs. (144) and ( 15lD ) have been shoved into a constant 



Ng, which was normalized in Eq. (|153D in order to obtain a unitary evolution in the stable case. 



The normalization of the short-distance oscillation formula ( |111D proceeded in the same way. 
Whereas these prefactors can be factorized from the sum on the mass eigenstates if 6mfj <C Eq, 
this is not possible for a mixing of a relativistic rn-j and a nonrelativistic nij mass eigenstate. 

Cardall has made an attempt to go all the way to a normalized event rate. Although 
he claims to resort to arbitrary external wave packets, his model is equivalent to the Gaussian 
external wave packet model, since he uses Gaussian approximations to perform the integrals. 
We shall explain how Cardall's argument can be applied to the external wave packet model 
developed in sections Q to |^. Only the large-distance case will be considered (L > Po/o"p); the 
short-distance case can be treated similarly. 

The expression .A(a— >/3, T, L).A*(a^/3, T, L) (see Eq. ([3^)) is proportional to a transition 
probability between one-particle states. In order to obtain a formula applicable to experiments, 
the external wave packets should be interpreted as densities of particles. Cardall proposes three 
rules of correspondence: 

1. Let us first examine the external wave packets normalization constants Np.^, Np^^^, Np).^ 
and Np^^^. They come from the overlap function and are included in the factor A^^ 
appearing in the transition probability ( [132| ). Eq. ( |52D shows that Np^^ is proportional to 
the modulus of the initial wave function at production. The constant Np.^ is interpreted 
in the same way. The constants Np_^^ and A^|,.^ are then interpreted as the initial state 
particle densities at production and detection, respectively: 

Nl, /(Q,xp) and Ar2 ^ fiQ',Xn), 

where / is the phase space density. The normalization constants Np^^^ and N^^^^ are 
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interpreted as particle densities for the final state particles at production and detection: 
2 dK , ..2 dK' 



Ni ^ - — and N 



(27r)32i?p„„, (27r)32i?z,„„/ 
Note that there may be more than one particle in the initial and final states. 

2. Let us now consider the production and detection widths. Recall that the factor A^^ in 
Eq. ( [13 21 ) contains a volume factor V"^ = (2~^7r^(Tpp CTepCTpn deo)'^ coming from the overlap 
function (see Eq. (|54|)). One factor V is interpreted as an integration on macroscopic space- 
time variables: 

V dxp dxp dx^, dx^jj . 

Since the production time is not measured, it is integrated over in the event rate. This 
integration is equivalent to the time average over the propagation time T done in the 
present report. 



3. Finally, the effective width cTpeff, appearing in front of the exponential in Eq. ( |132| ), is 
interpreted as the result of an integration over the energy spectrum of the oscillating 
particle. Furthermore, the second volume factor V, the geometrical factor 5^(1) and the 
term including am are interpreted, in the limit mj = rrij = mo, as a product of delta 
functions constraining the direction of L and the energy of the oscillating particle. More 
precisely, 



where p = {E, \J E"^ — tjiq 1), and with the approximations Pp = Pd = Po, Vp = Vo = 0, 
so that the energy and momentum decouple. The relationship is exact in the limit app = 
o'pD = (with the approximations for the momenta and velocities just mentioned). 

These three rules lead to the macroscopic event rate at the detector at time x^: 

dr{xl) = J d^pj d^^J ^/(Q,xp) J ^/(Q',x^)dr(Q,Q',Xp,x^), 

with the constraint x% = x^ — L/vq. The single particle event rate is given by 

dT{Q, Q', Xp, Xb) = j dE flux x Pmix x da . 

The flux is the flux of oscillating particles of energy E, produced at Xp and detected at Xp,- It 
includes a phase space factor for final state particles at production, J dK, as well as the delta 
function 6^^\p — pp), the interaction vertex \Mp{Q, K)\'^ , the geometrical factor l/L?, and a 
velocity factor \vq — Vp,\- 

The factor da is the cross section for the interaction of particle of mass rno in the detector. 
It includes a phase space factor for final state particles at production, dK.' (not integrated over 
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if these momenta are measured), as well as the delta function 5^^\p — pi:,), the interaction vertex 
|Md(Q',K')P and the M0ller factor \vo - Vd\''^. 

Finally, Pmix is identifiable as the flavor mixing transition probability: 

In comparison with Va^fsi^'), given by Eq. ( |152D , the geometrical decrease and the normalization 
factor, i.e. Ngg^{\)/L'^, are included in the flux and in da. The approximation tuq = {rrii + mj) /2 
has also been made in each interference term. 

In the case of a mixing of relativistic and nonrelativistic particles, the flavor-mixing proba- 
bility does not factorize from the amplitudes of production and detection. The interference terms 
are however negligible because decoherence occurs. The dependence on the mass rrij should be 
kept in the noninterference terms J" (iT|^jp, with the result that the prefactor vq in ( |132| ) should 
be replaced by Vj. Velocity-dependent prefactors were already derived in the intermediate wave 



packet model of Giunti, Kim and Lee |141], except that the one-dimensional treatment in that 



article leads to a peculiar dependence in v'J^ (in our case, the conversion of the l/T^ prefactor 
into a v'j/L?' prefactor yields an additional f ^ , leading to the expected linear dependence in Vj ) . 
Cardall's correspondence rules are still valid, although the flux and the cross-section now depend 
on the mass eigenstate. For example, mo and vq are replaced respectively by nrij and Vj in the 
third rule. The single particle event rate is replaced by 

dr(Q,Q',Xp,Xo) = X]l^J°l^l^/3i^l^ / dEfluxj X daj. 
j 

Mixings of relativistic and nonrelativistic neutrinos have been studied by Ahluwalia and Gold- 
man Q, who identify the third mass eigenstate with the 33.9 MeV particle suggested by the 
KARMEN experiment [ p6| . Note that the amplitude of production (and also of detection) 
should be computed separately for the relativistic and nonrelativistic neutrinos. 

9.2 Coherence, incoherence and decoherence 

Generally speaking, decoherence is said to occur in particle oscillations if the interference terms 
in the transition probability are averaged to zero by some mechanism. In that case, the tran- 
sition probability becomes independent of the distance. We have seen in section ^ that this 
phenomenon appears in the large-distance flavor- mixing transition probability ( 1521 ) through the 



localization term and through the coherence length. Localization conditions, such as L°*'^ > ax, 
determine whether there is decoherence from the start. If decoherence only occurs beyond a 
certain distance, the threshold is called the coherence length. These two phenomena are closely 
connected, since the coherence-length and localization conditions in the transversal-dispersion 
regime transform into the localization and coherence-length conditions in the longitudinal- 
dispersion regime, respectively. In the end, both mechanisms of decoherence originate in the 
wave packets widths of the external particles. 
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Other effects lead to similar constraints on the oscillations. They can be classified in coherent 
or incoherent eff'ects. In quantum field theory, a coherent efi'ect has to be taken into account in 
the amplitude, whereas an incoherent effect is incorporated into the computation only at the level 



of the probability |151]. Actually, this distinction has to be made only because approximations 
made in the computations often make the intrinsic decoherence effects disappear. For example, 
there are endless discussions about the energy and momentum coherence in the plane wave 
treatment of oscillations, since this approximation destroys all natural decoherence mechanisms. 
This should be contrasted with the fact that coherence-length and localization conditions appear 
explicitly in the oscillation formula ( |152| ) obtained in the external wave packet model. Moreover, 
decoherence between different energy components is automatic in the external wave packet model 
(see section ^.4| ). Note that one should be careful not to confuse the term '(in) coherent', referring 
here to a constraint applied to the amplitude or to the probability, with the term '(de)coherence', 
referring to the existence or disappearance of oscillations. 

9.2.1 Incoherent effects 



First, the energy-momentum spread of the beam has to be taken into account |37, 36, 34]. For 
example, the average of the oscillation term over a Gaussian momentum distribution of width 
Ap and mean value p gives 



/ 



/ L (p-p)M / L 2 

dp exp -2^2—— - ^TT-r-r^ ~ exp -27rz_ - 2tt 



L°f 2(Ap)2 J ^ I \^ p L«f 



where L°?'^ = 4:'Kp/5mf-. A new coherence length can be defined by 



T coh ^ P TOSC 

~ V2^ Ap ■ 



This new coherence length can be obtained from the coherence length ( |154| ), derived in the 
external wave packet model, by substituting the beam spread Ap for the effective width Cpe//. 
This result could be expected from the discussion of section by extending the incoherent 
sum over the energy from a Cpe// range to a AE = Ap/vQ range. 



Second, the macroscopic propagation distance L is not perfectly known |149, 66|. A Gaus- 
sian average over the macroscopic region of production gives 

L (L-Lf\ I L olAL''^^ 

dL exp -2m , ^ ' ~ exp -Itti 2t:^ 

^' 2(AL)2^ 1^L°- 

where L is the average propagation distance and AL is the size of the source. Thus a new 
localization condition has to be satisfied: 

This condition is similar to the localization condition (|166| ) which comes from a coherent effect. 
The only thing to do is to substitute the size of the source AL for the effective width (Jxeff- 
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As noted several times ( [|104 |94 f37| , |182| , |24q] , and references therein), it is impossible, in 
stationary cases, to distinguish wave packets from an incoherent plane wave superposition with 
the same energy-momentum spectrum. 

This observation leads directly to the generalization of the quantum field computation of 
the transition probability already done with external Gaussian wave packets. Arbitrary wave 
packets, with space width and momentum widths not minimizing the uncertainty, lead to an 
oscillation probability similar to Eq. ( |152| ) (or Eq. ( |187D if the particle is unstable) , except that 
the relation CxeffCpeff = 1/2 does not hold anymore. This result can be understood as follows. 
If the arbitrary external wave packets are decomposed in Gaussian wave packets, the amplitude 
Ai becomes a superposition of Gaussian amplitudes. These do not interfere with each other if 
their phases are very different. Thus each Gaussian amplitude will mostly interfere with itself, 
from which an incoherent superposition follows. 

In practice, the size of the region of production is usually much larger than the size of the 
wave packets. Similarly, the energy-momentum spread of the beam is usually much larger than 
the energy- momentum spread of the wave packet. Incoherent effects are thus very often domi- 
nant. At worst, the sizes of the regions of production and detection and of the energy-momentum 
spectrum are determined by the characteristics of the wave packet. Similar decoherence mech- 
anisms were found by Gabor |131| in connection with electron-interference experiments. This 
author obtains three decoherence factors, coming from an average over the size of the source, 
from an average over the energy spectrum, and from the path difference between the interfering 
beams. The two first effects have been explained in this section, and the last effect can be related 
to the wave packet separation in particle oscillations. 

9.2.2 Coherent effects 

Recall that coherent effects are constraints on the oscillation process, which should be applied 
at the level of the amplitude. Most coherent effects can be expressed through their influence 
on the sizes of the external wave packets associated to the production and detection of the 
oscillating particle. Unfortunately, the estimate of a wave packet size is not an easy matter and 
no consensus exists on whichever evaluation method is the best. 

Let us consider first a solar neutrinos. Only the cases of line spectra {pep or '^Be) de- 
serve careful thought, since the energy average always dominates coherent effects for continuum 
spectra (except if a detector with a extremely high energy resolution is invented). The most 
commonly discussed constraint on the wave packet size comes from the pressure broadening, 
that is, the interruption of coherent emission due to collisions of the emitting atoms. Nussinov 



|21£], Loeb [20C]and Kim and Pevsner [|185[| estimate the wave packet size of the parent nuclei 
at Ux ~ 10~® cm, whereas Krauss and Wilczek J188| | propose Ux ~ 10~^ cm. Kiers, Nussinov and 
Weiss [182] claim that the small wave packets of captured electrons give a stronger constraint 



(Tj; ~ 6 X 10 ^ cm. However, none of these estimates gives a momentum width (T„ as large as the 



IkeV energy spread of the solar neutrino line spectra [|l8q, 221, pq]. The main contribution to 



this energy spread comes from the thermal energy spread of the captured electron |188], as well 
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as from the Doppler shift due to the thermal motion of the emitting nucleus [221|. These effects 
cannot be modelized at the level of the amplitude and belong thus to incoherent effects. Finally, 



Malyshkin and Kulsrud |204] compute the effect of Coulomb collisions on the solar neutrino flux. 
They obtain a quantity cJa, equivalent to the wave packet width, and consistent with Nussinov 
and Loeb's estimates. While Malyshkin and Kulsrud agree that pressure broadening can be 
neglected with respect to the line width, they claim that it could be relevant to the decoherence 
in the case of the continuous solar spectra, for a detector of very high resolution and a very long 
oscillation length. This conclusion is in contradiction with the increase of the coherence length 
with the resolution of the detector: the coherence length is mainly determined by the energy 
resolution of the detector, if it is smaller than the wave packet size at the source. All these 
results show that coherent effects seem to be irrelevant to solar neutrinos. 

As regards neutrinos from supernovae, the source wave packet width has been estimated 
by Anada and Nishimura at ax ~ 10~^^ cm for neutrinos from the supernova core [^], and 
at ax ~ 10~^ cm for neutrinos from the neutrino sphere pO[| . Since p~10MeV, core neutrinos 
decohere before oscillating (p/(7p~0.1), whereas sphere neutrinos might oscillate if the oscillation 
length is not too short (p/fTp~ 10^). However, the incoherent momentum spread {p/ Apr-^l — 10) 
destroys any interference effect left, unless the oscillation length is, by chance, comparable to 
the supernova-Earth distance [ p36| ]. 

Another type of neutrino source is a radioactive nucleus in an atomic lattice, which is the 
case for reactor neutrinos. The wave packet size can be estimated by Rich [237| and Grimus and 
Stockinger jl50| ] at a^^lO"^'^ m. The wave packet size of the emitted electrons is larger and can 
be neglected [185]. With the neutrino energy around 1 MeV, decoherence occurs beyond ^ ~ 10^ 
oscillation lengths, so that there are no coherent effects relevant for laboratory experiments. Of 
course, oscillations may vanish because of the incoherent energy average. 

Coherent effects may also appear because of the finite lifetime of the source Tdecay, which 
interrupts the classical emission of the wave train and limits the size of the wave packet to 
~ CTdecay [182|. This effect could be relevant to atmospheric neutrinos and in accelerators. 
For quasi-stable sources, the constraint L°j*'^ > ax can be very stringent. For example, in the 
case of vr ^ /xi^, the length of the wave packet is bounded by CT^ecay ~ 7.8 m and could be 
macroscopically large. However this argument is only valid as long as the decay point of the 
source is not observed at all. As emphasized by Kayser |176|] , the detection of a final state, such 
as the muon in vr — > /xz/, can localize the decay point to a precision much better than either 
CTdecay or the macroscopic size of the production region. An extreme example is given by the 
neutrinos from the (3 decay of a nucleus with Tdecay ~ 1 sec. The observation of the /? particle 
allows to pin down the decay point to a precision much better than either CTdecay ~ 10^ km or the 
production region, for example a nuclear reactor. Thus the relation between the decay time of 
the source and the wave packet size of the oscillating particle is not direct. The decay time only 
puts an upper bound to the wave packet length. Some papers dealing with the finite lifetime of 
the source are reviewed in section |9.3| . 

Since the coherence length depends not only on the characteristics of the source, but also 
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on those of the detector, the minimal wave packet size in the detection process must also be 
checked. It can be roughly estimated at (Tx~10~^'^ m, that is, (Tp~ 10-3 MeV ^5^]. In all cases, 
this momentum width is not larger than the energy spread and can be neglected. 

As regards K and B mesons, particle decay takes place before decoherence occurs: 



coh 



PO 



■nin 



PO 



> 1, 



where Lf-^ and ]^'^'r'^°-y are given by Eqs. (154) and (186), respectively. 



To sum up, coherent effects are most likely irrelevant with respect to incoherent effects. 
This discussion also shows that, in practice, we are far from being able to increase the coherence 
length by more accurate energy measurements. 



9.3 Unstable source 

In principle, the instability of the source can be taken into account in a field-theoretical model 
by considering the source as another internal line of the global Feynman diagram describing the 
process. The difficulty, of technical nature, consists in integrating on both propagators, with a 
constraint on the decay point which can be either nonexistent, or very stringent. 

The case of an unstable source decaying in flight has been considered by Campagne in a 
field-theoretical model |84]. Besides the usual condition L"^'^ > (Txp.d^ this author obtains a new 



localization condition > Lp^™^, where L'jf'^^'y ig either the decay length p/mF of the unstable 
source, or the length of the decay tunnel, if the latter is shorter. However, this treatment is not 
a real improvement on the one proposed by Rich, who had derived earlier the same results in the 



framework of time perturbation theory in quantum mechanics | 237 |. It is not satisfactory that 
macroscopic quantities, such as the lengths of the source, of the target and of the pion decay 
tunnel, are treated in |84] on the same footing as microscopic quantities. 

The case of an unstable source at rest has been considered in detail by Grimus, Mohanty 



and Stockinger 151 , 152]. They use their previous field-theoretical model | 150 |, modified by a 
quantum-mechanical Wigner-Weisskopf approximation, in order to take into account the finite 
lifetime of the source. Besides the already known condition L^j"^ > axp,D-, they obtain a new 
localization condition: 

^ (197) 



< — L'? 



nipTp 

where (Tp, rup and Fp are respectively the wave packet momentum width, the mass, and the decay 
width of the unstable parent particle. In brief, the unstable source should not move on distances 
larger than the oscillation length during its lifetime, or else the oscillations vanish. They find 



that this condition is most likely satisfied in the experiments LSND |28, and KARMEN 
1 26, 27]. The same authors have also found a new coherence length, Lf-^ = AE"^ /dmfjFp, which 
is however completely irrelevant, since it is much larger than the coherence lengths discussed 
previously (Fp <^ Op for all weakly unstable particles). The widths axp,D can be estimated at 
10^2 MeV for LSND and KARMEN j|52|, so that the condition Ll^"" > Gxp^d is also satisfied. 
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A slightly different model for vr ^ /iz/ is proposed by Dolgov |105]. This author computes 
the transition amplitude by coupling the neutrino propagator to external wave packets at the 
source, but not at the detector (source-propagator model). The finite lifetime of the source is 
taken into account with a Wigner-Weisskopf approximation, like in the articles discussed above 



1 151, 152], but without the restriction of a source at rest. Dolgov considers first the case where 
the decay point of the muon is not registered (but its energy-momentum is perfectly known) 
and computes an amplitude of spatial width 7 = p^^ / m-,^T . The oscillating phase is given by 

(t>ij{T, L) = 27r— : + a -—^^ with a = 



TOSC TOSC -_2 ' 

^ij ^ij Vq - • Vo 

where v^^ and vq are the velocities of the pion and neutrino, respectively. With the help of the 
constraint < vqT — L the second term of the phase is seen to be negligible if L°^^ » 07. 
For a source at rest, q = and this localization condition vanishes. When the muon decay 
point is registered, Dolgov obtains an oscillation phase equal to the standard result as long 
as the localization condition a-,^ <^ L°^'^ is satisfied (where is the size of the wave packet 
of the initial pion). It would be interesting to study intermediate situations where the muon 
is registered with a space-time uncertainty and to compute the time-independent oscillation 
probability. 

This brief account of the treatment of oscillations from an unstable source shows that there 
is still work to do, especially regarding decays in flight. 

10 Correlated oscillations 
10.1 Introduction 

Experiments where two correlated mesons oscillate together are very interesting for the study 
of CP violation |82, ^] and can provide tests of the Einstein- Podolsky- Rosen (EPR) effect 



|2|, 1^ m, |TT|, |228|. At the present time, the process (/)(1020) K^K^ is studied by the 
experiment DA$NE at Frascati [^. The process T(4s) — > B^B^ is studied at B factories such 



as KEKB at Tsukuba [Q] or BABAR at Stanford |3^, 33 1, where asymmetric collisions allow 
the measurement of the B pathlengths. In principle, correlated oscillations could occur with 
particles not conjugated to each other, for example the lepton could oscillate between e and ^ 
and the neutrino between Ue and in the process tt^Iv. However, the charged lepton masses 
are too different for such oscillations to be observable (but note that if the masses were close 
enough to allow oscillations, it would be impossible to identify the flavor). 

In this section, we develop the formalism appropriate for correlated oscillations and we apply 
it to oscillations of correlated mesons. We then disprove a claim of a non-standard oscillation 
length. The question of the oscillation of recoil particles, like A in np^AK or fi in vr^/uz/, can 
be examined in the same framework, allowing us to answer the flfth question posed in section 
33. 
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The treatment of one-particle oscillations with the external wave packet model of section § 
can be easily extended by associating a propagator with each oscillating particle. Consider the 
two successive processes: 

(e+e-)(pp) ^ nX, (198) 

n ^ PlP2^ fl{Pn,)f2{pD,), (199) 

where pp and ^ are the average momenta of the source TZ and of the final states, respectively. 
We would like to apply the external wave packet model to the second process. In principle this 
model should be modified to include the decay width of the source TZ, since one external particle 
at an interaction point is not sufficient to localize it. However the decay widths of the (/)(1020) 
and T(4s) are large, so that their decay point will be close to their production point. The latter 
can be localized with the external wave packets. In the case of a small decay width, it will 
be seen below that the correlated oscillations become independent of the production point. 

Note that the initial flavors of the oscillating particles cannot be observed at the source and 
must be summed over. For example, B^B^ and B^B^ are both produced in the T(4s) decay. In 
contradistinction to the one-particle oscillation case, this sum does not destroy the oscillations, 
because of the flavor correlation between the oscillating particles. Indeed, two mesons of the 
same flavor cannot be observed in the above process at the same time in the center of mass 
frame (EPR effect). 

Let us first write the amplitude corresponding to the process ( |199D . The masses of the first 
particle are noted m-ia, etc., while the masses of the second particle are noted m2i, m2j etc. 
The two sets coincide if Pi and P2 are conjugated mesons, but it is not necessarily the case. The 
average production and detection points are noted Xp and y^,^, yn^-, respectively. The partial 
amplitude, corresponding to the propagation of the eigenstates of mass mia, from Xp to y^i) 
and of mass m2i, from Xp to yo^' '^^'^ written as 

with the overlap function given by 

i^iPl,P2) = A^V'p(rf +P2)V'Di(P?,Pi)V'd2(P2>P2) • 

The production overlap function is defined by Eq. (p4|). The detection overlap function ipp,^ 
is defined by the same equation, where po^, Ep,_^ and Vp,^ have been substituted to pp, Ep and 
Vp, respectively. The definition of iI^d^ is similar. 

The integration on the energies Pi 2 can be done with the help of the Jacob-Sachs theorem 



( |211| ), yielding 

j d'pi j d^P2Hzi,Z2,Pl,P2)e~'^/^^^'+'P'-'^' e-'V^T, + ^p,■U ^ (200) 

where Ti^2 = 2 ~ ■'-'I.2 = y^i 2 ~ average propagation time and distance of Pi, 

and Zi is the pole of the propagator Gi{p'^). 
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Since our aim is not to prove the conservation of energy-momentum between the initial 
state i;^)(1020) and the final states /i,2, we set 



Poi + = Pp and + E^^ = Ep . 

As before, it will be useful to define reference masses mi and m2 through 2 = 2 ~ 2 • 
The velocities vi^2 are defined by vi^2 = Pdi 2/ ^d-^ 2 (they should not be confused with v^,^^! 
which appear in the overlap function and refer to the velocities of the detection regions, as 
explained in section 4.5.2| ). 



10.2 Factorization 

The momentum integrations in the amplitude (|20C| ) can be easily evaluated, either with Laplace's 
method or with the stationary phase method, provided that the energy-momentum correlation 
is not too stringent at the source. More precisely, the integrations on pi and p2 can be done 
independently if the energy uncertainty at the source is larger than the energy uncertainties at 
the detection points (recall that the momentum width is always larger than the energy width). 
For example, the production widths for the (/>(1020) decay and the T(4s) decay can be estimated 
by their average decay widths, which are 4.26 MeV and 14MeV, respectively |156| ]. Both widths 
are larger than typical detector uncertainties. In that case, the energy-momentum width at the 
source can be neglected and the overlap function factorizes in one-particle overlap functions: 

'>P{zi,Z2,Pl,P2) = -01(^1, Pi) V'2(^;2,P2) • (201) 

Thus the resulting amplitude factorizes in one-particle oscillation amplitudes. The rest of the 
computation proceeds as in the one-particle case, except when the detection times are measured: 
the time average is then done only once, on the production time. 

If the energy-momentum uncertainty is smaller at the production point than at the detec- 
tion points, energy-momentum correlations are introduced at the source so that the amplitude 
cannot be factorized. However, the oscillation formula is not expected to be modified, as long 
as the energy-momentum uncertainty at the source is larger than the mass difference between 
the interfering mass eigenstates: app > 5mfj/pD^ 2- This condition resembles the localization 
condition derived in the one-particle oscillation case, which stated that oscillations vanish if the 
uncertainty on the position of the source is larger than the oscillation length. This similarity is 
misleading: it will be seen correlated oscillations do not vanish if Upp < 6mfj/pDi 2, since the 
knowledge of the source energy-momentum is not sufficient to ascertain which mass eigenstates 
are produced. This is because the energy-momentum of the source is shared between two inter- 
mediate particles, instead of one as in the one-particle oscillation model presented in section ^. 
Moreover, the position and time of the production process can be reconstructed, in principle, 
from the final states characteristics, so that the oscillation pattern is not washed out by an ill- 
defined production point. The discussion of the above constraint is somewhat academic, as it is 
satisfied for the experimentally studied processes, (p — > K^K^ and T(4s) — > B^B'^. Nevertheless, 
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the question of whether the violation of this constraint changes the oscihation formula will be 
examined, partly as a matter of principle, and partly because it will be useful for the discussion 
of the recoil oscillation conundrum. 



10.3 Energy-momentum correlation at the source 

There are two reasons to study more carefully energy-momentum correlations at the source. 
First, this analysis yields a correlated localization condition, arising from the source. In partic- 
ular, it allows to check explicitly that a zero momentum width at the source does not wash out 
the oscillations. Second, it will be useful for the examination of the so-called recoil oscillations. 
Unfortunately the treatment of the full 3-dimensional case is involved, because the correlation at 
the source might link a variation in the longitudinal momentum of one particle, with a variation 
in the transversal momentum of the other. If the transversal dispersion is neglected, the inte- 
grals can in principle be evaluated with Laplace's method in three dimensions, but the results 
are lengthy. For that reason, only the collinear case will be presented. Momenta, velocities and 
lengths can take positive or negative values, with the sign indicating the direction. 

Since the stationary phase and Laplace's method give the same results (see section |6|), we 
choose the latter for the longitudinal momenta integrations. Besides, the dispersion will be 
neglected. Recall that the masses of the first particle are noted niia, mib etc., while the masses 
of the second particle are noted m2i, m2j etc. Suppose that the overlap function ip is maximal 
for pi = pia and p2 = P2i- As before, these momenta can be computed by expanding the 
argument of the overlap function in small mass differences. The expansion parameters are given 
by Suii^ = mfjj — and (Jmgj = m^j — m2- At first order in mass differences, the average 
momenta of the particles are 

2 2 

Via = Pd, + , (202) 

V2i = Pd, + ^ , (203) 

where A = (^xeffi'^xeff2 ~ ^xeff3 — ^- momentum space, the eff'ective widths read 



1 , {Vl,2-Vpy , 1 , ivi,2-VD,^2f 



2 2~^2~^2~^2 ' 

^ = (204, 

"peff3 "pp "eP 

They are related to the effective widths in configuration space by o'pe//i,2Cxe//i,2 = 1/2- Their 
name is justified below by their appearance as second order coefficients in the expansion of 
the overlap function. These expressions can be compared to the effective width ( |89| ) in the 
one-particle oscillation case. The mass differences are included in the constants c^: 

Ck = 2 + 2 ('^1 + "^2) , 
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with 6i = 5m\^/2Eoi and 82 = bm\j2Ej^^. 

If the overlap function is noted = -^exp(— /(pi,p2)), the expansion of the argument of 
the exponential around its extremum can be written as 

fl \ fl \ , {Vl-Plaf , {P2-P2if . {Pl - Pla){P2 - P2i) 
f{Pl,P2) = f{Pla,P2i) + + + • 

The evaluation in one dimension of the integral ( |200| ) by Laplace's method yields 

Aa,i ~ exp {-f{pia,P2i)) exp -lEiaTi + ipiaLi — tE2iTi + ip2iLi — 

>< exp ^(^al^ff2i'"iaTi-Lif + al^j:f^{v2iT2-L2)'^ - 2al^ffs{viaTi-Li){^^^ ' ''^^^^ 

where via = Pia/Eia, with Eia = \/p\a + "T-a- The velocity V2i and the energy E2i are similarly 
defined. Once more 'wave packets' can be associated with the oscillating particles 1 and 2, but 
they are correlated by the finite width CTxeffs- The phase can be expanded around 2 
mi_2- For example, 

Spia + ^r-^ + ^r-^ Li , (206) 
2pi J 2pi 

where Spia = pia — Pi ■ The phase difference between two amplitudes vanishes but for the oscil- 
lation term already encountered in the previous sections, provided the wave packet is sufficiently 
localized, so that the term proportional to viTi — Li in (p06|) is negligible. 
At the extremum, the argument of the overlap function reads 

f[Pla,P2i) = ^ , (207) 

and yields a localization condition, that is, the mass differences cannot be too large compared 
with the energy-momentum widths. The coefficients a and 7 read 

1 1 {v2-VD,f {v2-VD2f\ 



^2' 



\'-''pD2 "^602 / y^^pp'^eDi ^eP^pD^ '^ep'^eDi J ^pD^'^eD^'^pef f2 

_ III {vi-Vd^){v2-Vd^) {v2-Vd^){vi-Vd^)\ {vi - V d^){v2 - V d^) 



+ - + 



^pp^ep ^eDj e^2 / ^ pP^ eDi'-' 

, _}_ ( _J_ , jvi - Vdi){'Vp - Vdi) \ ( _}_ {v2-Vd^){Vp -Vd^) \ 
0-2 i 0-2 C72 i I 0-2 cr2 I ' 

where 0"^^ = ^pi^ + f^'poj ■ The coefficient /3 is obtained from a by the exchange of the indices 
1 ^ 2. 

What happens when the energy-momentum uncertainty at the source goes to zero? Three 



features of the amplitude (|205D must be examined: the localization condition ( POTI) given by the 
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function f{piaiP2i), the value of the phase and the 'wave packet' effect on the time average. In 
the hmit app — > (so that aep — > too), the argument of the overlap function becomes at the 



extremum 



lim /pia,P2i = 77 ^2 + Tl \r~2 + "/!/ \2 2 

Thus a small, or even zero energy-momentum uncertainty at the source, does not destroy the 
oscillations. This result confirms our expectations, since a zero width at the source does not give 
information on which mass eigenstates propagate. On the other hand, the function f{pia^P2i) 
diverges when either Up^^ or CTpo^ Soes to zero, yielding the expected localization conditions at 
the detectors. 

In the limit of a zero energy-momentum width at the source, wave packets associated to 
oscillating particles are of infinite extent. This phenomenon was expected, since the time of 
production Tp, included in Ti and T2, becomes ill-defined. However the average on the time 
of production does not destroy the oscillations, at least if the detection times are measured. 



The explanation follows. In the limit Opp — > 0, Eqs. (202) and (203) show that the quantities 

^Via = Pia - Pd^ and 6p2i = P2i - Pd2 become 

lim 5pia = - lim 5p2i = _.^LlL^ . 

o-pP^O o-pP^O Vi — V2 

Note that these values can be directly obtained, in the one-dimensional case, from energy- 
momentum conservation at the source. The phase of the amplitude (see Eq. (|206[) ) can then be 
written as 

V (A ^A. \ X Li - L2 - V2{Ti - T2) , ^ Li - L2 - vijTi - T2) 

\\m [(pia + (p2i) = oi VO2 , (208) 

o-pP^O Vi — V2 Vi — V2 

where the term of order zero in the mass differences has been dropped. Since T1—T2 = tp,^ — tp,2 
and Li — L2 = Hb^ — yuj' the phase is independent of the production point. However the 
production point is implicit in the phase ( |208 ) when the classical limit is taken. The substitution 



of the classical relations Li^2 = 'i^i,27i,2 shows indeed that the expressions multiplying and 82 
are the propagation times Ti and T2, respectively. In other words the space-time coordinates of 
the production point can be reconstructed from the coordinates of detection and the velocities 
of the oscillating particles, with the exception of the vi = V2 case (recall that velocities of 
the same sign have the same direction). While complete computations in three dimensions are 
complicated, it can be seen that the phase difference becomes independent of the production 
coordinates, in the limit of zero energy-momentum uncertainty at production: 

lim {4)ia + (i)2i - 4>lb - (t>2j) = SEiab{Ti - T2) - (^Plafc • (Li - L2) . (209) 

(Tpp— >0 

10.4 No recoil oscillations 

It has been claimed that particles produced together with mixed states also oscillate because 
of the energy-momentum recoil. For example, the A baryon, in the process 7r~p IS 
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supposed to be in a superposition of two energy-momentum eigenstates, so that its detection 



probability should oscillate in space [244]. Furthermore the oscillation frequency of the kaon 
is modified by the momentum recoil against the A. In another example, muons produced in 
vr — > /iz^ are supposed to oscillate in space, because of the momentum recoil against the mixed 



state of the neutrino |24€, 247 1. As explained in section |2^ , recoil oscillations arise in quantum- 
mechanical models if different propagation times are associated to the different mass eigenstates. 
Recoil oscillations are unacceptable from basic principles, since the oscillation of the detection 
probability of a non- mixed state means that probability is not conserved at all distances. On 
the contrary, when several mass eigenstates are mixed, the sum of the detection probabilities of 
the different mass eigenstates is always equal to 1 for a given propagation distance. 

This assertion has been questioned in several papers in the framework of quantum mechan- 
ics. These refutations involve either the consideration of the different proper times associated 
with the oscillating particle and the recoil particle |201, 258, 30|, or the use of classical trajec- 
tories [107], or the use of energy-momentum conservation for the average energy-momenta of 



the propagating wave packets ]255, 212]. The source-propagator approach of Shtanov ]p42| is 
not a real improvement on the arguments of Ref. [107|, since Shtanov treats configuration space 
variables, which have a microscopic role, as classical macroscopic variables. 

Quantum field theory allows to understand which quantum-mechanical explanations are 
correct and important. The correlated oscillation model of the previous section is easily adapted 
to the treatment of recoil oscillations. Suppose that the first particle is in a superposition of 
several mass eigenstates mia, whereas the second particle has only one mass eigenstate m2 (it is 
the so-called recoil particle) . Since the phase of the amplitude ( p05D depends only on the average 
propagation times Ti^2 and distances Li,2) which are common to the different mass eigenstates, 
no recoil oscillations will be observed as long as these quantities are well-defined. This can be 
checked explicitly if app > Cp^^ ^, as the overlap function factorizes in that case. Thus the only 
case requiring a careful examination is when the uncertainty on the position of the source is 
larger than the oscillation length, i.e. when app < 5m\^f^/2pi. This condition is not satisfied 
in the process 7r~p KK^ . The momentum uncertainty of the proton, in this fixed-target 
experiment, is around 3 keV, i.e. much larger than the kaon mass difference of 3 x 10~^ keV 



[195]. Therefore the problem of recoil oscillations does not arise in the KK^ complex. 

In contradistinction to the case of the 7r~p collision, it is possible to consider a tt decay 
with a sharp momentum. Its theoretical minimal uncertainty is given by its 2.5 x 10~^ eV decay 
width, which is smaller than the neutrino mass differences. This case has been treated, with 
essentially correct quantum-mechanical arguments, by Dolgov, Morozov, Okun and Schepkin 
]107] and within a model coupling the ^ and v propagators with a source wave packet by 
Dolgov JlOS]. On the one hand, if only the neutrino is observed, the detailed study of the one- 
particle oscillation case in section |6| has shown that neutrino oscillations occur, provided that 
the momentum width associated with the muon is larger than the inverse oscillation length. On 
the other hand, if only the muon is observed, no oscillations in the probability of muon detection 
occur. The reason is that the sum over the neutrino flavors makes the transition probability 
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jl52| ) independent of the distance (apart from the geometrical decrease in L 



Furthermore oscillations also vanish after integration over the unknown source decay coordinates. 
This point was discussed in section |9.3| . 

Finally, the only case left is the detection of both muon and neutrino, coming from the decay 
of a pion endowed with a sharp momentum. As explained above, this process can be treated 
like a correlated oscillation, but with only one mass eigenstate m2 for the muon. The formula 
(|20g| ) shows that the oscillation of the probability is independent of the production point, so 
that the average over the production point does not have any effect on the phase. Moreover the 
oscillation term can be expressed as depending only on the detection point of the neutrino, and 
on the reconstructed decay point (see Eq. ( ^08|) with (52 = 0): 

T (A I A \ ^ L1-L2- V2{Ti - T2) 

hm [(pia + (p2) = di . 

o-pp-»0 Vi — V2 

The substitution of the classical relations Li^2 = ^^i,2?'i,2 in the above equation shows indeed 
that the phase difference depends only on the classical propagation time of the neutrino: 

lim (01, + 02) = = . 

In short the probability to detect both muon and neutrino can be interpreted as the result of 
a neutrino oscillation alone. Recall that the above relation can be obtained in a simple way in 
one dimension by energy-momentum conservation at the source. 

As regards the experimental data, oscillations of recoil particles have not been probed (in 
the case of the vr— >/iz^, it would only be possible for a neutrino oscillation length shorter than 
the muon decay length). However the neutral kaon mass difference obtained from strangeness 



oscillation experiments is consistent with the one obtained from regeneration experiments |201]. 



10.5 Oscillations of correlated mesons 

We compute here the amplitude associated to the correlated oscillations of a pair of neutral 
mesons produced in e"^e~ annihilation. The phenomenology of correlated K and B mesons is 



discussed in |T0|, |7|, |T5|, ^ and in |8|, g |lT|, |73l, respectively. 



Each meson oscillates between its components Pl — Ph, before decaying into final states 
fi{ki) and 72(^:2) at spacetime points yi and y2'- 

e+e- ^ n{q) ^ P^pO ^ h{ki)f2{k2) , 

where g, ki and k2 are the corresponding energy-momenta. 

Since each final state can be produced by either P^ or P^, the two amplitudes arising from 
the exchange of P^ and P^ as intermediate states must be coherently added. The resonance 
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quantum numbers of (f) and T(4s), J^*^ = 1 , are conserved by the strong interactions at the 
source, so that the pair 

pOpO 

is in an antisymmetric state under P and under C. Thus the 



relative sign of the two contributions to '7^— >/i/2 is negative |19(:] 



The total amplitude before antisymmetrization is given by Eq. (|20C|). In section 10.2, it was 



shown that the energy-momentum correlation at the source can be neglected for the resonances 



(j) and T(4s), so that the overlap function factorizes as in Eq. (201). The amplitude for correlated 
oscillations can thus be approximated by the product of two amplitudes, each corresponding to 
a single oscillation. Furthermore, it was shown in section |8.3| that the oscillation amplitude, in 
the case of nearly degenerate unstable particles, can be replaced by the effective amplitude ( |189| ) 
depending only on the distance L. In the end the antisymmetrized amplitude to detect a fi at 
distance Li and a /2 at distance L2 can be written as 

/ .,^mQL-i / 1 \ \ / ■?\/f"^0^2 / 

Tf,f, ~ Mp {Moi M-oi) V-'e-'''—V 1 (A^02 Mq^) V-^e-'^—y ' 



\ \ 1 



same expression with 








where A^P = M{n P°P^). The matrix M is defined by M = diag{mL - iTL/2,mH - iTH/2), 
niQ is the mass of the kaon in the degenerate limit and pj are the norms of the 3-momenta of 
the final states fj. The diagonalization matrix V is given by Eq. ( |192D . The amplitudes Mqj 
and Mqj stand for M{P^ ^ fj) and M{P^ ^ fj), respectively. 
The amplitude can be written as 

-'/1/2 ^ MpMoiMo2 



2a 



where Hf is defined by Eq. ( |194| ) and r/j = {l — a^f)/{l+afif). Note that the amplitude vanishes 
for identical final states (/i = /2 and pi = P2) and identical propagation distances (Li = L2), as 



expected [100|. In the center-of-mass frame, the equality of momenta Pi = P2 = P implies that 



the interference term oscillates like 

cos { p (Li - L2)j . (210) 

This oscillation formula coincides with the quantum-mechanical result obtained with the equal 
time prescription, whereas the different time prescription leads to an oscillation length shorter 



by a factor 2 |245]. As already explained in section 2.3, the quantum field treatment shows that 
neither prescription is meaningful, since the interference takes place over a space and time range, 
determined by the effective 'wave packet' width. We can only say that, in the end, the different 
mass eigenstates have the same average propagation time or distance, so that the oscillation 
formula will agree with the one obtained with the identical time prescription. As regards the 



experimental data, Kayser |178] has shown, in the case of the process T ^B^B^, that the mass 



difference obtained from the oscillation formula ( piOf ) is in agreement with the mass difference 
extracted from single B oscillations. 
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11 Summary and outlook 



Although the plane wave derivation of the vacuum oscillation formula has often been criticized 
in the literature, it is still used in most articles and textbooks. There are two reasons for this. 
First, the oscillation formula obtained in this way is believed to be correct by most physicists, in 
spite of the numerous inconsistencies present in its derivation. Second, the other approaches are 
not felt to be completely satisfying, with the consequence that more sophisticated treatments 
are not considered worth the effort. 

Let us first point out that we do not claim that the oscillation formula obtained with the 
plane wave approach should be revised. However we have argued that this approach becomes 
unacceptable when all its inconsistencies are added up: the perfect knowledge of the momentum 
precludes spatial oscillations, observability conditions (such as ax < L°^'^ or L < L™'^) are 
not taken into account, flavor states are ill-defined, unstable oscillating particles cannot be 
consistently described in that model. Furthermore, the plane wave derivation requires two 
prescriptions: the classical propagation condition — j;| <C t) and the equal time prescription 
{ti = tj). Whereas the first prescription can be justified in a quantum- mechanical wave packet 
approach, the second prescription can only be proved right in a field-theoretical treatment. This 
is an important argument in favor of the latter approach, since the diff'erent time prescription 
{ti 7^ tj) leads to oscillation formulas differing significantly from the standard result. We have 
also shown that an energy- momentum prescription (such as Ei = Ej or pi = pj ) is not necessary 
to derive the oscillation formula. 

Once the necessity of a more sophisticated approach is understood, we are faced with 
the choice between two main methods: on the one hand, a quantum-mechanical treatment 
associating wave packets with the propagating mass eigenstates and, on the other hand, a field- 
theoretical treatment where the oscillating particle is considered as an internal line of a Feynman 
diagram. The first possibility is rich in physical insights (such as the existence of a coherence 
length) but not satisfactory in many respects: flavor states are still ill-deflned, the equal time 
prescription is still needed, the nonrelativistic limit is problematic if the mass eigenstates are 
not nearly degenerate, the size of the wave packet is hard to estimate, the coherence length is 
not well-deflned in the nonrelativistic limit, the case of a vanishing energy uncertainty is not 
included, and flnally the treatment is inadequate for unstable particles. Whereas most of these 



problems are solved in the interacting wave packet model of Giunti [138|, this last model requires 
quantum held theory to compute the interactions and is as complex as the external wave packet 
model. 

Derivations of the oscillation formula resorting to fleld-theoretical methods are not very 
popular, although they seem the only way out. The flrst reason is that they are thought to be 
very complicated. The second reason is that the existing quantum held computations of the 
oscillation formula do not agree in all respects. The aim of our report was to counter both 
objections. The first of these is easily refuted by noting that the oscillation formula for a stable 
particle can be derived in a very simple field-theoretical model, the Kobzarev et al. model 1 187 ] , 



106 



in which the source and the detector are approximated by infinitely heavy nuclei. This is the 
simplest model of all in which the oscillation formula can be consistently derived. The second 
objection could only be countered by a detailed and complete computation of the oscillation 
formula with all approximations carefully considered. 

In this report, we have shown that all existing field-theoretical treatments (with the note- 
worthy exception of the Blasone-Vitiello approach [^]) can be included in the so-called external 
wave packet model, where the oscillating particle is described as an internal line of a Feynman 
diagram and propagates between a source and a detector localized with in- and outgoing wave 
packets. In particular, we have paid attention to the case of a vanishing energy uncertainty 
(stationary limit) and its wave packet interpretation, so that we could prove that the Grimus- 
Stockinger model [|150(| is a subcase of the external wave packet model. Two other limits have 
been considered. First, we have proved that oscillations near the source can be described by 
the standard oscillation formula, contrary to what was claimed by loannisian and Pilaftsis |163| ] . 
Second, we have also shown that oscillations far from a well-localized source are also described 
by the standard oscillation formula, contrary to Shtanov's claim [p42[| . In order to include all 
these cases in our formalism, it has been necessary to evaluate the amplitude with two approx- 
imation schemes: Laplace's and stationary phase methods. As a result, the propagation range 
was divided into three regimes, distinguished by the dispersion (or spreading in space-time) of 
the amplitude. In the first regime, near the source, the dispersion is negligible and Laplace's 
method is sufficient. This case corresponds to loannisian and Pilaftsis' limit. In the second 
regime, the transversal dispersion becomes important, so that both approximation methods 
need to be used. In the third regime, far from the source, the longitudinal dispersion becomes 
important and the stationary phase method is sufficient. This case corresponds to Shtanov's 
limit. We have also discussed another method of computation (in energy-momentum space) 
valid for stable oscillating particles, as it clears up in which sense interference occurs between 
same energy states. 

The computation of the intrinsic decoherence at the wave packet level has yielded two well- 
known observability conditions for oscillations. However, our physical explanation differs from 
the standard one by taking into account the dispersion of the amplitude, so that the origin of 
decoherence depends on the distance at which it occurs. The first condition of observability of 
oscillations, L°j^ » ax, is either due to the initial spread of the associated 'wave packets', or to 
their constant overlap at large distance. Similarly the existence of a coherence length (yielding 
the second observability condition L <C L'^"^) is either due to the separation of 'wave packets', or 
to their dispersion. Our 3-dimensional treatment has also yielded a third observability condition, 
Ipo 1| ^ <7p (modified into |vo x L| < v^ax at short distance) which has the obvious geometrical 
interpretation of constraining the propagation within a cone (respectively a cylinder at short 
distance). It is now clear that the coherence length tends to infinity in the stationary limit, as 
noted by Kiers, Nussinov and Weiss |182| ]. 

The oscillation of unstable particles, whether mesons or neutrinos, has been analyzed in the 
same framework, with the result that the oscillation formula is only modified, as expected, by 
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the well-known exponential decay term. Nonexponential corrections to the oscillation formula 
have been estimated for the first time in quantum field theory. Finally we have given a field- 
theoretical treatment of correlated oscillations, with the aim of disproving the existence of recoil 
oscillations. 

Our work confirms that the standard vacuum oscillation formula given by Eqs. and 



(14) is correct if the observability conditions mentioned above are satisfied. This fact is rather 
remarkable, in the light of the numerous inconsistencies which were present in the plane wave 
derivation of this formula. It can be explained by the following observations: 

• the energy and momentum present in the phase are linked by the on-shell condition, and 
the average propagation time and length are linked by a constraint equivalent to a wave 
packet localization in space-time. For these reasons, the EjT and pjL components of the 
phase cancel but for the usual oscillating factor ; 

• only same energy components interfere because of the time average, so that the time- 
averaged transition probability computed in the external wave packet model is equivalent 
to the energy- integrated transition probability computed with plane waves; 

• from the previous argument, decoherence effects arising from external wave packets cannot 
be distinguished from decoherence effects arising from averages over the production (and 
detection) region and over the energy spectrum; since the former effects are dominated by 
the latter, they have no influence on the oscillation formula. 

For ten years, new experimental results in neutrino physics have stimulated the research on the 
theory of particle oscillations, leading to the development of numerous (and sometimes conflict- 
ing) field-theoretical approaches. In this report, we have tried to unify these treatments and to 
extend them to oscillations of K and B mesons. We indeed believe that the neutrino oscillation 
formula is much strengthened by considering it in the same framework as the mesonic oscilla- 
tion formula, whose parameters have been confirmed by other methods such as regeneration in 
matter. We hope to have convinced the reader that field-theoretical models provide, on the one 
hand, the most secure foundation to the oscillation formula and give, on the other hand, a good 
physical understanding of oscillation phenomena, whatever the particles involved. 

One case needing further investigation is the decay in flight of the source of the oscillating 
particle, which is important for atmospheric neutrinos and some laboratory neutrino experi- 
ments. It would also be interesting to apply a field-theoretical treatment to cascade decays, 
such as B ^ J/ipK ^ J/ipfx, where double-flavor oscillations occur. These processes are useful 



for the determination of the signs of the BjjBl mass and lifetime differences [ p2| , 179, 77, p33[ |. 

The field-theoretical approach faces its biggest challenge in the description of neutrino 
oscillations in matter. A first task consists in the rederivation of the MSW oscillation formula 



(see [205, 214, 223, 87 1 for a few attempts). Another difficult task is to describe the nonstationary 
evolution of neutrinos in supernovae or in the early universe, where interaction rates are in 
competition with the flavor oscillation period |234, p^. 
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The turn of the century marks a new golden age for CP violation experiments, thanks to 
the B factories, as well as the golden era of neutrino oscillation experiments. These favorable 
auspices should be a great stimulation for further research on the theoretical foundations of 
oscillations. 
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Appendix: the Jacob-Sachs theorem 

Let 'ip{E,p) be a function which is distinct from zero only within certain bounds: 



with = E"^ — p^. On this interval ^(£^,p) is taken to be infinitely differentiable. Let the 
function I{T) be defined by 



il^iE, p) / for < Ml < / < M| and E>Q 




where the function GijP') is the complete scalar propagator in momentum space. 

We are going to prove that the function I{T) has the following asymptotic behavior: 



I{T) Z ^i V^^, p) e-^V^ 



(211) 



where zq is the pole of the integrand and Z is the residue. This section follows | 166 |. 



It will be easier to work with an overlap function depending explicitly on , so we write 
■0(p^,p). We choose to implement the assumption of compact domain of '(/'(p^iP) by 

p) = _ MfY (/ - M|)" p) for < < / < (n > 0) , 
V'(p^)P) = otherwise. 



The symmetrical behavior at Mf and M| is chosen only to simplify the algebra. The function 
r2 is taken to be analytic over the domain of interest. 



109 



Under the change of variable z = p , the function I{T) becomes 

I(T) = - / dz (z + p2)- 2 p) G{z) e-V-+p'r . 
2 Jul 

All functions in the integrand are analytically continued in the complex plane. For example, the 
real propagator is considered as the limit of an analytic function of the complex variable z: 

G{p^) = lim G{z) . 

The full propagator associated to a scalar field can be written in the Schwinger's spectral 
representation [241] as 



where Mq is the bare mass of the scalar field and n(z) is an analytic function defined by 

U{z) = z r ds , (213) 
Jyi z - s 

which can be interpreted as the self-energy of the scalar field. The positive real function (t{x) 
satisfies cr{x) = for x < 6^. is the invariant mass of the lightest multi-particle state in 
interaction with the scalar field corresponding to the propagator. This spectral representation 
has the advantage of being in direct correspondence with the full propagator (Eq. (p6|)), obtained 
in perturbation theory by an infinite sum over the self-energy insertions. If x is on the real axis, 
it can be checked that 

I\.{x + ie) — n(x — ie) = —2Trixa{x) . 

The self-energy can be written as 

n(x ± ie) = 5M'^ + u{x) =F iv{x) , (214) 

where v{x) = TTxa{x) and 6M'^ is an infinite constant which is absorbed in the renormalized 
mass M^: 

= + 6M'^ . 

Note that M is not the physical mass since it is not the pole of the propagator. 

The propagator G{z) has branch points at the thresholds corresponding to the multi-particle 
states. The first branch point is at z = b'^. We assume that the other branch points are above 
M|. The analytic continuation of G{z) into the second Riemann sheet |78] is such that the 
analytically continued function Gjj{z) just below the real axis is equal to the original function 
just above the real axis. The analytically continued function Hu^z) of the function n(z) is 
defined in the same way, Ilii{x — ie) = Il{x + ie), so that 

Gniz) = \ (215) 

z- -nii[z) 

Gii{z) has a pole at ZQ = m? — imV. 

Different cases have to be considered, according to the stability or instability of the particle, 
and according to the positions of the pole zq and threshold b^. The pole is assumed to be inside 
the contour (otherwise the particle cannot be emitted). 
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Figure 9: Integration contour if 6 < Mi < m < M2. 



Unstable particle, with h < Mi < m < M2 

The path of integration is shown on Fig. 9. The analytic expressions of the paths Vj read 

: z= (-iu; + ^M2 + p2) - p2 (j = l,2), 

where uj ranges from zero to Wqo on Vi and from a; 00 to zero on 1^2- The analytic expression of 
the path T" reads 

where M ranges from Mi to M2, with the limit lOoo 00 understood. 
The integral I{T) is equal to 

/(T) = J + J1 + J2 + J' . 

J is the contribution of the pole zq, whereas Ji, J2 and J' are the contributions of Vi, V2 and 
V' , respectively. Their analytic expressions read 

1. Contribution of the pole: 

Setting AM = \m — Mi 2I, a rough estimate of J is given by J ~ w?''^~^ (AM)^" e"'"'^/^. 

2. Contribution of the sickle-shaped paths: 

J,- =i(-iy 6"'^^^"+^'^ r duj'iP{z{Lo),p) Gn{z{uj)) e-"^^. 

Jo 

At large T, the dominant contribution to the integral comes from the uj values near 
zero, because of the decreasing exponential. Since the integrand tends to zero with to 
(■0(M?,p) = 0)), the asymptotic behavior of the integral depends on the way the inte- 
grand tends to zero. Setting y = ujT and expanding in 1/T, we obtain: 

J, ~ (Ml-Mf)" (M|+p2)"/' n (M|, p) Gn (M|) T-^^+D e" . 
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The corrections to this formula are of order (AMT)~^'^^'^\ where AM « Mj — m, that 
is, they are of the order of the uncertainty on the particle mass. It is assumed that the 
function ^}{z, p) diverges more slowly than exp(— ojT) on the paths Vj as uj tends to infinity. 
This assumption is true for a wide class of functions, in particular for Gaussians, whereas 
their path integral diverges on half-circles at infinity. A rough estimate of Jj is given by 
Jj ~ (AM)"-i r-("+^). 

3. Contribution from the path at infinity: 

J/ ^ ^-^^T / ^ {z{M)) Gh {z{M)) e- V*^'+p' t _ 

J Ml VM2 + p2 

If Q. satisfies the same conditions at infinity as above, J'~exp(— Woo?") and tends to zero 

as Woo — ^ co- 
in conclusion, the contribution of the pole is a decreasing exponential in T, whereas the con- 
tributions due to the bounded character of the energy spectrum decrease in inverse powers of 
T. 

At small T, Jj/J ~ (AMT)~("+^), which is nonneghgible for AM T < 1. Below that value, 
the asymptotic evaluation of Jj is not valid anymore, because terms in (AMT)"*^""*"^) have been 
neglected. 

First, let us consider weakly decaying particles. In the case of the X^, the mass is measured 
with a precision of AM ~ 10~^ MeV. Thus nonexponential corrections will be important for 
T < 10~^^ s, which is not observable since the lifetime is 0.89 x 10"^'' s. In the case of the 
, AM « 2 MeV, so that nonexponential corrections will be important for T < 10"^'^ s, which 
is not observable since the lifetime is 1.29 x 10~^'^s. Next, let us next consider resonances. 
In the example of A(1232), AM ~ 2 MeV, so that nonexponential corrections are important 
for T < 10^^^ s, which is large compared to the inverse width equal to 5 x 10^^^ s. Thus the 
propagation of resonances can never be modelized by the contribution of the pole alone. 

In the case of the weakly decaying particles, power law corrections are also important at 
large times and dominate if 

rr > 2(n + 1) In(AMr) = 2(n + 1) (In(rr) + ln(AM/r)) . 

For ET^, the two-particle threshold is characterized by n = 1/2 and AM/F ~ 0(10^'^), so that 
the threshold of the nonexponential behavior is given by FT > 69, at which time the amplitude 
will be much too small to be observable. For B^, taking n = 1/2 and AM/F ~ O(IO^), we 
obtain FT > 62, with the same conclusion as above. 

To sum up the analysis of the propagation of a weakly decaying particle (with no multi- 
particle thresholds included in the spectrum), the function I{T) is very well approximated in 
the observable time domain by Eq. ( pill ). 
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Figure 10: Integration contour if M\ < b < m < M2. 



Unstable particle, with Mi < b < m < M2 

In this case, the multi-particle threshold is included in the energy spectrum. It generates new 
power law corrections to the amplitude. 

The path of integration is shown on Fig. 10. The contributions of the paths Vj are the 
same as in the previous case, except that G{z) replaces Gii{z) on Vi, since this path lies on the 
first Riemann sheet. The asymptotic value of Jj as T— >-oo does not change. The contributions 
of the paths Vbj do not cancel each other out, since they do not lie on the same Riemann sheet. 
Their sum reads 



\ f dz {zWr-- ^{z, p) {Gu[z) - Giz)) e~^^^ 

_i^-iV^+^T r dLoi^{z{uj),p) {Gn{z)-G{z)) e-'^^. (216) 



The analytic expression of the paths Vb reads 



Vb-- z= + v'&^ + p2 -p\ (217) 



where uj ranges from zero to tOoo, with the limit Wqo ^ 00 understood. The function I{T) is 
equal to 

I{T) = J + Ji + J2 + Jb + J' ■ 

The asymptotic behavior of Jb is studied in the same way as for Jj. As before, the dominant 
contribution to Jb, for large T, comes from oj values near zero. Again, we set y = uoT and expand 
the integrand in 1/T. Recall that Gii{z) is defined by the analytic continuation of G{x + ie) 
below the cut, Gii{x — ie) = G{x + ie). Its discontinuity through the cut can be computed with 
Eqs. (I2T2I ), (HI) and ( ^ : 

n I ■ \ n( ■ \ ^^^^^ 
Gn\x — le) — G{x — le) 



{x — M'^ — u{x)) +v'^{x) 



113 



In the example of K^, the self-energy is given at first order in by the pion bubble diagram 
(where g is the coupling constant between a kaon and two pions), yielding 

where b = 2m7r is the two-pion threshold. The function u{z) is also of order g^. 
The evaluation of Jb with the same asymptotic method as used for Jj gives 

2 //,2 _|_ 2\^/'^ 

J, = _i (_2.vr) I A_ +P ) P) e-^V^T 

The corrections to this formula are of order (QT)-^/'^, where Q = M - 6 is the energy release 
on decay of the unstable particle. Noting that cutting rules give v{m?) = mV, the coupling 
constant g"^ can be replaced by its expression in function of F, m and Q. A rough estimate of 
Jh is then given by 

Q m 

At small r, the ratio Jh/ J is of the order of 

which is much smaller than 1 if QT ^ 1 and T /Q <^ 1. This is the case for weakly decaying 
particles. For example, T/Q ~ 10~^^ for Kg. It is not true for resonances. For example, 
T/Q ~ 0.8 for A(1232). Below the value QT ~ 1, the asymptotic computation of Jb is not valid 
anymore since corrections in (QT)^^^'^ have been neglected. 

At large T, the power law contribution Jb dominates the pole contribution if 

rr-31n(Fr) >51n(Q/F), 

so that the nonexponential time thresholds are given by F5T > 165 for Kg, T^T > 202 for K^ 
and rL,HT > 157 for (in the last case, we have taken Q = IMeV). Thus nonexponential 

effects are not observable at large times for weakly decaying particles. In the case of resonances, 
the contribution of Jb is always of the same order or larger than the pole contribution. Reso- 
nances do not propagate macroscopically, since their width is of the same order as the typical 
energy of the processes. 

To sum up the analysis of the propagation of a weakly decaying particle, including the 
multi-particle thresholds, the function I{T) is very well approximated in the observable time 
domain by the Eq. ( pill ). 

Stable particle, with Mi < m < M2 < b 

This case can be examined like the first one, except that the pole is real, so that the result is 
given by Eq. (|211| ), with F = 0. There is no exponential decrease of the amplitude. Power law 
corrections are negligible at large T. 
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stable particle, with Mi < m < b < M2 



This case can be examined like the second one, except that the pole is real. The result is again 
given by Eq. (211), with T = 0. The ratio Jb/J is of the order of Jb/J ~ (mT)~^/^(7^/Q^, where 



g is the coupling constant with the particles produced at the threshold. The contribution of Jb 
is too small to be observed at small T in the current experiments, and is negligible at large T. 
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